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>  to  model  circular  fuselages  and  item.  The  noncircular  fuselage  and  elliptic 
•tom  surfacoi  am  modeled  with  constant  source  panels.  Nonlinear  corrections 
am  aade  to  tbs  wing,  fuselage,  rack,  store  and  fuselage  inlet  Models  to 
simulate  shocks.  Aw  program  also  calculates  the  trajectory  of  the  store  as 
it  separates  fron  the  aircraft.  Ibis  report  describes  the  program,  presents 
instructions  for  preparing  input  for  the  program,  describee  the  output  from 
the  program,  and  presents  a  sample  case.  The  program  represents  an  extension 
of  an  earlier  program  restricted  to  circular  bodies  at  supersonic  speeds 
written  by  the  present  authors  and  described  in  AFFDL-TR- 76-41. 

This  volume  presents  the  instructions  for  preparing  input  for  eac&  of 
two  programs,  a  sample  case  for  each,  and  the  descriptions  of  the  output 
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This  report  "Prediction  of  Supersonic  Store  Separation  Charac¬ 
teristics  Including  Fuselage  and  Stores  of  Noncircular  Cross  Section," 
describes  a  combined  theoretical-experimental  program  directed 
toward  developing  a  computer  program  for  predicting  the  trajectory 
of  an  external  store  separated  from  an  aircraft  flying  at  super¬ 
sonic  speed.  It  represents  an  extension  of  previous  work  covered 
in  AFFDL-TR- 76-41  to  include  more  realistic  modeling  of  fuselage 
shapes  including  noncircular  cross  sections  and  ramp  type  engine 
air  inlets,  and  to  include  modeling  store  shapes  with  elliptic 
cross  section  with  multiple  sets  of  arbitrarily  oriented  fins. 

Volume  I,  "Theoretical  Methods  and  Comparisons  with  Experiment," 
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sons  with  experimental  data.  This  volume.  Volume  II,  "User's  Manual 
for  the  Computer  Program,"  presents  detailed  instructions  on  the 
use  of  the  computer  program  with  emphasis  on  preparation  of  input 
data  and  interpretation  of  output.  Volume  III,  "Appendices  A 
and  B,  Details  of  Program  I,"  provides  additional  descriptions  of 
the  individual  subroutines  and  program  variables  passed  between 
modules  in  the  first  of  two  programs.  Volume  IV,  "Appendices  C 
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the  individual  subroutines  and  program  variables  passed  between 
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PREDICTION  OF  SUPERSONIC  STORE  SEPARATION 
CHARACTERISTICS  INCLUDING  FUSELAGE  AND 
STORES  OF  NONCIRCULAR  CROSS  SECTION 


Volume  II.-  User's  Manual  for  the  Computer  Program 


1 .  INTRODUCTION 

The  purpose  of  this  volume  of  the  report  is  to  describe  and 
present  instructions  for  using  the  two  computer  programs  developed 
in  conjunction  with  the  theoretical  work  presented  in  Reference  1. 
This  work  represents  four  year's  effort  on  the  development  of  an 
analytical  method  for  predicting  the  trajectory  of  a  store 
separated  from  an  aircraft  flying  at  supersonic  speeds.  The 
present  work  includes  the  extension  of  the  work  of  References  2 
and  3  to  include  the  generality  of  noncircular  parent  aircraft 
configurations  and  an  elliptic  store  with  two  sets  of  arbitrary 
oriented  fins.  This  extends  the  present  work  to  match  the 
capabilities  of  the  subsonic  store  separation  prediction  method 
and  program  in  References  4  and  5  though  different  theoretical 
approaches  are  taken.  Preliminary  versions  of  that  method  and 
program  are  described  in  References  6,  7,  and  8.  Because  of  the 
program's  size,  the  program  has  been  split  into  two  parts  to 
reduce  core  requirements  and  avoid  unnecessary  repetition  of 
calculations. 

In  extending  the  work  of  References  2  and  3  to  noncircular 
body  shapes  every  attempt  has  been  made  to  maintain  the  circular 
shape  options  in  their  original  form  and  to  parallel  their 
development  with  new  calculations.  Wherever .possible  two 
separate  and  distinct  options  exist  for  each  of  the  methods  with 
both  producing  the  same  information. 

The  descriptions  of  the  user  information  describing  the  input 
and  execution  of  the  various  routines  is  split  into  three  volumes. 
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The  present  volume  contains  the  descriptions  of  the  capabilities 
of  the  programs,  the  descriptions  and  formats  of  the  input  param¬ 
eters,  the  descriptions  of  the  output  quantities,  and  examples  of 
input  and  output.  The  second  volume  contains  Appendices  A  and  B 
which  provide  the  detailed  descriptions  of  the  operations  and 
equations  implemented  in  each  of  the  subroutines  and  the 
definitions  of  all  parameters  passed  through  both  labeled  and 
blank  common  in  Program  I.  The  third  volume  similarly  contains 
Appendices  C  and  D  with  the  detailed  descriptions  of  subroutine 
operations  and  parameters  passed  through  common  in  Program  II. 

All  information  necessary  to  specify  the  input  and  interpret  the 
output  should  reside  in  the  present  volume. 

The  two  programs  as  described  in  this  report  have  been  run  on 
the  CYBER- 17 5  computer.  The  programs  should  run  with  very  little, 
if  any,  modification  on  other  computers  with  large  enough  memory 
capacity.  The  first  program  requires  approximately  250,000  octal 
words  of  core  storage  on  the  CYBER-175.  In  addition  to  INPUT  and 
OUTPUT  files  six  external  files  are  used  internally  by  the  first 
program.  Five  of  these  are  temporary  files  in  which  all  informa¬ 
tion  is  generated  and  used  internally  without  direction  by  the 
user.  The  last  file  TAPE12,  is  used  to  store  the  dataset 
generated  by  Program  I  to  be  passed  to  Program  II. 

Program  II  requires  approximately  372,000  octal  words  of  core 
storage  to  run  the  maximum  size  case  which  can  be  specified  by  the 
input  on  the  CYBER-175.  In  addition  to  INPUT  and  OUTPUT  files 
seven  external  files  are  used  internally  by  the  second  program. 

Six  of  these  are  temporary  files  used  internally  by  the  program  for 
temporary  storage  of  data  only  during  the  duration  of  the  run. 

The  seventh,  TAPE12,  is  an  input  file  containing  the  data  generated 
in  Program  I  to  be  used  in  Program  II.  The  information  on  TAPE12 
is  written  unformatted  and  may  be  stored  on  any  external  device 
specified  by  the  user.  The  remaining  input  to  both  programs  is 
from  cards  and  all  output  is  printed. 
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The  next  sections  of  this  report  will  describe  in  more  detail 
the  general  use  of  the  program  and  then  present  instructions  for 
preparing  the  input  data  and  interpreting  the  output.  A  sample 
case  is  presented.  Details  of  the  program  are  contained  in  the 
appendices  in  Volumes  III  and  IV. 

The  method  of  describing  the  aircraft  components,  the  aircraft 
operational  parameters,  and  the  dynamical  characteristics  of  the 
store  are  of  importance.  The  various  parameters  which  are 
included  in  the  computer  program  are: 

Wing  Panels 

Thickness  distribution:  Specified  at  a  large  number  of 
chordwise  and  spanwise  locations  independent  of  the  constant 
u-velocity  panel  layout  used  to  represent  the  wing  as  a  lifting 
surface. 

Mean  camber  surface:  May  have  both  twist  and  camber. 

Leading-edge  shape:  Represented  by  straight  line  segments 
of  differing  sweep. 

Trailing-edge  shape:  Also  represented  by  straight  line 
segments  of  differing  sweep. 

Dihedral:  Represented  by  straight  spanwise  segments  of 
varying  dihedral. 

Fuselage 

Two  optional  shapes: 

1)  Body  of  revolution  represented  by  line  sources 
and  doublets 

i 
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2)  Noncircular  general  shaped  body  represented  by 
source  panels 

Inlet  shape:  Ramp  inlet  using  oblique  shock  external 
compression. 

Pylon 

Thickness  distribution:  Same  method  of  description  as  for 
wing  panel. 

Mean  camber  surface:  Planar. 

Orientation:  Vertical  and  streamwise. 

Leading-edge  shape:  Represented  by  line  segments  of  differing 
sweep. 

Trail ing-edge  shape:  Also  represented  by  straight  line 
segments  of  differing  sweep. 

Tip:  Parallel  to  pylon  root  chord. 

Number:  One  pylon  per  wing  panel  or  one  under  fuselage 
centerline. 

Rack 

Shape:  Body  of  revolution. 

Number  of  racks:  One  on  pylon. 


* 
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Store 
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Two  optional  shapes: 

1)  Body  of  revolution  represented  by  line  sources 
and  doublets 

2)  Elliptic  body  represented  by  source  panels. 

Number  of  stores:  Seven  of  any  combination  of  shapes. 

Number  of  shapes:  Maximum  of  seven  bodies  of  revolution 
or  maximum  of  two  elliptic  shapes. 

Attached  orientation:  Initial  pitch  angle  arbitrary, 
initial  yaw  angle  zero,  and  initial  roll  orientation  arbitrary. 

Empennages: 

1)  Body  of  revolution  may  have  one  set  of  planar  or 
cruciform  fins 

2)  Elliptic  body  may  have  one  or  two  sets  of  fins; 
each  set  may  have  one  through  four  fins  at  arbitrary 
orientation. 

Airplane  Operating  Characteristics 

Flight  path:  Straight  but  not  necessarily  horizontal. 

Flight  velocity:  Constant. 

Density:  Constant. 

Angle  of  attack:  Constant. 

Yaw  angle:  Zero. 
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Store  Inertial  Characte  ristics 


Moments  of  inertia:  Constant. 

Products  of  inertia:  Constant. 

Center  of  gravity  position:  Not  necessarily  on  store 
longitudinal  axis. 

Store  Ejection  Conditions 

Initial  translational  velocities:  Arbitrary. 

Initial  pitching  velocity:  Arbitrary. 

Initial  yawing  velocity:  Arbitrary. 

Initial  rolling  velocity:  Arbitrary. 

Power:  Option  of  specifying  a  thrust  time  history. 

Ejection:  Option  of  specifying  an  ejector  force  versus  time 
or  displacement  history. 

Two  analytical  methods  may  optionally  be  used  to  compute  flow 
fields  around  bodies  in  these  programs.  For  purposes  of  simplic¬ 
ity  of  nomenclature  the  technique  employing  line  singularities 
will  be  referred  to  as  the  circular  body  method.  Similarly,  the 
use  of  the  source  paneling  technique  for  fuselages  will  be  referred 
to  as  the  noncircular  fuselage  method.  The  use  of  source  paneling 
techniques  for  store  configurations  will  be  referred  to  as  the 
elliptical  store  method.  In  describing  the  input  to  the  non¬ 
circular  and  elliptic  body  methods  the  same  variables  and  input 
routines  have  been  used  for  both  fuselage  and  stores.  All  data 
required  for  further  calculation  are  saved  on  external  files  via 
internal  program  logic. 


6 


2.  GENERAL  DESCRIPTION  OP  THE  USE  OF  THE  PROGRAMS 


The  store  separation  trajectory  computation  at  supersonic 
speeds  requires  the  execution  of  two  programs.  The  first  program, 
LDCALC,  consists  of  three  sections.  The  first  reads  the  input 
and  control  parameters  required  to  specify  the  parent  aircraft 
and  store  geometries.  The  next  solves  for  the  singularity  dis¬ 
tributions  which  represent  the  wing-pylon  loading  including 
interference  on  the  fuselage,  the  circular  or  noncircular  fuselage, 
the  rack,  and  the  circular  or  elliptic  stores  excluding  the 
empennages.  The  third  section  takes  the  singularity  and  geometric 
data  required  for  continuing  the  calculations  and  saves  it  on  an 
external  file. 

The  second  program,  TRJTRY,  consists  of  tliree  sections.  The 
first  reads  the  file  from  Program  I  and  the  additional  data  to  be 
used  to  specify  the  separating  store  and  its  empennages  if  any. 

The  second  section  generates  the  influence  coefficient  matrices  for 
any  empennages  that  are  present.  The  third  section  calculates  the 
forces  and  moments  on  the  separating  store  and  computes  the 
trajectory.  This  division  into  two  programs  allows  the  user  to 
generate  the  parent  aircraft  geometry  and  singularity  strengths  once 
and  then  compute  several  trajectories. 

2.1  Descriptions  of  the  Use  of  Program  I 

The  first  program  is  organized  in  three  sections  to  read 
parent  aircraft  and  store  data,  to  generate  individual  body 
and  wing-pylon  singularity  strengths,  and  to  save  required 
results  on  an  external  dataset.  A  brief  description  of  the 
information  to  be  input  or  tasks  to  be  performed  in  each  of 
the  sections  follows. 
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2.1.1  Input  Section 


The  input  section  of  the  program  reads  in  the  following 
information. 

(1)  Parent  aircraft  flight  conditions 

(2)  Indices  specifying  what  aircraft  components  are  present 

(3)  Fuselage  data 

(4)  Wing  data 

(5)  Pylon  data 

(6)  Rack  data 

(7)  Store  data 

The  aircraft  flight  conditions,  item  (1)  above,  which  can  be 
varied  are  the  angle  of  attack  and  Mach  number  of  the  parent 
configuration. 

Item  (2)  consists  of  four  indices.  The  first  specifies 
whether  a  fuselage  is  present  and,  if  so,  what  type.  The  second 
specifies  whether  a  pylon  is  present  and  the  third  index  specifies 
whether  a  rack  is  present  on  the  pylon.  The  fourth  index 
specifies  how  many  stores  are  present.  The  program  in  its  present 
form  will  handle  up  to  seven  stores. 

The  fuselage  input  data,  item  (3)  above,  is  different  for 
each  of  the  fuselage  modeling  methods.  The  circular  fuselage 
method  requires  the  length,  maximum  radius,  and  polynomials 
specifying  the  body  radius  as  a  function  of  axial  location.  Other 
input  data  specify  the  number  of  sources  and  doublets  to  be  used 
in  representing  the  fuselage  volume  and  angle  of  attack  effects, 
the  fuselage  interference  panel  layout,  and  locate  the  nonlinear 
shock  wave  from  the  nose.  The  noncircular  fuselage  method  requires 
the  elliptic  axes  or  cross-sectional  shape  distribution  as  a  func¬ 
tion  of  axial  location.  Other  noncircular  fuselage  data  specify 
the  number  and  locations  of  panels  to  be  used,  the  interference 
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shell  layout,  and  generate  the  nonlinear  shock  shape.  The  panels 
used  to  represent  the  engine  inlets,  if  present,  are  indicated. 

The  wing  input  data,  item  (4) ,  locate  the  wing  relative  to 
the  fuselage  and  supply  information  required  to  lay  out  separate 
u-velocity  panel  and  thickness  panel  distributions.  These  data 
include  the  twist  and  camber  distribution  and  the  slope  distribu¬ 
tion  of  the  thickness  envelope.  The  wing  leading  and  trailing 
edges  can  have  breaks  in  sweep  and  dihedral. 

Similarly,  the  pylon  input  data,  item  (5) ,  locate  the  pylon 
and  provide  information  required  to  lay  out  the  u-velocity  panel 
and  thickness  panel  distributions.  The  pylon  is  located  laterally 
relative  to  the  fuselage  centerline  and  longitudinally  relative 
to  the  leading  edge  of  the  local  wing  chord.  The  data  also 
include  the  slope  distribution  of  the  thickness  envelope.  The 
pylon  cannot  have  twist  or  camber.  The  leading  and  trailing 
edges  can  have  breaks  in  sweep. 

The  rack  input  data,  item  (6) ,  are  the  length,  maximum 
radius,  and  location  relative  to  local  wing  chord.  Rack  incidence 
relative  to  the  wing  root  chord  is  also  specified.  Also  required 
are  the  polynomials  describing  the  radius  distribution  as  a  func¬ 
tion  of  axial  station.  Only  the  line  singularity  method,  circular 
body  method,  is  available  to  represent  rack  volume  and  angle  of 
attack  effects. 

Store  data,  item  (7),  are  input  which  assign  a  store  number 
and  specify  a  shape  number,  length,  and  maximum  radius.  Each 
store  is  located  by  specifying  its  lateral  position  relative  to 
the  fuselage  centerline  and  the  longitudinal  and  vertical  position 
of  the  store  nose  relative  to  the  wing  chord  immediately  above 
the  store.  The  incidence  of  the  store  relative  to  the  wing  root 
chord  is  also  specified.  The  store  must  be  under  the  fuselage 
centerline  or  to  the  left  of  this  line  as  seen  by  the  pilot. 


9 


Shape  numbers  less  than  or  equal  to  50  are  used  to  denote 
circular  store  input,  while  shape  numbers  greater  than  50  are 
used  to  denote  elliptical  store  input.  For  the  circular  store, 
data  used  to  determine  the  source  and  doublet  distributions  to 
represent  the  store  volume  and  angle  of  attack  effects  are 
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input.  These  data  include  the  number  of  singularities  and  poly¬ 
nomials  specifying  the  store  shape.  For  the  elliptic  store,  data 
used  to  determine  the  body  source  panel  layout  are  input.  These 
data  include  the  number  and  locations  of  the  panels  covering  the 
body  surface.  The  data  used  to  generate  the  nonlinear  shock  shape 
associated  with  the  store  nose  are  also  input. 

The  strengths  of  the  singularities  and  the  nonlinear  shock 
shapes  associated  with  the  presence  of  the  fuselage,  rack  and 
stores  are  computed  at  the  time  the  data  for  that  particular  body 
are  read.  Only  the  geometry  variables,  singularity  strengths, 
and  shock  shapes  are  saved  for  use  in  Program  II. 

2.1.2  Wing-Pylon  Loading  Section 

The  second  main  section  of  the  program  calculates  the  singu¬ 
larity  distribution  representing  the  wing-pylon  loading  including 
interference  on  the  fuselage.  This  is  done  using  the  method 
described  in  Section  3.3  of  Reference  1.  The  right  hand  side  of 
the  matrix  equations  used  to  determine  the  singularity  strengths  is 
generated  first  in  order  to  take  advantage  of  arrays  already 
present  in  blank  common.  In  determining  the  influence  of  the  air¬ 
craft  components  at  u-velocity  panel  control  points  for  purposes  of 
the  right  hand  side  computation,  only  the  fuselage  and  wing  and 
pylon  thickness  effects  are  included.  Their  influences  based  on 
linear  theory  are  considered,  that  is,  there  is  no  modification  due 
to  the  presence  of  nonlinear  shocks.  No  influence  due  to  the  rack 
or  stores  are  considered  in  determining  the  wing-pylon  boundary 
conditions. 
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The  coefficient  matrix  is  computed  next.  The  generation  of 
the  matrix  of  coefficients  overwrites  the  existing  data  in  blank 
common.  The  solution  of  the  set  of  simultaneous  equations  is 
computed  by  triangularization  of  the  matrix.  The  results  are  the 
strengths  of  the  constant  u-velocity  panels  divided  by  tt. 

Condensation  of  panel  strengths  associated  with  individual 
panels  to  net  strengths  computed  for  each  corner  is  done  last. 

The  net  corner  strengs  are  used  in  Program  II  to  compute  wing- 
pylon  influence.  For  panel  layouts  where  multiple  panels  lay  in 
the  same  plane,  up  to  a  factor  of  four  reduction  in  computational 
effort  is  achieved. 

2.1.3  Dataset  Generation  Section 

The  third  and  last  main  section  in  Program  I  saves  the  data 
required  for  the  calculation  of  the  store  trajectory  in  the 
second  program.  The  geometry  and  singularity  strengths  necessary 
to  compute  the  influence  on  the  separated  store  of  the  fuselage, 
rack  and  store  bodies,  of  the  wing-pylon-interference  shell 
u-velocity  panels,  and  of  the  wing-pylon  thickness  panels  are 
written  on  TAPE12.  The  data  to  be  transferred  for  the  noncircular 
fuselage  or  elliptic  stores  is  transferred  from  TAPEll  where  it  is 
temporarily  saved.  The  entire  dataset  on  TAPE12  must  be  saved  by 
the  user  for  use  in  Program  II  where  it  may  be  used  for  a  series 
of  trajectory  calculations.  The  dataset  is  for  a  particular  air¬ 
craft  configuration,  Mach  number,  and  angle  of  attack.  In 
Program  II  calculations  can  be  made  for  various  altitudes  and 
flight  path  angles.  If  multiple  stores  are  present,  the 
trajectories  of  the  various  stores  can  be  studied.  The  fin 
arrangement  on  the  separating  store  can  also  be  varied  as  can  the 
store  mass  and  inertia  characteristics.  If  the  thrust  and/or 
ejector  force  options  are  used,  these  two  force  histories  nay  also 
be  varied. 


11 


2.2  Descriptions  of  Use  of  Program  II 

The  second  program  is  organized  in  three  sections  to  read 
the  data  file,  to  input  additional  data  for  the  separating 
store,  and  compute  the  separating  store  forces  and  moments  and 
calculate  the  resulting  trajectory. 

2.2.1  Input  Separating  Store  Section 

The  first  section  inputs  the  information  necessary  to  define 
the  properties  of  the  separating  store.  This  includes  reading 
the  dataset  generated  in  the  first  program  containing  the  geometry 
and  singularity  strengths  modeling  the  parent  aircraft  components. 
The  input  section  consists  of  the  following  steps: 

(1)  Input  additional  titles. 

(2)  Read  parent  aircraft  and  stores  file  on  TAPE12 ;  if 
either  noncircular  fuselage  or  elliptic  store  data  is 
included,  the  fuselage  data  is  written  to  TAPEll  and 
the  elliptic  store  data  is  written  to  TAPE10  for 
future  use. 

(3)  Read  separating  store  flight  conditions. 

(4)  Read  additional  information  to  describe  separating 
store. 

The  aircraft  flight  conditions,  item  (3)  above,  which  can  be 
varied  are  the  angle  of  attack,  flight  path  angle,  Mach  number, 
free-stream  air  density,  and  flight  velocity.  The  angle  of 
attack  and  Mach  number  must  be  identical  to  those  used  to 
compute  parent  aircraft  results  read  in  item  (2) .  This  is  checked 
in  the  program  and,  if  it  is  not  true,  error  messages  are  written 
and  the  program  stops.  The  remaining  flight  properties  may  be 
varied  independently  in  Program  II. 


The  additional  data  describing  the  store  to  be  separated 
include  indices  specifying  the  store  number,  the  number  of  seg¬ 
ments  into  which  a  circular  body  is  to  be  broken  for  the  force 
calculation  and  the  number  of  sources  and  doublets,  whether  the 
store  has  empennages,  whether  it  is  powered,  and  whether  an 
ejector  is  present.  Also,  the  store  mass  and  inertia  character¬ 
istics  are  read  in  along  with  the  location  of  the  point  about 
which  the  aerodynamic  moments  are  to  be  calculated.  This  is  the 
point  about  which  the  moments  and  products  of  inertia  are  also 
calculated.  The  location  of  the  store  center  of  mass  relative  to 
this  point  is  also  specified  as  is  the  store  axial-force  coef¬ 
ficient  to  be  used  in  the  trajectory  calculation. 

Three  other  indices  are  input  which  pertain  to  options 
included  in  the  computer  program.  Provision  has  been  made  to 
include  or  exclude  the  damping  terms  in  the  velocity  field  calcula¬ 
tions.  Also,  for  a  store  with  an  empennage,  rolling  moment  may  or 
may  not  be  included  in  the  acceleration  determination.  The  third 
option  pertains  to  the  calculation  of  free-flight  trajectories 
as  opposed  to  captive-store  trajectories  as  obtained  in  the  wind 
tunnel.  In  wind-tunnel  captive-store  testing  it  is  customary  to 
change  the  store  pitch  and  yaw  angles  to  account  for  transla¬ 
tional  motion  only  while  measuring  the  aerodynamic  force  and 
moments.  This  changes  its  position  in  the  nonuniform  flow  field. 
Provision  has  been  made  in  the  computer  program  to  simulate  this. 

If  an  elliptic  store  or  noncircular  fuselage  is  included, 
the  data  for  an  elliptic  separated  store,  noncircular  fuselage, 
and  one  additional  elliptic  store  shape  is  arranged  sequentially 
in  blank  common.  Up  to  seven  stores  of  elliptic  shape  may  exist, 
but  there  may  be  only  two  independent  shapes. 

2.2.2  Empennage  Input  and  Loading  Section 

For  a  store  with  an  empennage,  additional  quantities  must  be 
specified  in  this  section.  Separate  analytical  models  for  the 
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empennages  are  used  with  each  of  the  circular  and  elliptical 
store  methods.  For  the  circular  store  method  a  single  set  of 
planar  or  cruciform  fins  is  modeled  using  slender  body  theory 
with  reverse  flow  theorems.  The  additional  data  are  an  index 
indicating  whether  the  empennage  is  planar  or  cruciform,  the 
tail-fin  semispan,  the  average  body  radius  in  the  tail-fin 
region,  the  initial  roll  orientation  of  the  fins,  and  the  lift- 
curve  slope  of  the  fins  alone.  In  addition,  the  axial  position  at 
which  the  tail  forces  are  assumed  to  act  is  specified. 

When  the  elliptic  store  method  is  used,  two  empennages  may 
be  modeled  using  three  dimensional  u-velocity  panels.  Each 
empennage  consists  of  from  one  to  four  fins  arbitrarily  oriented 
around  the  body.  The  additional  data  required  to  lay  out  the 
panels  on  the  fins  and  interference  shell  include  the  number  of 
spanwise  and  chordwise  rows  of  panels  on  each  fin  and  the  inter¬ 
ference  shell;  the  root  chord,  the  leading  and  trailing  edge 
sweep  angles,  and  exposed  semispan  for  each  fin;  and  the  meridional 
angle  and  dihedral  angle  for  attachment  of  the  fin  to  the  body. 

When  breaks  of  sweep  are  desired,  the  additional  data  for  each 
segment  must  be  supplied. 

2.2.3  Trajectory  Section 

The  third  and  last  main  section  of  the  trajectory  program  is 
the  trajectory  calculation  which  consists  of  the  following  steps: 

(1)  Input  thrust  and  ejector  force  time  histories  if 
powered  or  ejected. 

(2)  Initialize  for  trajectory  calculation. 

(3)  Calculate  aerodynamic  forces  and  moments. 

(4)  Calculate  accelerations  and  rates  of  change  of 
orientation  angles. 


14 


(5)  Integrate  equations  of  motion. 

(6)  Repeat  steps  (3),  (4) ,  and  (5)  to  end  of  trajectory. 

For  a  powered  store,  a  series  of  polynomials  is  used  to 
specify  the  thrust  history.  If  a  forced  separation  is  selected, 
a  similar  input  is  prescribed  for  the  ejector  force. 

In  the  trajectory  initialization  certain  store  separation 
conditions  are  specified.  These  are  the  initial  translational 
velocities  and  rotational  velocities. 

The  integration  of  the  equations  of  motion  is  done  by  a 
standard  numerical  integration  technique  with  the  aerodynamic 
forces  and  moments  calculated  at  each  point  required  by  the 
integration  scheme.  The  calculation  of  the  nonuniform  velocity 
field  and  the  resulting  forces  and  moments  is  described  in 
Sections  5  and  6  of  Reference  1. 

For  a  given  aircraft-store  combination  and  Mach  number,  it 
can  be  seen  that  a  series  of  trajectories  can  be  run  with  only 
minor  changes  to  the  input  data  deck.  For  example,  the  aircraft 
flight  path  angle  can  be  varied  by  changing  one  number  on  one 
card.  The  altitude  can  be  varied  by  changing  the  free-stream 
density  and  possibly  the  free-stream  velocity  to  account  for 
changes  in  the  speed  of  sound.  Among  other  things  easily  varied 
are  the  store  mass  and  inertia  properties,  center  of  gravity 
location,  and  ejection  conditions. 

Provision  has  also  been  made  for  restarting  a  trajectory. 
This  is  accomplished  by  changing  one  card  which  specifies  the 
initial  and  final  times  and  adding  two  cards  specifying  the 
current  values  of  the  dependent  variables.  These  are  tabulated 
in  the  output  at  each  integration  step. 


3.  DETAILS  OF  THE  USE  OF  PROGRAM  I 

The  first  program  used  to  generate  the  parent  aircraft 
geometry  and  singularity  strengths  consists  of  a  main  program 
and  59  subroutines.  Appendix  A  of  Volume  III  presents  a 
detailed  description  of  the  calculations  in  each  of  these 
routines.  Table  A-l  in  Appendix  A  lists  these  subroutines  in 
alphabetical  order  and  gives  a  one-sentence  description  of  what 
each  subroutine  does.  A  listing  of  the  routines  in  Program  I 
is  presented  in  Figure  A-l.  A  general  flow  chart  of  the  main 
program  for  Program  I  (LDCALC)  is  presented  in  Figure  1  of 
this  volume.  The  program  as  written  in  Figure  A-l  adheres  to 
ANSI  FORTRAN  standards.  Only  the  first  card  is  specifically  for 
the  CYBER-175  series  machines.  The  program  should  run  with  only 
its  removal  on  other  machines.  The  following  sections  outline 
the  flow  of  calculations  as  presented  in  the  flow  chart  of  the  main 
program,  LDCALC,  the  preparation  of  input  cards,  the  descriptions 
of  the  output,  and  any  special  conditions  and  messages  encountered 
within  the  program. 

3.1  General  Flow  Chart  for  Program  I 

The  purpose  of  the  first  program  is  to  read  the  geometric 
input  necessary  to  define  the  parent  aircraft  components  and  store 
bodies.  It  computes  the  strengths  of  the  singularities  associated 
with  volume  and  angle  of  attack  effects  for  each  of  the  components. 

A  general  flow  chart  of  the  first  program  is  presented  in 
Figure  1.  Constants  are  defined  and  heading  information  is  read 
and  printed.  The  aircraft  flight  conditions  are  input  as  are 
indices  specifying  what  aircraft  components  are  present.  If  a 
circular  fuselage  is  present,  the  fuselage  data  are  read  and  printed 
in  subroutine  FUSEIO.  This  routine  also  calls  BDYGEN  and  SHKSHP  to 
calculate  the  source  and  doublet  distributions  which  represent  the 
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fuselage  volume  and  angle  of  attack  effects  and  generate  the 
nonlinear  nose  shock  shape.  If  a  noncircular  fuselage  is  present, 
the  fuselage  data  is  read  and  printed  in  subroutine  WDYBDY.  This 
routine  also  lays  out  the  panel  geometry  in  GEOM,  solves  for 
their  strengts  in  VELCMP  and  SOLVE,  and  generates  the  nonlinear 
nose  shock  shape  in  BSHOCK. 

The  next  steps  in  the  program  read  in  the  data  required  to 
model  the  wing.  The  data  locating  the  wing  are  first  read  and 
then  subroutine  SWNGIN  is  called.  This  routine  reads  in  the  data 
required  to  lay  out  the  constant  u-velocity  panels  which  will 
represent  the  loaded  wing.  In  addition  it  reads  in  the  twist  and 
camber  distribution  at  the  panel  control  points.  Subroutine 
WLYOUT  is  called  to  lay  out  the  panels.  The  last  wing  input  data 
is  the  thickness  distribution.  These  data  are  read  in  by  sub¬ 
routine  WITHIN. 

A  check  is  made  in  the  program  to  determine  whether  or  not  a 
pylon  is  present.  If  one  is,  subroutine  PLYOUT  is  called  to  read 
in  the  data  required  to  lay  out  the  constant  u-velocity  panels 
which  will  represent  the  pylon  loading  This  routine  also  lays 
out  the  panels.  Subroutine  PYTHIN  is  next  called  to  read  in  the 
pylon  thickness  data. 

The  next  two  steps  in  the  program  are  calls  to  subroutines 
THKOUT  and  THKLYT .  These  two  routines,  respectively,  print  the 
input  data  for  the  Wing  and  pylon  thickness  and  lay  out  the 
thickness  panels. 

A  check  is  then  made  to  see  if  a  fuselage  is  present.  If 
the  circular  fuselage  is  present,  subroutine  BLYOUT  is  called  to 
lay  out  the  constant  u-velocity  panels  on  the  fuselage  interference 
shell.  If  the  noncircular  fuselage  is  present,  subroutine  BLY0T2 
is  called  to  lay  out  the  constant  u-velocity  panels  on  the  fuselage 
interference  shell. 
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If  a  rack  is  present,  subroutine  RACKIO  is  called  to  read 
the  rack  input,  compute  its  volume  and  angle  of  attack  effects, 
and  nonlinear  nose  shock  shape.  The  rack  is  then  located  in 
the  fuselage  coordinate  system. 

Provision  has  been  made  for  not  including  a  store  in  the 
input  data.  This  has  been  done  so  that  the  program  can  be  used  to 
determine  the  coordinates  of  the  points  at  which  the  wing  twist 
and  camber  distribution  must  be  input.  These  points  are  the  con¬ 
trol  points  of  the  wing  constant  u-velocity  panels.  Their 
locations  in  the  wing  coordinate  system  are  calculated  and  output 
by  the  computer  program.  Some  computer  time  can  be  saved  by 
making  a  run  without  any  stores  present.  The  next  step  in  the 
program  checks  to  see  if  there  are  stores.  If  there  are  none, 
the  store  input  section  of  the  program  is  bypassed.  If  there  is 
at  least  one  store,  subroutine  STORIO  on  the  second  page  of 
Figure  1  is  called  to  read  in  and  print  all  of  the  store  data. 

The  type  of  data  used  to  model  the  store  is  determined  from  the 
store  shape  number  read.  If  the  shape  number,  NSHAPE ,  is  less 
than  51,  the  circular  body  method  is  used  and  the  data  consist 
of  the  location  and  incidence  of  the  store  and  polynomials 
specifying  its  shape.  For  this  shape  subroutines  BDYGEN  and  SHKSHP 
are  called  to  compute  the  source  and  doublet  distributions  which 
represent  the  store  volume  and  angle  of  attack  effects  and  the 
nonlinear  nose  shock  shape.  If  the  shape  number  is  greater  than 
50,  the  elliptical  store  method  is  used  and  the  data  consist  of 
elliptic  semi-axes  versus  axial  station  and  the  additional  control 
integers  and  variables  to  lay  out  the  source  panels  on  the  store 
surface.  Subroutine  GEOM  is  called  only  once  for  each  different 
elliptic  shape  read.  The  source  panel  strengtns  are  computed 
for  all  elliptic  stores.  The  nonlinear  nose  shock  shape  is  computed 
only  at  zero  degrees  angle  of  attack  for  each  different  store 
shape.  Following  this  all  the  stores  are  located  in  the  fuselage 
coordinate  system. 
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The  next  five  boxes  of  Figure  1(b)  are  associated  with  the 
calculation  of  the  strengths  of  the  constant  u-velocity  panels  to 
represent  the  wing-pylon  loading  including  interference  effects 
on  the  fuselage.  The  equations  used  to  solve  for  the  strengths 
are  given  in  Section  3.3  of  Reference  1.  The  first  step  is  to 
call  subroutine  DPRHS  which  calculates  the  right-hand  side  of  this 
set  of  simultaneous  equations.  The  next  step  is  to  call  sub¬ 
routine  DPCOEF  which  calculates  the  coefficient  matrix,  that  is, 
the  coefficients  multiplying  the  unknown  panel  strengths  divided 
by  it,  (u+/V0O)/ir.  The  resulting  set  of  equations  is  then  solved 
by  calling  subroutine  INVER1.  The  next  section  of  the  program 
prints  the  resulting  panel  strengths,  control  point  coordinates, 
and  interference  velocities  at  these  points.  Since  up  to  four 
panels  may  share  the  same  corner,  subroutine  NULYT  is  called  to 
sum  the  net  strengths  at  panel  corners.  This  increases  the  cal- 
culative  efficiency  of  the  u-velocity  panel  influences  in  the 
second  program.  The  same  summing  up  is  also  done  for  the  thickness 
panels. 

The  last  step  in  the  first  program  is  to  save  the  data  file 
required  to  run  Program  II.  The  call  to  subroutine  WRFILE  writes 
all  parent  aircraft  geometry,  singularity  strengths,  and  control 
variables,  including  the  store  body  data  onto  TAPE12.  The  user 
is  responsible  for  disposition  of  this  dataset.  All  information 
is  written  unformatted. 

3 . 2  Input  Data 

This  section  of  the  report  will  describe  in  detail  the 
preparation  of  the  input  data  deck  for  Program  I,  LDCALC. 

Program  I  reads  the  data  required  to  compute  the  singularity 
strengths  modeling  all  components  of  the  parent  aircraft  con¬ 
figuration.  Data  which  are  repeated  in  both  Programs  I  and  II 
will  be  noted.  All  remaining  data  are  passed  from  Program  I  to 
Program  II  via  an  external  file  to  be  described  later. 
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Two  analytical  methods  may  optionally  be  used  to  compute 
flow  fields  around  bodies  in  these  programs.  For  purposes  of 
simplicity  of  nomenclature  the  technique  employing  line 
singularities  will  be  referred  to  as  the  circular  body  method. 
Similarly,  the  use  of  the  source  paneling  technique  for  the 
fuselage  will  be  referred  to  as  the  noncircular  fuselage  method. 

The  use  of  source  paneling  techniques  for  store  configurations 
will  be  referred  to  as  the  elliptical  store  method.  In  describing 
the  input  to  the  noncircular  and  elliptic  body  methods  the  same 
variables  and  input  routines  have  been  used  for  both  fuselage 
and  stores.  All  data  required  for  further  calculations  in 
Program  I  are  saved  on  external  files  via  internal  program  logic. 
The  data  required  by  Program  II  are  written  to  TAPE12  at  the  end 
of  Program  I.  ' 

3.2.1  Input  Formats 

The  format  for  the  input  data  for  Program  I  is  shown  in 
Figure  2.  Four  lines  of  information  are  shown  for  each  item. 

The  first  line  gives  the  item  number,  how  many  and  when  the  cards 
are  read,  and  the  routine  in  the  program  in  which  data  are  to  be 
read.  The  second  line  gives  the  program  variable  names,  the  third 
line  shows  the  card  column  fields  into  which  the  data  are  to  be 
punched,  and  the  fourth  line  shows  the  FORTRAN  format  type.  Data 
punched  in  I  and  E  formats  are  right  justified  in  the  fields 
whereas  data  in  F  format  can  be  punched  anywhere  in  the  field. 

A  decimal  point  should  be  included  in  both  E-  and  F-type  data. 

The  C  format  designation  in  columns  73  through  80  indicates  that 
the  fields  in  the  second  line  are  comments  not  read  by  the  program. 

3.2.2  Input  Descriptions  for  Program  I 

The  following  are  the  descriptions  of  the  variables  to  be 
read  in  Program  I  for  each  of  the  items  shown  in  Figure  2. 
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Item  number  1  is  an  index  NCARDS  which  indicates  how  many 
cards  of  information  are  to  follow  to  identify  the  run.  The 
value  of  NCARDS  must  be  one  or  greater. 

Item  number  2  is  a  set  of  NCARDS  cards  containing  hollerith 
information  identifying  the  run  and  may  start  and  end  anywhere 
on  the  card.  The  cards  are  reproduced  in  the  output  just  as 
they  are  read  in. 

Item  number  3  consists  of  one  card  and  contains  the  following 
information : 

ALFAC  fuselage  angle  of  attack,  degrees 

FMACH  Mach  number 

These  are  the  values  associated  with  the  parent  aircraft  used 
in  both  Programs  I  and  II.  These  parameters  are  again  read  in 
Program  II  and  compared  with  the  values  read  here.  If  they  do  not 
agree,  error  messages  are  written  and  Program  II  stops.  The  Mach 
number  should  be  between  1.2  and  3.0  and  the  angle  of  attack 
should  not  exceed  10°.  The  wing-fuselage  flow  model  is  valid 
within  these  limits. 

Item  number  4  contains  four  indices  which  specify  what  air¬ 
craft  components  are  present.  They  are 

NFU  fuselage?  NFU=0,  no 

NFU=1,  yes  -  use  circular  fuselage  method 
NFU=2,  yes  -  use  noncircular  fuselage  method 

NPY  pylon?  NPY=0,  no 

NPY»1,  yes 

NRACK  store  rack?  NRACK**0,  no 

NRACK-1,  yes 
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NSTRS  store?  NSTRS«0,  no 

NSTRS=1  to  7,  yes  -  total  number  of  stores 

Provision  has  been  made  for  omitting  the  fuselage.  For  cases 
with  no  fuselage  the  reference  coordinate  system  is  fixed  at  the 
wing  root  chord  leading  edge  and  ALFAC  should  be  the  wing  angle 
of  attack.  The  present  version  of  the  program  is  limited  to 
seven  circular  stores  of  different  shapes.  For  elliptic  store 
shapes,  however,  the  program  is  limited  to  two  different  elliptic 
shapes  for  which  there  may  be  multiple  stores  of  the  same  size 
and  shape.  A  total  of  seven  stores  are  permitted.  When  elliptic 
and  circular  stores  are  mixed,  two  elliptic  shapes  may  be  present 
with  as  many  circular  shapes  as  available  up  to  a  total  of  seven 
stores.  The  number  of  stores  can  be  zero.  In  addition,  the  program 
can  also  be  used  as  an  aid  in  determining  the  points  on  the  wing 
and  pylon  at  which  the  slopes  of  the  camberline  and  thickness 
distribution  are  to  be  input.  This  was  discussed  in  Section  3.1. 

Items  5  through  10  are  input  only  if  the  circular  fuselage 
method  is  selected,  that  is,  NFU«1.  For  all  other  options,  NFU=0 
or  2,  these  items  are  omitted  from  the  input  deck. 

Item  number  5  consists  of  three  fuselage  quantities  which 
follow  and  are  read  only  if  NFU=1. 

FLTHC  length  of  fuselage,  L,  feet 

FRMAX  maximum  fuselage  radius,  feet 

FTHSHK  maximum  nose  shock  wave  angle,  degrees 

The  latter  angle  is  used  as  the  limiting  value  at  the  tip  of  the 
nose  in  generating  the  modified  shock  wave  shape.  It  is  obtained 
from  Chart  5,  Reference  9,  as  a  function  of  the  semiapex  angle  of 
the  nose  tip.  Should  the  cone  angle  result  in  a  detached  shock, 
the  maximum  attached  shock  angle  should  be  used. 
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The  next  three  items  of  input,  items  6 ,  7,  and  8  of  Figure  2, 
describe  the  circular  fuselage  shape.  They  are 

Item  6: 

NFPOLY  number  of  polynomials  specifying  circular  fuselage 

shape,  1  <  NFPOLY  <  7 

Item  7; 

FXEND(J)  x/l  of  end  points  of  polynomials  specifying  fuselage 
shape,  NFPOLY  values 

Item  8 : 

FCOEF (J,K)  coefficients  of  polynomials  specifying  shape, 

K  *  1,2, ...7 

These  data  specify  the  radius  distribution  of  the  fuselage  and 
are  used  in  the  calculation  of  the  source-sink  distribution 
which  represents  the  fuselage  volume  and  the  doublet  distribu¬ 
tion  which  represents  the  fuselage  angle  of  attack  effects.  Up 
to  seven  polynomials  may  be  used.  The  polynomial  programmed  is 


where  through  C?  are  the  coefficients,  r  is  the  local  fuselage 
radius,  and  l  is  the  fuselage  length  of  Item  5.  The  polynomials 
must  be  input  for  shapes  which  are  made  dimensionless  by  the  body 
length  since  the  trajectory  program  is  written  assuming  this  to  be 
the  case. 

Item  6  specifies  the  number  of  polynomials.  Item  7  consists  of 
one  card  which  contains  the  NFPOLY  values  of  the  end  points  of  the 
polynomials  describing  the  shape.  The  decimal  point  can  be  placed 
anywhere  in  the  ten-column  field. 
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Item  number  8  is  a  set  of  NFPOLY  cards  specifying  the 
values  of  the  coefficients  of  the  polynomials,  Equation  (1).  All 
seven  coefficients  are  input  even  though  some  of  them  may  be  zero. 

Items  9  and  10  specify  the  body  interference  panel  layout  and 
the  number  of  line  sources  and  line  doublets  to  be  used  to  model 
the  fuselage  volume  and  angle  of  attack  effects.  The  various 
quantities  are 

Item  9: 

NCWB  number  of  rings  of  body  interference  panels 

NBDCRl  number  of  panels  in  a  ring  lying  above  the  wing 

on  the  left  half  of  the  fuselage 

NBDCR2  number  of  panels  in  a  ring  lying  below  the  wing  on 

the  left  half  of  the  fuselage 

NFSOR  number  of  line  sources  and  line  doublets  to  be 

used  over  the  fuselage  length 

Item  10; 

BODYPL  length  of  fuselage  over  which  body  interference 

panels  are  to  be  placed,  feet 

The  definitions  of  some  of  the  above  quantities  can  be  clarified 
by  the  use  of  Figure  3.  In  this  simplified  layout  there  are  five 
rings  of  body  interference  panels,  NCWB=5,  laid  out  over  the 
length  BODYPL.  In  each  of  the  rings  one  panel,  NBDCR1=1,  lies 
above  the  wing  zw  **  0  plane  and  three  panels,  NBDCR2=3,  lie  below 
this  plane.  Some  general  rules  for  determining  the  input  values 
of  these  quantities  will  now  be  given. 

First,  consider  the  length  BODYPL.  If  the  wing  trailing  edge 
is  supersonic,  that  is,  the  component  of  the  free-stream  Mach 
number  perpendicular  to  the  trailing  edge  is  greater  than  one, 
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body  interference  panels  should  be  laid  out  over  the  wing  root- 
chord  length,  CRW,  in  Figure  3.  For  a  subsonic  wing  trailing 
edge  {perpendicular  Mach  number  less  than  one) ,  body  interference 
panels  should  extend  to  the  xw  coordinate  of  the  wing  tip-chord 
trailing  edge.  If  the  base  of  the  store  being  separated  is  behind 
the  point  selected  using  the  preceeding  rules  then  the  body  inter¬ 
ference  panels  should  be  extended  to  the  xw  location  of  the  store 
base.  When  a  large  rearward  axial  motion  of  the  store  is  expected 
during  its  trajectory,  the  body  panels  should  be  extended  aft  to 
cover  this  motion. 

The  number  of  rings  of  body  interference  panels  and  the  number 
of  panels  in  each  ring  are  defined  as  follows.  The  number  of 
rings  is  dependent  on  the  number  of  panels  in  a  chordwise  row  on 
the  wing,  NCW.  This  quantity  is  input  as  part  of  item  36  and  its 
selection  will  be  discussed  later.  Since,  if  at  all  possible, 
the  same  number  of  chordwise  panels  should  be  used  on  both  the 
body  and  the  wing,  over  the  wing  root-chord  region,  CRW  in 
Figure  3,  the  number  of  rings  of  body  panels  can  be  found  using 
the  following  relation 

NCWB  =  NCW  (2) 

The  value  of  CRW  is  determined  as  described  under  the  item  35  in¬ 
put  data.  In  the  above  expression,  the  value  of  BODYPL  should  be 
adjusted  so  that  NCWB  is  an  integer. 

Two  rules  can  be  given  for  selecting  the  number  of  panels  in  a 
ring  on  the  left  side  of  the  body.  If  the  store  to  be  separated 
is  under  the  fuselage  at  least  eight  panels  should  be  used.  If 
the  store  is  under  the  wing,  six  panels  is  probably  sufficient. 

The  two  input  parameters  specifying  the  number  in  a  ring  are 
NBDCR1  and  NBDCR2.  These  are  the  number  of  panels  above  and 
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below  the  wing  zw  =  0  plane  of  Figure  3.  Assume  that  six  panels 
are  to  be  used  in  a  ring,  for  example.  Then,  the  values  of 
NBDCRl  and  NBDCR2  for  three  specific  wing  positions  are: 

Mid-wing 

NBDCRl  =  NBDCR2  =  3 

Wing  tangent  with  top  of  fuselage 
NBDCRl  =  0,  NBDCR2  =  6 

Wing  tangent  with  bottom  of  fuselage 
NBDCRl  =  6,  NBDCR2  =  0 

For  wings  located  intermediate  between  the  mid-wing  and  high  or 
low  wing  positions,  the  panels  should  be  divided  above  and  below 
the  wing  so  that  the  panel  widths  are  as  equal  as  possible. 

The  computer  program  as  documented  in  this  report  limits  the 
total  number  of  constant  u-velocity  panels,  including  interference 
panels,  on  the  fuselage,  wing  and  pylon  in  two  ways.  It  is  first 
limited  by  array  dimensions  to  a  maximum  of  200  panels  with  no 
more  than  20  in  a  ring  on  the  fuselage.  Thus,  the  following 
relations  must  be  satisfied. 

(NBDCRl  +  NBDCR2)  <20 

NCWB*  (NBDCRl  +  NBDCR2 )  +  NCW*MSW  +  NCP*MSP  <_  200  (3) 

Secondly,  the  technique  used  to  consolidate  panel  strengths 

at  corner  points  adds  additional  constraints  which  are  a  function 

of  the  number  of  breaks  in  sweep  or  dihedral  on  the  wing,  kw, 

and  the  number  of  breaks  in  sweep  on  the  pylon,  k  .  Thus,  the 

P 

following  relations  must  also  be  satisfied. 

MSW  +  k  +  MSP  +  k  +  2 (NBDCRl  +  NBDCR2 )  +  2  <  100  (4) 

w  p  — 
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and 


(MSW+k  +1) * (NCW+1)  +  (MSP+k  +1)*(NCP+1) 
w  p 

+  2 (NBDCR1+NBDCR2) * (NCWB+1)  £  500  (5) 

The  variables  NCW  and  MSW  are  input  in  item  36  and  the  variables 
NCP  and  MSP  are  input  in  item  44  and  have  individual  limitations 
defined  there.  The  number  of  wing  and  pylon  breaks  may  be 
determined  from  the  number  of  changes  in  sweep  or  dihedral  in 
items  37  and  45,  respectively. 

The  last  input  variable  in  item  number  9  is  NSFOR,  the  number 
of  line  sources  and  line  doublets  to  be  used  to  model  the  fuselage 
volume  and  angle  of  attack  effects.  The  general  rule  to  use  in 
determining  NFSOR  is  that  the  distance  between  the  origins  of 
successive  sources  or  doublets  should  be  approximately  equal  to 
the  chordwise  length  of  the  wing  constant  u-velocity  panels  at 
the  root  chord.  Therefore 


NFSOR  =  NCW 


FLTHC 

CRW 


<  100 


(6) 


The  maximum  value  of  100  is  imposed  by  dimension  statements  in 
the  program.  The  fuselage  length,  FLTHC,  was  input  as  part  of 
item  5.  The  wing  root  chord  length,  CRW,  will  be  input  as  part 
of  item  35  and  the  number  of  chordwise  panels,  NCW,  will  be  input 
as  part  of  item  36. 

Items  11  through  33  describe  the  fuselage  when  the  noncircular 
fuselage  method  is  selected  and  are  included  only  if  NFU=2. 

The  first  five  items  define  general  information  about  the 
body  size,  the  interference  shell,  and  set  execution  and  print 
options.  Two  sets  of  geometric  information  are  input  to  define 
the  fuselage  shape  and  lay  out  the  panels  on  the  surface.  The 
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first,  items  16  through  25,  define  the  external  shape  of  the 
fuselage  to  as  great  a  resolution  as  the  user  desires.  Items  26 
through  32  are  then  used  to  specify  the  subdivision  of  the 
external  shape  into  discrete  surface  panels.  The  panel  coordinates 
are  generated  from  interpolation  in  the  coordinates  specified 
or  computed  for  the  external  shape.  The  number  of  panels 
meridionally  or  axially  may  be  less  than,  equal  to,  or  greater 
than  specified  in  items  16  through  25. 

The  coordinates  used  in  the  definition  of  the  source  panels 
are  different  from  the  fuselage  body  coordinates  shown  in 
Figures  3  and  5.  For  the  purposes  of  laying  out  the  fuselage 
source  panels  all  input  will  be  defined  for  the  right  half  of  the 
configuration  as  shown  in  Figure  6  with  the  axial  direction 
positive  aft  from  the  nose,  the  y-direction  positive  out  the  right 
wing,  and  the  z-direction  positive  up. 

Items  11  and  12  are  used  to  determine  the  modified  shock  wave 
shape  and  to  lay  out  the  body  interference  panels. 

Item  11: 

FLTHC  length  of  fuselage,  feet 

FRMAX  maximum  lateral  dimension  of  fuselage,  feet;  used  as 

radial  location  for  first  shock  wave  shape  traverse 

FTHSHK  maximum  nose  shock  wave  angle,  degrees 

BODYPL  length  of  fuselage  over  which  body  interference 

panels  are  to  be  placed,  feet 

Item  12; 

NCWB  number  of  rings  of  body  interference  panels 

See  items  5,  9,  and  10  for  further  definitions  of  these  variables, 
except  FRMAX. 
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Item  13  consists  of  one  card  with  the  alphanumeric  descrip¬ 
tion  for  the  noncircular  fuselage  definition. 

In  items  14  through  33  columns  73-80  have  been  left  open  for 
user  comments  and  will  be  used  here  to  designate  card  names. 

Item  14  defines  some  of  the  optional  print  and  execution 
parameters  associated  with  the  source  paneling  of  the  noncircular 
fuselage. 

IXZSYM  symmetry  option;  for  the  fuselage  input  IXZSYM=0 

=0,  configuration  is  symmetric;  only  half  of  the 
body  is  paneled 

=1,  configuration  is  symmetric;  panel  both  halves 
of  body  from  data  of  first  half 
=-l,  configuration  is  asymmetric;  data  for  both 
halves  are  input 

IPRT  optional  print  control  parameter 

IPRT (1) =0 ,  do  not  print  copy  of  input  data 
IPRT (1) =1 ,  print  formatted  copy  of  input  data 

IPRT (2) =0,  do  not  print  supplementary  geometry 
data 

IPRT (2) =1 ,  print  panel  geometry,  angles,  and  areas 

IPRT(3)=0,  do  not  print  modified  shock  wave  shape 
IPRT (3) =1,  print  modified  shock  wave  shape  in 
BSHOCK 

IPRT (4 ) =0 ,  not  used 

IPRT (5 ) =0 ,  do  not  print  velocity  and  source  strength 
output 

IPRT (5) =1,  print  source  strengths  and  velocities 
used  in  boundary  condition  associated 
with  each  panel 
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IUVW 


component  velocity  influence  coefficient  calcula- 


NSHOCK 


MAXSHK 

NINLET 

NINBLK 


tion  option. 
IUVW=0,  yes 
IUVW=1,  no 


controls  number  and  where  modified  shock  wave 
shape  traverses  are  made 

NSHOCK=0,  a  single  traverse  is  calculated  at  90°, 
the  side  of  the  fuselage;  circular 
symmetry  is  assumed 

NSHOCK=positive,  NSHOCK  equally  spaced  shock 

traverses  are  computed  between  0°  and  90° 
no  symmetry  assumed 

NSHOCK=negative ,  NSHOCK  traverses  are  computed  at 
input  values  of  PHIS (K) ,  K=l, | NSHOCK |  in 
item  33;  no  symmetry  assumed 

| NSHOCK |  <  100/MAXSHK 

maximum  number  of  integration  steps  in  radial 
direction  used  to  generate  nonlinear  shock  shapes 

number  of  open  inlet  panels;  NINLET=0  if  no  inlet 

number  of  fully  blocked  inlet  panels;  NINBLK*0  if 
no  inlet 


NINVEL  number  of  additional  panels  used  in  inlet  velocity 

calculation;  NINVEL=0  if  no  inlet 
NINLET+NINBLK+NINVEL<25 


The  parameter  IXZSYM  is  to  allow  the  user  to  model  only 
half  of  the  body,  the  full  body  from  the  geometry  of  only  the 
right  hand  side,  or  a  full  nonsymmetric  configuration.  Since  the 
fuselage  assumes  symmetry  about  the  centerline,  only  the  right 
hand  side  is  modeled.  The  print  controls,  IPRT,  are  used  to 
control  output  during  execution.  The  output  from  the  input  phase 
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is  used  to  echo  the  data  read.  The  geometry  output  includes  the 
panel  coordinates,  control  points,  areas,  and  inclination  angles 
and  may  be  necessary  in  laying  out  inlet  input  to  confirm  panel 
index  numbers.  The  shock  output  includes  the  nose  and  inlet  shock 
shape  output.  The  velocity  print  control  is  used  only  to  examine 
the  influence  coefficient  matrix.  The  last  print  control  is  used 
to  display  the  panel  strengths  and  boundary  conditions.  The 
parameter  IUVW  is  used  to  suppress  the  storage  of  the  component 
velocity  coefficient  arrays.  These  are  only  required  during  the 
calculation  of  panel  pressures  and  loads.  These  are  not  needed 
for  the  fuselage  but  are  required  to  calculate  loadings  on  the 
separating  store.  The  number  of  nose  shock  shape  traverses, 

NSHOCK,  allows  the  user  to  specify  the  traverses  which  best 
describe  the  nonlinear  shape  of  the  shock  in  the  meridional  direc¬ 
tion  or  allow  the  program  to  specify  an  even  spacing.  MAXSHK 
limits  the  number  of  points  in  each  radial  traverse  to  minimize 
computing  time  and  also  acts  as  the  internal  first  dimension  in 
the  X  and  R  arrays  containing  the  shapes  of  the  nose  and  inlet 
shocks.  It  should  be  sufficient  to  allow  enough  points  to  be 
computed  to  describe  the  shape. 

For  the  fuselage,  the  suggested  values  of  the  above  parameters 
to  use  are  IXZSYM=0,  IPRT(1)=IPRT(2)=IPRT(3)=IPRT(5)*=1,  IPRT(4)=0, 
NSH0CK=3,  and  MAXSHK=8 .  The  parameter  IUVW  is  set  internally  to 
one  to  avoid  unnecessary  calculations.  The  card  is  designated 
the  IOPTS  in  columns  73-80. 

The  number  of  open  inlet  panels,  NINLET,  is  determined  from 
the  layout  over  the  face  of  the  inlet  for  panels  which  allow 
unimpeded  flow.  This  is  the  number  of  panels  assigned  special 
properties  in  items  28  and  29  during  the  solution  for  panel 
strengths.  NINBLK  is  the  number  of  fully  blocked  panels  laid 
over  the  face  of  the  inlet.  This  is  the  number  of  panels  assigned 
special  properties  in  items  28  and  29  which  are  used  to  model  the 
portion  of  the  flow  blocked  for  mass  flow  ratios  less  than  one. 
NINVEL  is  the  number  of  additional  panels  adjacent  to  the  inlet 
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panels  which  are  to  be  used  in  the  flow  field  calculations  used 
to  locate  the  inlet  shock  shape.  The  rules  for  selection  of 
these  panels  are  given  in  item  28.  The  combined  numbers  of 
NINLET+NINBLK+NINVEL  must  be  less  than  or  equal  to  25. 

Item  15  specifies  several  variables  which  locate  the  body 
interference  shell  and  modified  shock  wave  shape. 

XPIB  the  fuselage  x-station  at  which  the  fuselage  cross- 

sectional  shape  defines  the  body  interference  shell 
shape . 

the  fuselage  x-station  which  is  considered  to  be  the 
"shoulder"  of  the  body  nose.  If  set  equal  to  zero, 
the  value  is  computed  from  the  location  of  the  first 
occurrence  of  a  maximum  in  the  cross-sectional  area. 

angle  of  attack  correction  used  in  rotating  the  body 
nose  shock  wave  relative  to  the  fuselage  axis. 

The  body  interference  shell  lies  on  a  cylinder  (not 
necessarily  circular  in  cross  section)  of  constant  cross  section 
over  its  axial  length.  The  Y,Z  geometry  defining  the  surface  of 
the  shell  is  taken  to  be  equal  to  the  Y,Z  body  coordinates  at  a 
specified  x  station.  That  station  is  designated  by  XBIP.  It 
is  necessary  that  the  shell  be  either  on  or  outside  the  body 
source  panels  so  as  to  avoid  using  the  panel  solution  associated 
with  the  inside  of  the  source  panel.  Therefore,  the  body  x 
station  for  XBIP  should  be  that  at  which  the  largest  lateral 
dimensions  of  the  source  panels  occur. 

The  specification  of  the  shoulder  x-location,  XSHLDR,  is 
only  meant  to  provide  a  reasonable  estimate  of  the  point  where 
the  nose  stops  increasing  in  area.  It  is  used  only  to  estimate 
the  aft  point  where  the  modification  of  the  shock  wave  shape 


XSHLDR 


EALPHA 
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returns  to  the  original  linear  theory  value.  Typical  locations 
of  XBIP  and  XSHLDR  are  shown  pictorially  in  Figure  6. 

The  parameter  EALPHA  defines  the  amount  of  deviation  from 
a  straight  rigid  rotation  of  the  body  nose  shock  wave  about  the 
nose.  The  angle  from  the  axis  of  symmetry  of  the  nose  shock  to 
the  body  centerline  is  thus  defined  as  ag  =  ac (1-EALPHA) .  EALPHA 
is  limited  to  values  from  0.0  to  1.0  representing  the  range  of 
rigid  shock  rotation  to  no  shock  rotation  with  angle  of  attack. 
The  value  of  EALPHA  for  circular  cones  may  be  determined  as 
a  function  of  Mach  number  and  nose  semi -vertex  angle,  o,  from 
Table  14  in  Reference  10.  A  plot  of  the  information  in  Table  14 
is  also  found  in  Figure  4  of  this  report.  In  lieu  of  better 
information  for  noncircular  shapes  use  the  vertex  angle  of  an 
equivalent  circular  cone.  This  card  is  designated  XBIP  in 
columns  73-80. 

Item  16  contains  the  control  integers  used  to  specify  the 
type  and  the  amount  of  input  used  to  define  the  fuselage  external 
geometry.  The  variables  defined  here  and  in  the  next  nine  items 
control  the  specification  of  the  external  body  shape. 

JO  Read  reference  area  card;  for  fuselage,  J0=0 

J0=0,  do  not  read  item  17 
J0=1,  yes,  read  item  17 

J2  integer  flag  indicating  type  of  data  describing 

the  fuselage  external  geometric  shape 
J2=+l,  data  for  arbitrary  shaped  body  to  be  entered; 

include  section  data  in  items  20  and  21 
J2=-l,  data  for  circular  body  in  form  of  cross- 

sectional  areas  versus  XFUS  (item  18)  to  be 
entered;  include  section  data  in  item  22 
J2=-2,  data  for  circular  body  in  form  of  radius 

versus  XFUS  to  be  entered;  include  section 
data  in  item  23 
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J2=-3,  data  for  elliptic  body  to  be  entered  in 
form  of  semi-axes  in  y-  and  z-directions; 
include  items  23  and  24 

J2=-4 ,  data  for  elliptic  body  to  be  entered  in 

form  of  vertical  semi-axis.  A,  and  elliptic 
ratio,  B/A;  include  items  24  and  25 
J2=-5,  data  for  elliptic  body  to  be  entered  in 
form  of  horizontal  semi-axis,  B,  and 
elliptic  ratio,  A/B;  include  items  23  and  25 

J6  integer  indicating  whether  body  centerline  is 

cambered 

J6=0,  body  is  cambered;  if  J2<0,  include  item  19; 

if  J2=+l,  J6  must  be  input  as  0 
J6=+l,  body  is  uncambered;  omit  item  19 

NFUS  number  of  axial  body  segments;  1  _<  NFUS  jc  5 

NRADX(I)  number  of  points  in  Ith  body  segment  used  to  specify 

points  on  cross  section  about  the  periphery; 

1  <  I  _<  NFUS.  If  the  configuration  is  symmetric 
(IXZSYM=0, 1) ,  NRADX  is  input  for  the  half  section. 

If  the  entire  configuration  is  input  (IXZSYM=-1) , 
NRADX  is  input  for  the  full  section.  If  the  body 
is  circular  or  elliptic,  the  program  computes 
NRADX  Y-  and  Z-ordinates  about  the  section. 

NRADX (I)  <33. 

NFORX(I)  number  of  axial  stations  in  Ith  body  segment  at 

which  external  geometry  is  specified; 

NFUS 

1  £  I  <  NFUS,  l  NFORX(I)  <  51. 

1=1 

The  quantities  NRADX (I)  and  NFORX(I)  must  be  input  in  pairs. 
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The  reference  area  read  under  J0=1  may  also  be  read  later 
for  K0=1  in  item  30.  For  the  fuselage,  J0=0  and  K0=0  should  be 
used  since  no  force  or  moment  coefficients  are  computed.  The 
input  items  required  for  the  combinations  of  J2  and  J6  are 
shown  in  the  following  table. 


J2 

Required  Inputs  for  J2 

Required  Variable 
J6  Quantities 

and  J6 

Options 

Item  Numbers 

1 

0 

XFUS , 

Y,  Z 

H* 

00 

20, 

21 

-1 

±1 

XFUS , 

FUSARD 

18, 

22 

0 

XFUS, 

ZFUS ,  FUSARD 

18, 

19, 

22 

-2 

±1 

XFUS, 

FUSBY 

18, 

23 

0 

XFUS, 

ZFUS,  FUSBY 

18, 

19, 

23 

-3 

±1 

XFUS, 

FUSBY,  FUSAZ 

18, 

23, 

24 

0 

XFUS, 

ZFUS,  FUSBY, 

FUSAZ 

18, 

19, 

23, 

24 

-4 

±1 

XFUS, 

FUSAZ,  B/A 

18, 

24, 

25 

0 

XFUS, 

ZFUS,  FUSAZ, 

B/A 

18, 

19, 

24, 

25 

-5 

±1 

XFUS, 

FUSBY,  A/B 

18, 

23, 

25 

0 

XFUS, 

ZFUS,  FUSBY, 

A/B 

18, 

19, 

23, 

25 

The  quantities  NRADX  and  NFORX  specify  the  number  of  points 
in  the  meridional  and  axial  directions  per  segment  to  be  input  or 
computed  to  define  the  external  fuselage  shape. 

The  value  NRADX  for  the  Ith  segment  remains  constant  for  all 
sections  within  that  segment.  For  J2=l,  NRADX  represents  the 
number  of  y,z  pairs  per  section  which  are  input.  Long  constant 
geometry  sections  need  only  be  defined  once  at  each  end.  The 
values  of  NRADX  and  NFORX  may  vary  from  segment  to  segment.  NRADX 
and  NFORX  should  be  greater  than  or  equal  to  KRADX  and  KFORX 
defined  later  to  avoid  paneling  that  lies  inside  the  actual 
external  geometry  shape.  When  used  to  define  the  actual  panel 
layout  they  should  folio*’  the  rules  for  KRADX  and  KFORX.  KRADX 
and  KFORX  are  defined  in  item  27.  The  choice  of  NFORX (I)  must 
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further  satisfy  the  constraint  that  their  sum  for  all  segments 
must  also  be  less  than  or  equal  to  51.  In  selecting  the  number 
of  axial  stations  for  multiple  segment  layouts,  the  sections  at 
the  segment  interfaces  are  counted  twice.  They  are  counted 
as  the  trailing  stations  of  the  forward  segments  and  as  the  first 
station  of  the  aft  segments. 

This  card  is  designated  JCARD  in  columns  73-80. 

Item  17  contains  the  fuselage  reference  area.  It  is  input 
only  of  J0=1  in  item  16.  It  may  be  redefined  by  item  30.  If  not 
defined,  REFA  is  set  equal  to  1.0.  Since  fuselage  forces  and 
moments  are  not  calculated,  JO  in  item  16  should  be  input  as 
J 0=0.  If  this  is  done  this  input  item  is  skipped.  This  card 
is  designated  REFA  in  columns  73-80. 

The  values  in  items  18  through  25  define  the  external  shape 
of  the  fuselage  as  dictated  by  the  previous  control  variables. 

This  complete  set  of  items  is  repeated  for  each  of  NFUS  fuselage 
segments  as  required.  Segment  descriptions  must  be  given  in 
order  of  increasing  values  of  x. 

Items  18  and  19  specify  the  x  and  z  values  of  each  axial 
station  in  the  Kth  segment  as  follows.  NFORX(K)  values  are 
input  for  each  item. 

Item  18 : 

XFUS (I ,K)  axial  station  at  which  body  cross  section  data 

is  included,  feet;  K=1,NFUS;  I=l,NFORX(K) 

Item  19; 

ZFUS(I,K)  input  this  item  only  if  J2<0  and  J6=0  (see  item  16). 

z-value  of  camber  offset,  feet;  K=1,NFUS; 

1=1 , NFORX (K) 
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The  values  of  XFUS  define  the  axial  stations  at  which  subsequent 
section  data  are  input.  The  first  value  of  XFUS  (=0.0)  must  start 
at  the  body  nose.  All  values  must  increase  monotonically  to  the 
last  value  within  a  segment,  XFUS (NFORX (K) )  for  the  Kth  segment. 
These  values  are  shown  pictorially  in  Figure  6.  When  multiple 
segments  exist,  NFUS>1,  the  values  at  the  last  x  station  of 
NFORX (K-l)  segment  must  be  repeated  as  the  first  station  of  the 
NFORX (K)  segment.  For  circular  or  elliptic  sections,  the  body 
external  shape  at  a  given  x  station  is  described  first  by  the  y,z 
coordinates  of  the  shape  with  respect  to  a  local  "centerline" 
and  then  by  the  offset  of  this  local  "centerline"  with  respect 
to  the  body  axis.  This  is  the  camber  offset,  ZFUS.  The  cards 
are  designated  the  XFUS  and  ZFUS  cards  in  columns  73-80. 

If  J2,  item  16,  is  equal  to  one,  NFORX (K)  sets  of  items  20 
and  21  are  included  in  pairs.  They  define  I  descrete  points  on  the 
surface  of  an  arbitrary  body  as  follows: 

Item  20: 

Y (I , J , K)  y ^-ordinates  of  arbitrary  body  at  Jth  section  in 

Kth  segment,  feet;  1=1 ,NRADX (K) ,  J=l, NFORX (K) , 

K=1 ,NFUS 

Item  21: 

Z (I , J,K)  Zg-ordinates  of  arbitrary  body  at  Jth  section  in 

Kth  segment,  feet;  1=1 ,NRADX (K) ,  J=1 , NFORX (K) , 

K=1 ,NFUS 

NRADX(K)  values  of  first  Y,  then  Z,  are  input  to  define  the  half 
configuration  (IXZSYM=0  in  item  14)  as  required  for  each  station. 
When  ordering  the  points,  the  first  value  starts  on  the  vertical 
plane  of  symmetry  at  the  bottom  (Y=0,  Z-negative)  and  proceeds 
counterclockwise  about  the  body  when  looking  upstream.  If  the 
full  section  is  input,  the  first  and  last  values  must  be  identical 
in  order  to  close  the  body.  A  typical  section  is  shown  in  Figure  6 
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for  J2=l.  The  cards  are  designated  YJ  and  ZJ  where  J  is  the 
axial  station  number  in  columns  73-80. 


If  J2,  item  15,  is  less  than  zero,  user  requested  combinations 
of  items  22  through  25  are  used  to  define  the  circular  or  elliptic 
shape  for  the  Kth  segment  as  follows: 


Item  22: 

FUSARD  (I ,K)  cross-sectional  area  of  circular  body,  square  feet; 
1+1 , NFORX (K) ;  K=1,NFUS 


Item  23: 

FUSBY (I ,K)  for  J2=-2,  radius  of  circular  body,  feet;  for 

J2=-3  or  -5,  length  of  horizontal  semi-axis,  B, 
of  an  ellipse,  feet;  1=1 , NFORX (K) ;  K=1,NFUS 

Item  24: 

FUSAZ ( I , K)  length  of  vertical  semi-axis.  A,  of  an  ellipse, 

feet;  1=1, NFORX (K) ;  K=1,NFUS 


Item  25: 

ERATIO ( I , K)  for  J2=-4,  elliptic  ratio  B/A;  for  J2=-5,  elliptic 
ratio  A/B;  1=1, NFORX 


These  combinations  of  variables  define  the  possibilities  of 
specifying  the  standard  circular  and  elliptic  shapes  at  each 
axial  station.  The  permitted  combinations  of  these  items  are 
given  in  the  table  for  item  16.  When  defining  an  elliptic  shape 
(J2=-4  or  -5)  using  one  of  the  semi-axes  and  the  elliptic  ratio, 
only  the  first  value  of  the  elliptic  ratio  has  to  be  defined. 

All  subsequent  values  equal  to  zero  (not  defined)  are  set  equal 
to  the  previous  value.  The  elliptic  ratio  may  assume  a  new  value 
at  any  or  all  stations.  A  typical  section  is  shown  in  Figure  6 
for  J2=-3.  These  cards  are  designated  FUSARD,  FUSBY,  FUSAZ,  or 
A/B  in  columns  73-80,  respectively. 
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Item  26  contains  one  card  with  any  desired  identifying 
information  relating  the  paneling  sequence  to  the  previous 
external  geometry  shape. 


Item  27  defines  the  control  variables  used  to  specify  the 
resolution  of  the  configuration  into  source  panels.  This  card 
is  similar  to  item  16  with  the  difference  that  it  defines  the 
layout  of  the  source  panels  themselves  as  follows: 


KO  read  reference  area  and  length  card;  for  fuselage, 

K0=0 

K0=0,  no  exclude  item  30 

K0=1 ,  yes,  include  item  30 

KRADX(I)  number  of  meridian  lines  in  Ith  segment  used  to 

define  panel  edges  (I=1,NFUS).  There  are  three 
options  for  specifying  the  number  of  panel  edges. 

- 3  3 <KRADX<3  3 . 

KRADX (I) =0 ,  number  of  meridians  is  set  equal  to 
NRADX(I) 

KRADX ( I ) ^positive  number;  the  meridian  lines  are 
calculated  at  KRADX (I)  equally  spaced 
angles  about  the  body 

KRADX (I) “negative  number;  the  locations  of  meridian 
lines  are  input  at  KRADX (I)  values  of 
PHIK  in  item  31 

KFORX(I)  number  of  axial  stations  used  to  define  leading 

and  trailing  edges  of  panels  in  Ith  body  segment. 
Three  options  are  available 

KFORX(I)=0,  the  number  of  axial  stations  is  set 
equal  to  NFORX(I) 

KFORX (I) “positive  number;  panel  edges  are  defined 
at  KFORX(I)  values  of  XJ  defined  in 
item  32 
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KFORX (I) =negative  number;  |NFORX(I)|  equally  spaced 
panel  edges  are  defined  within  Ith  segment 

NFUS 

l  | KFORX (I) |  <  51 
1  =  1 

The  purpose  of  the  last  two  K-option  parameters  on  this  card 
is  to  specify  the  options  used  to  resolve  the  external  body  shape 
defined  in  items  16  through  25  into  source  panels.  The  user  is 
allowed  to  either  use  the  same  geometry  or  modify  the  original 
coordinates  in  laying  out  the  corners  of  the  source  panels.  This 
alows  the  user  to  easily  modify  the  panel  layout  without  changing 
the  original  shape  int-  it.  For  the  fuselage  KO  should  be  set  to 
zero,  since  no  loads  are  computed. 

For  symmetric  configurations  (IXZSYM=0 , 1) ,  KRADX  is  the 
number  of  meridians  on  the  half-body  including  meridians  at  0°  and 
180°.  For  full  configurations  (IXZSYM=-1) ,  KRADX  is  the  number  of 
meridians  on  the  full  body  including  the  meridians  at  0°  and  360°. 
However,  for  IXZSYM=1,  after  generating  the  panel  geometry, 

KRADX (1)  is  redefined  equal  to  the  number  of  merdiars  on  the  full 
section  for  internal  program  use.  The  sum  of  KFORX(I),  1=1, NFUS, 
must  be  less  than  or  equal  to  51.  The  number  of  panels  in  the 
circumferential  and  axial  directions  in  each  segment  is  thus  one 
less  than  KRADX (I)  and  KFORX (I),  respectively.  For  the  fuselage 
KRADX (I)  should  be  seven  or  greater  for  practical  results.  In 
selecting  the  axial  spacing  keep  the  ring  thickness  ratio 
(Ax/diameter)  less  than  1.0  for  the  fuselage. 

This  card  is  designated  the  KCARD  in  the  data  deck  in 
columns  73-80. 

Items  28  and  29  define  the  inlet  and  are  included  for  NINLET+ 
NINBLK+NINVEL  greater  than  zero.  Item  28  defines  the  number  of 
traverses  to  be  used  to  generate  the  inlet  shock  shape  and  defines 
the  panel  numbers  designated  as  inlet  panels  or  to  be  used  in  the 
inlet  velocity  calculations. 
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Item  28 : 


number  of  inlet  shock  shape  traverses  to  be 
computed  in  INLSHK;  1  <  NIS  £  80/MAXSHK  where  MAXSHK 
was  input  in  item  14 

panel  number  of  Ith  open  inlet  panel,  blocked  inlet 
panel,  or  additional  panel  used  in  inlet  velocity 
calculation.  First  NINLET  values  define  the  open 
inlet  panels;  the  next  NINBLK  panel  numbers  define 
the  blocked  inlet  panels  and  last,  the  next  NINVEL 
panel  numbers  specify  the  additional  panels  which  are 
used  in  inlet  velocity  calculations.  I=1,IMAX  where 
IV  X=NINLET4-NINBLK+NINVEL  <  25. 

The  JINLT (I :  panel  numbers  must  be  compatible  with  the  source 
panel  numbers  generated  by  the  program.  These  numbers  can  be 
det»  rmined  by  running  the  program  with  NINLET=NINBLK=NINVEL=0  and 
examining  the  output  obtained  by  setting  IPRT(2)=1  in  item  14. 

Item  29  defines  the  P  to  be  associated  with  the  inlet  panel 
special  properties. 

rojxj'  um  ;  associated  with  inlet  panel  solution 

The  user  i  required  to  lay  out  three  types  of  panels  in  the 
vicinity  of  the  inlet  for  use  in  the  program.  Each  of  these 
panels  have  special  properties  ascribed  to  them  during  the 
execution  of  the  program.  The  inlet  panels  as  defined  here  are 
those  placed  in  the  open  face  of  the  fuselage  inlet.  These  panels 
are  further  subdivided  into  open  and  blocked  inlet  panels  in 
accordance  with  the  intake  mass  flow  ratio  (MFR) .  The  first,  the 
open  inlet  panels,  are  placed  to  model  the  unimpeded  flow  into  the 
intake  and  enforce  the  condition  that  the  flow  must  be  equal  to 
the  free-stream  velocity.  The  second,  the  blocked  inlet  panels, 
are  placed  to  model  that  portion  of  the  mass  flow  (MFR< 1 )  which  is 
spilled,  or  blocked  from  entering  the  inlet.  Both  panels 


NIS 


JINLT (I) 
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propagate  their  linear  theory  influence  with  $  =  BTINLT.  This  is 
required  to  keep  the  panel  surface  aft  of  the  panel  Mach  cone  in 
order  to  keep  the  solution  numerically  tractable.  The  ratio  of 
open  panel  frontal  area  to  the  total  of  open  and  blocked  frontal 
area  defines  the  mass  flow  ratio  (MFR)  into  the  intake.  It  is 
the  responsibility  of  the  user  when  paneling  the  inlet  to  ensure 
that  the  physical  layout  of  open  and  blocked  panels  match  the 
desired  inlet  mass  flow  ratio  for  the  engine  and  flight  conditions 
at  hand. 

The  third  cype  of  panels  required  in  the  vicinity  of  the  inlet 
are  the  additional  panels  on  the  inlet  cowl  adjacent  to  and  aft 
of  the  intake.  These  panels  are  only  required  in  order  to  define 
the  subset  of  the  total  configuration  to  be  used  in  the  generation 
of  the  inlet  shock  shape.  Together  with  the  panels  covering  the 
face  of  the  intake,  these  additional  panels  form  the  set  of  panels 
used  to  compute  the  location  of  the  shock  emanating  from  the  inlet. 
During  this  calculation  they  are  assigned  the  property  that 
B  =  BTINLT.  For  the  remaining  computations,  their  influence 
is  propagated  at  B^.  The  minimum  set  of  panels  should  include  all 
panels  on  the  cowl  aft  of  the  intake  face  extending  past  the  inlet 
one  ring  of  panels.  The  quantities  defined  here  are  used  to  lay 
out  and  assign  the  special  properties  to  the  inlet  panels  in  the 
generation  of  the  nonlinear  shock  shape  associated  with  the 
oblique  shock  ramp  type  inlet  shown  on  the  next  page. 

The  first  value  in  item  28,  NIS,  determines  the  number  and 
therefore  the  location  of  the  traverses  defining  the  shock.  The 
traverses  are  positioned  in  a  sequence  predetermined  by  the  number 
NIS  as  follows.  The  first  traverse  is  positioned  at  y=YINLT  and 
(>  =  0°.  This  determines  the  shock  shape  in  the  y=YINLT  plane.  The 
remaining  traverses  are  computed  alternately  in  regions  I  and  II. 
The  second  traverse  propagates  from  y=YINLT  at  $=PHISHK.  The 
third  propagates  from  y=0  at  $  =  0°.  The  fourth  is  at  y=YINLT  and 
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REGION  III  ^ 


(XINLT.YINLT, , 
I  /ZINLT) 


_ ty.z)  | _ ^ ■ 


REGION  I 


REGION  II 


O  Open  inlet  panel 
□  Blocked  inlet  panel 

*  Additional  panels 
used  in  inlet 
calculations 


<f>=PHISHK/2  and  the  fifth  is  at  y=YCPI  and  <J>=0°.  Additional 
traverses  are  laid  out  by  subdividing  regions  I  and  II  alternately 
with  even  increments.  The  value  of  PHISHK  is  currently  internally 
set  to  90°  to  approximate  the  entire  region  under  the  wing  of 
an  aircraft  with  the  wing  tangent  to  the  top  of  the  inlet. 


The  inlet  panel  sequence  numbers  to  be  defined  here  are  those 
indicated  in  the  listing  of  the  panel  geometry.  The  panel  numbers 
in  the  first  NINLET  values  in  the  JINLT  array  are  the  numbers  of 
the  open  inlet  panels  placed  in  the  intake  face  through  which  air 
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may  flow.  The  next  NINBLK  values  in  the  JINLT  array  are  the 
numbers  of  the  blocked  inlet  panels  placed  on  the  intake  face 
which  have  been  positioned  according  to  the  engine  mass  flow 
ration  to  block  the  flow.  The  third  set  of  panel  numbers  in 
the  JINLT  array  are  the  additional  panels  covering  the  inlet 
cowling  immediately  aft  of  the  inlet  face  that  are  necessary 
for  use  in  the  inlet  velocity  calculation.  A  typical  set  is 
indicated  in  the  preceeding  sketch.  They  are  used  to  maintain 
the  continuity  of  the  linear  theory  velocity  calculations 
immediately  aft  of  the  inlet  during  the  determination  of  the 
shape  of  the  nonlinear  shock  emanating  from  the  inlet. 

The  maximum  S(BTINLT)  to  be  used  for  the  inlet  panels  is 

determined  as  a  function  of  and  the  2-D  inlet  ramp  angle,  6, 

or  deflection  angle,  from  Chart  2  of  Reference  9.  The  6  is 

computed  from  the  shock  turning  angle  0  as  BTINLT* tan (9O°-0  ) 

s  s 

as  shown  in  the  following  sketch.  This  slope  must  be  selected 
such  that  it  does  not  exceed  the  slope  of  the  panels  on  the  inlet 
face.  BTINLT  should  typically  be  2-3%  less  than  (ahead  of)  the 
corresponding  slope  of  the  inlet  face. 


Inlet  Side  View 
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For  additional  details  on  the  use  of  the  above  parameters 
and  the  methods  implemented  in  the  program  to  model  the  inlet 
see  Sections  3.2.3  and  4.2  of  Reference  1. 

Item  30  defines  the  reference  lengths  and  areas  used  in  the 
force  and  moment  coefficient  definitions.  This  item  is  omitted 
for  the  fuselage  since  K 0  in  item  27  should  be  input  as  zero. 

REFAR  reference  area,  square  feet.  If  greater  than  zero, 

redefines  REFA  in  item  17 

REFD  reference  length  or  diameter  used  to  nondimen- 

sionalize  moments,  feet.  If  blank,  program  sets 
value  equal  to  1.0 

REFL  body  length,  feet.  If  blank,  program  sets  value 

equal  to  XFUS (NFORX (NFUS) ) -XFUS (1 ) 

REFX  Xg- location  of  moment  center,  feet;  see  Figure  6 

REFZ  z^-location  of  moment  center,  feet;  see  Figure  6 

Item  31  defines  the  meridian  angles  for  the  edges  of  panels 
as  optionally  requested  by  KRADX(J)  for  each  body  segment.  This 
item  is  included  only  if  KRADX(J)  is  negative. 

PHIK(I,J)  meridian  angle  of  panel  edges  expressed  in  degrees; 
1=1, |KRADX(J) i , J=1,KFUS 

The  convention  is  observed  that  PHIK=0°  is  at  the  bottom  of  the 
body  and  PHIK=180°  is  at  the  top  of  the  body.  The  value  of 
PHIK (1 , J) =0. 0°  and  that  of  PHIK ( | KRADX (J) | ,J) =180°  or  360° 
depending  on  the  value  of  IXZSYM  input  in  item  14.  This  option 
should  be  used  to  ensure  that  panel  edges  in  the  presence  of 
a  wing  match  exactly  the  meridian  angle  at  the  wing-body 
junction.  Repeat  this  item  for  each  fuselage  segment.  This 
card  is  designated  PHIK  in  columns  73-80. 
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Item  32  defines  the  axial  stations  for  the  edges  of  panels. 
This  item  is  included  only  if  KFORX(J)  >  0. 

XJ(I,J)  x-stations  along  fuselage  defining  panel  leading 

and  trailing  edges  for  the  Jth  body  segment,  feet 

These  values  are  used  to  define  the  body  leading  and  trailing 
edges  of  adjacent  rings  of  panels.  When  selecting  values  in 
the  vicinity  of  "the  body  interference  panels,  values  of  XJ 
should  be  specified  which  correspond  exactly  to  the  leading  and 
trailing  edge  of  the  interference  shell.  This  card  is  repeated 
for  each  body  segment  and  is  designated  XJ  in  columns  73-80. 

Item  33  defines  the  location  of  the  traverses  which  form 
the  table  describing  the  modified  nose  shock  wave  shape.  These 
are  input  only  for  NSHOCK  <  0  in  item  14 . 

PHIS (I)  specified  angles  at  which  integration  or  modified 

shock  wave  shape  is  calculated,  degrees 

The  convention  is  used  such  that  PHIS(l)  is  zero  at  the 
centerline  on  the  bottom,  and  increases  counterclockwise. 

PHIS (NSHOCK)  only  needs  to  span  the  region  in  which  the  store 
may  travel.  Only  values  which  highlight  nonlinearities  in  the 
Rghock  versus  ♦  need  be  defined.  If  NSHOCK  >  0,  the  program 
generates  values  from  0°  to  the  angle  formed  by  a  radial  line 
through  the  wing  tip,  PHISHK.  The  angle  PHISHK  is  internally 
defined  as  <J>SHK  =  90°. 

Items  34  through  42  are  used  in  the  definition  of  the  wing 
and  body  interference  u-velocity  panels  and  the  wing  thickness 
panels. 

Item  34  contains  three  parameters  which  specify  the  wing 
location  relative  to  the  fuselage  nose  and  the  wing  incidence 
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angle.  The  three  parameters  are  shown  pictori&lly  in 
Figure  5. 

XBWOC  xfi  location  of  intersection  of  wing  leading  edge 

with  fuselage,  feet;  negative  as  shown  in  Figure  5 

ZBWO  location  of  intersection  of  wing  leading  edge 

with  fuselage,  feet;  negative  as  shown  in  Figure  5 

WIC  wing  incidence  angle  relative  to  fuselage  longi¬ 

tudinal  axis,  degrees;  positive  as  shown  in 
Figure  5 

Item  35  contains 

CRW  wing  root  chord 

SSPAN  wing  semispan,  feet 

The  definitions  are  shown  in  Figure  7 (a) .  The  wing  root 
chord  is  the  wing  chord  at  the  spanwise  station,  Y(l),  at  which 
the  wing  leading  edge  intersects  the  fuselage.  Both  quantities 
are  input  as  positive  quantities. 

Items  36  and  37  are  input  data  describing  the  left  wing 
panel  which  are  used  to  lay  out  the  constant  u-velocity  panels. 
The  quantities  are 


number  of  u-velocity  panels  in  a  chordwise  row 
on  wing;  NCW  >  4 

number  of  u-velocity  panels  in  a  spanwise  row 
on  wing;  MSW  £  19 

Item  37 

I  wing  u-velocity  panel  side-edge  number; 

I  ■  1  to  MSW+1 


Item  36 : 
NCW 

MSW 


Y(I) 


yw  location  of  Ith  side  edge  on  the  left  wing 
panel,  feet  (1=1  value  shown  in  Figure  7(a); 
negative  for  all  I's  since  on  the  left  wing  panel; 
measured  in  wing  planform  plane) 

PSIWLE(I)  -  leading-edge  sweep  of  wing  section  to  the  right 

of  the  Ith  side  edge,  degrees;  positive  swept  back 
(measured  in  wing  planform  plane).  PSIWLE (1) =0 . 0 

PSIWTE(I)  trailing-edge  sweep  of  wing  section  to  the  right 

of  the  Ith  side  edge,  degrees;  positive  swept  back 
(measured  in  wing  planform  plane).  PSIWTE (1) =0 . 0 

PHID(I)  dihedral  angle  of  wing  section  to  the  right  of  the 

Ith  trailing  leg,  degrees;  positive  up. 

PHID (1) =0 . 0 

Based  on  these  input  data,  the  wing  is  divided  chordwise  and 
spanwise  into  trapezoidal  shaped  constant  u-velocity  panels.  All 
the  NCW  panels  in  a  chordwise  row  have  equal  side  edge  chords 
and  spans,  the  spans  being  determined  by  the  Y(I)'s. 

The  question  arises  as  to  the  values  to  use  for  NCW  and 
MSW.  No  specific  rules  can  be  given  since  the  number  needed  is 
determined  to  some  extent  by  the  wing  planform  shape  and  the 
location  of  the  store  being  separated.  Fewer  panels  can  be  used 
on  the  wing  if  the  store  is  under  the  fuselage  than  if  the  store 
is  under  the  wing.  The  number  of  panels  in  a  chordwise  row, 

NCW,  is  also  determined  to  some  extent  by  the  camber  distribution, 
if  any,  since  the  camber,  item  39,  is  specified  at  the  panel 
control  points.  If  the  wing  is  uncambered  except  near  the 
leading  edge  then  a  fairly  large  number  of  panels  is  required 
in  a  chordwise  row  if  this  effect  is  to  be  included.  In 
general,  a  minimum  of  eight  panels,  NCW  =  8,  should  be  used. 

The  number  of  panels  in  a  spanwise  row,  MSW,  is  controlled 
to  some  extent  by  the  wing.  A  panel  side  edge  must  coincide 
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with  the  root  chord  and  each  break  in  leading-edge  sweep  angle, 
trailing-edge  sweep  angle  and  dihedral  angle,  and  if  a  pylon  is 
present,  a  trailing  leg  must  coincide  with  the  spanwise  location 
of  the  pylon.  One  must  also  coincide  with  the  wing  tip.  Con¬ 
sider  the  example  wing  alone  in  Figure  7 (b) .  There  is  a  break 
in  trailing-edge  sweep  at  yw/s  =  -0.2,  a  break  in  leading-edae 
sweep  at  -0.4,  a  pylon  at  -0.6,  and  breaks  in  dihedral  at  -0.6 
and  -0.8.  To  place  a  panel  side  edge  at  each  of  these  positions 
plus  the  wing  tip  and  the  root  chord  requires  five  panels  across 
the  semispan.  This  is  the  minimum  number  which  can  be  used  for 
this  wing  and  is  probably  not  sufficient.  Experience  with  the 
program  has  shown  that  in  some  cases  for  simple  wings  six 
panels  in  a  spanwise  row  have  given  good  results. 

The  only  sure  way  of  determining  convergence  with  number  of 
panels,  both  NCW  and  MSW,  is  to  examine  the  results  obtained 
from  the  trajectory  program.  For  a  particular  wing  various 
panel  layouts  should  be  tried  to  assure  convergence.  The 
minimum  number  of  panels,  consistent  with  the  desired  accuracy, 
should  be  used  in  order  to  minimize  the  trajectory  calculation 
time. 


The  maximum  number  of  panels  that  can  be  placed  on  the  left 
wing  panel,  pylon,  and  fuselage,  is  200.  This  limit  is  imposed 
by  dimension  statements  in  the  computer  program.  Thus,  limits 
may  be  imposed  on  NCW  and  MSW  in  order  to  satisfy  Equations  (3) , 
(4) ,  and  (5) .  These  equations  are  presented  in  the  discussion  of 
items  9  and  10. 

Item  number  37  consists  of  a  deck  of  MSW+1  cards.  The 
index  I  is  the  panel  side-edge  number.  The  side  edges  are 
numbered  from  the  root  chord,  1=1,  to  the  tip,  I=MSW+1. 

Associated  with  each  I  are  the  spanwise  location  of  the  side  edge, 
Y ( I ) ,  the  sweep  angles  of  the  leading  and  trailing  edges  of  the 
wing  segment  to  the  right  of  the  side  edge,  PSIWLE(I)  and 
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PSIWTE(I),  and  the  dihedral  angle  of  the  section  to  the  right, 
PHID(I).  When  1=1  these  angles  are  input  as  zero. 

The  two  indices  of  item  number  38  are  associated  with  the 
wing  twist  and  camber  distribution. 

Item  38 ; 

NTAC  twist  and/or  camber? 

NTAC=0 ,  no 
NTAC=1,  yes 

NUNI  if  wing  has  no  twist,  and  the  camber  distribution 

is  similar  at  all  spanwise  stations,  NUNI=1;  for 
all  other  cases,  NUNI=0  (omit  if  NTAC=0) 

If  NTAC=1,  item  number  39  is  included  in  the  input  data  deck. 
These  data  specify  the  wing  twist  and/or  camber  distribution  in 
terms  of  the  tangent  of  the  local  angle  of  attack  of  the  camber¬ 
line  for  a  wing  root  chord  angle  of  attack  of  zero  degrees.  The 
function  of  the  index  NUNI  is  explained  below. 

Item  39: 

ALP HAL ( J)  tan  of  the  wing  camberline  at  the  constant 

u-velocity  panel  control  points.  If  NUNI-1,  only 
data  for  the  chordwise  row  adjacent  to  the  root 
chord  are  input.  The  first  value  is  for  the  control 
point  nearest  the  leading  edge.  If  NUNI=0,  data  for 
all  chordwise  rows  must  be  input  starting  nearest 
the  root  chord  and  working  outboard.  Data  for  each 
.  row  start  on  a  new  card  (omit  if  NTAC=0) 

The  constant  u-velocity  panel  control  points  are  at  95  percent 
of  the  chord  which  passes  through  the  centroid  of  area  (see 
Figure  5,  Reference  1)  of  each  elemental  panel  laid  out  by  NCW, 
MSW,  and  Y ( I ) ' s  of  items  36  and  37. 
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The  values  of  ALP HAL { J)  are  obtained  as  follows.  Consider 
the  following  sketch  which  shows  the  cambered  and  twisted  section 


of  the  lifting  surface  at  some  spanwise  station  for  zero  wing 
angle  of  attack.  At  the  control  point  P,  a  tangent  to  the 
camberline  is  constructed,  which  makes  an  angle  a ^  with  the  root 
chord  (the  xw  axis) .  The  positive  sense  of  is  shown.  The 
input  value  required  is  ALP HAL ( J)  =  tan  a^.  For  wings  which  have 
the  same  camber  distribution  at  all  spanwise  stations  and  no 
twist,  NUNI=1,  data  are  only  input  for  the  row  of  control  points 
closest  to  the  root  chord.  The  program  assigns  these  values 
to  all  other  rows. 

The  three  indices  of  item  number  40  are  associated  with 
the  specification  of  the  wing  thickness  distribution.  They  are 

Item  40: 

NCWS  number  of  thickness  panels  in  a  chordwise  row  on 

the  wing 

MSWS  number  of  thickness  panels  in  a  spanwise  row  on 

the  wing;  MSWS  1  19 

NUNIS  if  wing  has  a  similar  thickness  distribution  at 

all  spanwise  stations,  NUNIS=1;  if  not,  NUNIS=*0 
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The  thickness  panels  are  laid  out  on  the  wing  independently 
from  the  u-velocity  panels.  Three  constraints  limit  the  choice 
of  the  number  of  panels  used.  First,  a  total  of  400  thickness 
panels  can  be  used  on  the  wing  and  pylon; 

MSWS*NCWS  +  MSPS*NCPS  <  400  (7) 

Similar  to  the  wing,  the  technique  of  consolidation  of  strengths 
at  panel  corners  imposes  the  two  additional  limits  based  on 
dimensions  as  follows 


MSWS  +  k,  +  MSPS  +  k  +  2  <  60  (8) 

w  p  — 

and 

(MSWS+k  +1) (NCWS+1)  +  (MSPS+k  +1) (NCPS+1)  <  1000  (9) 

W  p 

where  k  is  the  number  of  breaks  in  sweep  and  dihedral  on  the 
w 

wing  and  kp  is  the  number  of  breaks  in  sweep  on  the  pylon.  The 
number  of  spanwise  panels  required  is  a  function  of  the  spanwise 
variation  of  the  thickness  distribution  and  the  number  of  breaks 
in  sweep  or  dihedral.  If  the  distribution  of  t/c  with  x/f.  is 
constant  in  the  spanwise  direction  (NUNIS=1)  and  there  are  no 
breaks  in  sweep  or  dihedral,  only  one  row  of  panels  is  required. 
If  the  thickness  is  not  similar  and  there  are  breaks  in  sweep 
and/or  dihedral,  more  spanwise  panel  rows  should  be  used  as 
necessary.  Each  of  these  chordwise  strips  is  divided  into  NCWS 
equal  chord  panels.  Usually  more  thickness  panels  should  be  used 
in  a  chordwise  row  than  u-velocity  panels.  Experience  with  the 
program  has  shown  12  to  14  thickness  panels  usually  to  be 
sufficient.  Again,  this  can  only  be  checked  by  varying  the 
number  and  examining  the  resulting  store  load  distributions. 

The  minimum,  consistent  with  the  desired  accuracy,  should  be  used 
to  minimize  the  trajectory  calculation  time. 
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Item  41  defines  the  wing  planform  information  for  the 
thickness  panel  layout  in  a  manner  similar  to  that  in  item  37. 
The  variables  are 


Item  41: 
I 

YS(I) 


PSWSLE (I) 


PSWSTE(I) 


PHIS (I) 


wing  thickness  panel  side-edge  number;  1=1  to  MSWS+1 

yw  location  of  Ith  side  edge  on  the  left  wing  panel, 
feet  (1=1  values  same  as  for  item  37  shown  in 
Figure  7(a),  negative  for  all  I's  since  on  left 
wing  panel;  measured  in  wing  planform  plane) 

leading-edge  sweep  of  wing  section  of  thickness 
panels  to  the  right  of  the  Ith  side  edge,  degrees. 
Positive  swept  back  (measured  in  wing  planform 
plane) .  PSWSLE (1) =0.0 

trailing-edge  sweep  of  wing  section  of  thickness 
panels  to  the  right  of  the  Ith  side  edge,  degrees. 
Positive  swept  back  (measured  in  wing  planform 
plane).  PSWSTE  (1) =0 . 0 

dihedral  angle  of  wing  section  of  thickness  panels 
to  the  right  of  the  Ith  trailing  leg,  degrees; 
positive  up.  PHIS (1) =0.0 


This  information  is  identical  in  form  to  that  used  to  define 
the  u-velocity  panels  on  the  wing.  Though  the  identical  informa¬ 
tion  could  be  used  for  both  layouts,  it  is  recommended  that  more 
panels  be  used  in  the  chordwise  direction  and  only  as  many  in 
the  spanwise  direction  as  required  to  account  for  breaks  in 
sweep,  dihedral,  or  changes  in  thickness  distribution.  A 
single  strip  of  panels  is  acceptable  when  no  such  changes  occur. 


Item  number  42  contains  the  thickness  distribut  on. 
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Item  42: 


THETAL(J)  slope  of  the  wing  thickness  distribution  at  the 
centroid  of  area  of  the  thickness  panels.  If 
NUNIS=1  only  data  for  the  chordwise  row  adjacent 
to  the  root  chord  are  input.  The  first  value  is 
for  the  panel  at  the  leading  edge.  If  NUNIS=0, 
data  for  all  chordwise  rows  must  be  input  starting 
at  the  root  chord  and  working  outboard.  Data  for 
each  row  start  on  a  new  card 

Note  that  the  values  of  the  thickness  slopes  are  input  for  the 
centroid  of  area  of  each  panel.  Also,  for  wings  with  similar 
distributions  at  all  spanwise  stations,  NUNIS=1,  data  are  only 
input  for  the  row  of  panels  adjacent  to  the  root  chord.  The 
program  assigns  these  values  to  all  other  rows. 

The  values  of  THETAL(J)  are  obtained  as  follows.  Consider 
the  following  sketch  which  shows  the  thickness  envelope  at  a 


Thickness 

envelope 


Mean  camber 
surface 


Tangent  to  thickness 
envelope  at  point  P 


spanwise  station  where  the  slopes  are  to  be  determined.  At  a 
point  P  at  which  the  slope  is  required,  a  tangent  to  the  surface 
of  the  thickness  envelope  is  constructed  which  makes  an  angle  9 
with  the  line  connecting  the  leading  and  trailing  edges  of  the 
envelope.  The  input  value  required  is  THETAL(J)  =  tan  9. 

Forward  of  the  point  of  maximum  thickness  9  is  positive  and  aft 
of  this  point  it  is  negative.  In  some  cases  for  a  wing  with  a 
blunt  leadina  edae,  the  thickness  slope  should  not  exceed  a 
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certain  value.  The  determination  of  this  value  is  discussed  in 
detail  in  Section  3.4.1  of  Reference  2. 


The  next  six  items  of  input  data,  items  43  through  48,  are 
associated  with  the  pylon  if  one  is  present,  NPY=1.  If  NPY=0 
items  43  through  48  are  omitted  from  the  input  deck.  Both  the 
u-velocity  and  thickness  panels  on  the  pylon  follow  the  format 
and  limitations  prescribed  for  the  wing  panels.  Item  43  consists 
of  one  card  and  contains  the  following  quantities: 

Item  43: 

IP  index  of  the  y^  location  of  the  pylon.  This  must 

be  one  of  the  Y{I)'s  read  in  for  the  wing,  item  37, 
if  under  the  wing.  If  the  pylon  is  under  the 
fuselage  centerline,  IP=0 

CRP  length  of  pylon  root  chord,  feet 

HP  height  of  pylon  measured  from  wing  chordal  plane, 

feet 

XPLE  location  of  pylon  root  chord  leading  edge  measured 

from  local  wing  chord  leading  edge,  feet  (negative 
behind) 

The  pylon  location  can  be  under  the  fuselage,  IP=0,  or  at  one  of 
the  wing-panel  side-edge  locations  outboard  of  the  wing-fuselage 
juncture.  The  remaining  variables  are  shown  in  Figure  8  for  the 
two  cases  with  only  one  spanwise  panel  on  the  pylon.  Note  from 
Figure  8  that  for  a  pylon  under  the  fuselage,  XPLE  is  measured 
from  the  wing  root  chord  leading  edge  and  HP  is  measured  from  the 
bottom  of  the  fuselage. 

The  next  two  items  of  input  data  are  associated  with  the 
constant  u-velocity  panels  which  represent  the  pylon  loading. 

Item  number  44  is  one  card  and  contains: 


Item  44: 


NCP 

MSP 


number  of  panels  in  a  chordwise  row  on  the  pylon 

number  of  panels  in  a  spanwise  row  on  the  pylon 
MSP  <  19 


For  a  typical  pylon,  two  or  three  panels  in  a  spanwise  row  are 
sufficient,  MSP=2  or  3.  The  number  chordwise,  NCP,  depends  on 
the  length  of  the  pylon  root  chord.  The  chordwise  dimensions 
of  the  trapezoidal  shaped  area  elements  on  the  pylon  should  be 
approximately  the  same  as  those  on  the  wing  immediately  above 
the  pylon.  That  is,  the  local  wing  chord  length  divided  by  NCW 
of  item  36  should  be  approximately  equal  to  the  pylon  root  chord 
length  divided  by  NCP. 

The  program  is  limited  to  200  constant  u-velocity  panels  on 
the  wing-pylon-fuselage  combination.  NCP  and  MSP  must  thus 
follow  the  restrictions  imposed  by  program  dimensions  as  specified 
in  Equations  (3) ,  (4) ,  and  (5) . 

Item  number  45  consists  of  a  deck  of  MSP+1  cards  which 
contain  the  following  information. 


Item  45: 

K 

Z  (K) 

PSIPLE(K) 


PSIPTE(K) 


pylon-panel  side-edge  number;  K  =  1  to  MSP+1 

z  location  of  the  Kth  side  edge,  feet;  measured 
from  local  wing  chord 

sweep  angle  of  the  pylon  leading  edge  above  the  Kth 
side  edge  in  the  wing  coordinate  system,  degrees. 
Positive  sweep  is  swept  back.  PSIPLE (1) =0 . 0 

sweep  angle  of  the  pylon  trailing  edge  above  the 
Kth  side  edge  in  the  wing  coordinate  system,  degrees. 
Positive  sweep  is  swept  back.  PSIPTE (1) =0. 0 
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The  first  side  edge  should  be  placed  on  the  pylon  root  chord 
and  the  last  on  the  tip  chord.  The  remaining  side  edges  are 
spaced  between  these  two.  For  a  pylon  located  under  the 
fuselage,  Figure  8 (a),  Z(l)  should  be  the  zw  location  of  the 
pylon  root  chord,  the  bottom  of  the  fuselage,  and  Z(MSP+1)  should 
equal  Z(l)  +  HP.  For  a  pylon  under  the  wing,  Z(l)  =  0  and 
Z (MSP+1 )  =  HP. 

Items  46,  47,  and  48  provide  data  required  to  model  the 
pylon  thickness  distribution.  These  items  follow  the  same  format, 
basic  descriptions,  and  limitations  imposed  on  the  wing  thickness 
panels  in  items  40  through  42  with  the  exception  that  no  provision 
is  made  for  dihedral  effects. 

Item  46 : 

NCPS  number  of  thickness  panels  in  a  chordwise  row  on 

the  pylon 

MSPS  number  of  thickness  panels  in  a  spanwise  row  on  the 

pylon;  MSPS  £  10 

NUNIP  if  pylon  has  a  similar  thickness  distribution  at 

all  spanwise  stations  and  there  are  no  breaks  in 
sweep,  NUNIP=1,  if  not,  NUNIP=0 

Item  47  consists  of  a  deck  of  MSPS+1  cards  which  contain 
the  following  information. 

Item  47 ; 

K  pylon  thickness  panel  side-edge  number;  K=1  to 

MSPS+1 

ZS  (K)  z  location  of  the  Kth  thickness  panel  side  edge, 

feet;  measured  from  local  wing  chord 
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PSPSLE (K) 


sweep  angle  of  the  pylon  leading  edge  above  the  Kth 
thickness  panel  side  edge  in  the  wing  coordinate 
system,  degrees.  Positive  sweep  is  swept  back. 

PSPSLE (1 ) =0 . 0 . 

PSPSTE(K)  sweep  angle  of  the  pylon  trailing  edge  above  the 

Kth  thickness  panel  side  edge  in  the  wing  coordinate 
system,  degrees.  Positive  sweep  is  swept  back. 

PSPSTE (1) =0 . 0 

The  placement  of  the  panel  side  edges  follows  the  same  rules  as 
for  locating  the  u-velocity  panels.  Panel  edges  are  required  at 
all  breaks  in  sweep  and  changes  in  thickness  distribution. 

Item  48  defines  the  thickness  envelope  of  the  pylon  in  the 
same  manner  shown  for  the  wing  in  item  42  as  follows. 

I tem  4  8 : 

THETPL(J)  slope  of  the  pylon  thickness  distribution  at  the 
centroids  of  the  thickness  panels.  If  NUNIP=1 
only  data  for  the  chordwise  row  adjacent  to  the  root 
chord  are  input.  The  first  value  is  for  the  panel 
at  the  leading  edge.  If  NUNIP=0,  data  for  all 
chordwise  rows  must  be  input  starting  at  the  root 
chord  and  working  outboard .  Data  for  each  row 
start  on  a  new  card. 

These  data  are  prepared  in  the  same  manner  as  were  the  correspond¬ 
ing  data  for  the  wing  thickness,  item  42.  Comments  made  there 
concerning  a  blunt  leading  edge  also  apply  to  the  pylon. 

Choice  of  the  number  of  thickness  panels  on  the  pylon  is 
still  limited  by  the  number  of  panels  on  the  wing-pylon  combina¬ 
tion  and  should  always  be  kept  at  the  minimum  consistent  with 
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the  accuracy  of  the  results.  NCPS  and  MSPS  must  thus  follow  the 
restrictions  imposed  by  the  program  dimensions  as  specified  in 
Equations  (7) ,  (8) ,  and  (9) . 

The  next  five  items  of  input  data  locate  and  describe  the 
ejection  rack  if  one  is  present  (NRACK  of  item  4  is  equal  to  one) . 
If  NRACK=0,  these  five  items  are  omitted.  If  a  rack  is  present 
the  preceding  data  for  the  pylon  must  have  been  input  since  the 
program  has  been  written  assuming  that  if  there  is  a  rack  there 
is  also  a  pylon.  The  rack  data  to  be  input  only  model  the  body 
of  the  rack. 

Item  49  is  one  card  and  locates  and  sizes  the  rack  from  the 
following  information. 

Item  49; 

RLTHC  length  of  rack,  feet 

RRMAX  maximum  rack  radius,  feet 

XRNC  x^  location  of  rack  nose  measured  from  local  wing 

chord  leading  edge,  feet;  positive  ahead.  The  x^, 
z^  coordinate  system  is  shown  in  Figure  8 

ZRN  z^  location  of  rack  nose  measured  from  local  wing 

chord  leading  edge,  feet;  positive  below.  The  x^, 
z^  coordinate  system  is  shown  in  Figure  8 

RIC  rack  incidence  angle  measured  relative  to  wing  root 

chord,  degrees;  positive  nose  up 

Item  numbers  50,  51,  and  52  define  the  shape  of  an  equivalent 
circular  body  used  to  model  the  rack. 

Item  50; 

NRPOLY  number  of  polynomials  specifying  the  circular  rack 

shape;  1  <  NRPOLY  <  7 
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Item  51: 


RXEND(J)  X/l  of  end  points  of  polynomials  specifying  rack 

shape,  NRPOLY  values 

Item  52: 

RCOEF ( J , K)  coefficients  of  polynomials  specifying  shape 

These  data  specify  the  radius  distribution  of  the  equivalent 
circular  body  rack  and  are  used  in  the  calculation  of  the  source- 
sink  distribution  which  represents  the  rack  volume  and  the 
doublet  distribution  which  represents  the  rack  angle  of  attack 
effects.  The  polynomials,  given  by  Equation  (1) ,  follow  the  same 
rules  and  definitions  as  used  to  define  the  circular  body  fuselage 
in  items  6  through  8 . 

Item  53  is  one  card  and  contains  the  following: 

Item  53: 

NRSOR  number  of  rack  sources;  NRSOR  £  100 

RTHSHK  maximum  shock  wave  angle  at  nose,  degrees 

The  number  of  sources,  NRSOR,  follows  the  same  general  rules  used 
to  specify  the  source  distribution  of  the  fuselage  in  item  9.  The 
angle,  RTHSHK,  is  used  as  the  limiting  value  at  nose  in  generating 
the  modified  shock  wave  shape.  It  is  obtained  from  Reference  9, 
and  follows  the  same  rules  used  to  define  FTHSHK  for  the  fuselage 
in  item  11. 

If  ro  store  is  present  (NSTRS  of  item  4  is  equal  to  zero) 
this  concludes  the  input  data  deck.  If  NSTRS  is  not  zero,  the 
remaining  items  in  the  data  deck  for  Program  I  are  used  to  describe 
each  of  up  to  seven  circular  or  elliptic  shaped  stores.  The 
information  included  here  describes  the  stores  in  their  carriage 
position.  Information  pertaining  to  the  separated  store  will  be 
included  in  Program  II. 
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The  number  of  stores  present  is  defined  in  item  4.  The  stores 
are  limited  to  two  types,  circular  and  elliptic.  The  circular 
store  method  employs  the  line  source  and  doublet  analysis  as  used 
for  the  circular  fuselage.  The  method  is  the  most  efficient 
numerically,  but  is  limited  to  the  consideration  of  only  one  set 
of  monoplane  or  cruciform  fins  in  Program  II.  The  elliptic 
store  method  uses  the  source  paneling  method  of  analysis.  A 
circular  store  can  also  be  modeled  by  this  method.  Though  many 
routines  used  are  common  with  the  fuselage  description  for  a 
noncircular  body,  a  few  subroutines  still  exist  restricting  this 
portion  of  the  analysis  to  consideration  of  only  elliptic  store 
body  shapes.  This  store  option  also  must  be  used  when  considering 
multiple  sets  of  wing  or  fin  combinations. 

Item  54  consists  of  NSTRS  cards,  one  for  each  store,  con¬ 
taining  the  following  information  for  J=1  to  NSTRS. 


Item  54 : 

NUMSTR(J) 

NSHAPE(J) 


SLTHC(J) 

SRMAX(J) 

XSNC(J) 


YSN(J) 

ZSN(J) 


store  number;  <99 

shape  number  of  store.  If  NSHAPE (J)£50,  circular 
store  body  option  is  used;  if  5l£NSHAPE ( J)£99 , 
elliptic  store  body  option  is  used 

length  of  store,  feet 

maximum  radius  of  store  feet;  SRMAX=0.0  for  elliptic 
store  option 

x^  location  of  store  nose  measured  from  wing  chord 
leading  edge  immediately  above  store,  feet;  positive 
ahead 

yw  location  of  store  nose  measured  from  fuselage 
centerline,  feet;  positive  to  the  right 

location  of  store  nose  measured  from  wing  chord 
leading  edge  immediately  above  store,  feet;  positive 
below 
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SIC(J) 


store  incidence  angle  measured  relative  to  wing 
root  chord,  degrees;  positive  nose  up 

SPHIR(J)  roll  angle  of  store  body  coordinate  system  relative 

to  inertial  or  fuselage  coordinates,  degrees;  posi¬ 
tive  right  wing  down  viewed  from  the  rear 

The  shape  number  is  used  to  distinguish  between  circular  and 
elliptic  analysis  options.  Shape  numbers  greater  than  50  are 
arbitrarily  reserved  to  designate  stores  to  be  modeled  using  the 
elliptic  store  option.  A  total  of  seven  stores  and  seven  shapes 
may  be  used.  Only  two  of  those  shapes  may  be  elliptic,  though 
multiple  stores  of  those  shapes  may  exist.  The  above  coordinates 
and  orientations  define  the  locations  of  the  stores  in  their 
carriage  positions.  For  the  ejected  store  in  Program  II,  the 
initial  values  for  angles  0  and  4>  are  set  to  SIC  and  SPHIR.  The 
roll  orientation  of  the  store  coordinates  should  typically  be 
zero  for  circular  body  option.  It  may  be  nonzero  for  the  elliptic 
option  where  alignment  of  the  store  coordinates  with  the  particular 
carriage  position  is  required. 

To  avoid  redundant  input  data  associated  with  multiple  stores 
with  the  same  shape,  the  remaining  store  data  is  included  only 
once  for  each  shape  appearing  in  item  54.  Item  number  55  is  one 
card  which  contains: 

Item  55: 

NSHPT  number  of  different  values  of  NSHAPE(J)  from  item 

number  54;  1  £  NSHPT  <  7.  Only  two  shapes  may  have 
NSHAPE(J)  greater  than  50 

One  data  deck  of  items  56  through  76  follows  for  each  of  the 
NSHPT  store  shapes  as  dictated  by  NSHAPE(J),  J=1,NSTRS.  For 
stores  modeled  using  the  circular  store  option,  NSHAPE(J) <50,  only 
items  56  through  59  are  required.  For  stores  modeled  using  the 
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elliptic  store  option,  51<NSHAPE (J)<99,  items  56  and  60  through 
76  are  input  as  required. 

Item  number  56  is  the  first  shape  card  and  contains  the 
following  s 

Item  56: 

MSHAPE(J)  shape  number  of  Jth  store  shape; 

if  MSHAPE { J) <50 ,  circular  store  option  and  data  to 
be  used; 

if  £51MSHAPE {J)£99 ,  elliptic  option  and  data  to 
be  used 

MSOR  number  of  line  sources  and  line  doublets  to  be  used 

to  model  the  store  volume  and  angle  of  attack 
effects;  MSOR£lOO.  Value  required  only  if 
MSHAPE (J) <50 

STHSHK  maximum  nose  shock  wave  angle,  degrees 

The  variable  MSOR  is  required  only  for  the  circular  store  option. 
A  value  of  MSOR  of  30  to  40  should  be  sufficient  to  model  a 
circular  store.  The  value  can  be  varied  to  determine  whether 
increasing  the  number  affects  the  trajectory.  The  angle,  STHSHK, 
is  used  as  the  limiting  value  at  nose  in  generating  the  modified 
shock  wave  shape.  It  is  obtained  from  Reference  9,  and  follows 
the  same  rules  used  to  define  FTHSHK  for  the  fuselage  in  item  11. 
For  elliptic  configurations  either  an  equivalent  circular  cone 
value  should  be  used,  or  if  better  information  is  known,  the 
limiting  angles  associated  with  each  radial  traverse  should  be 
included  in  item  76. 

Item  numbers  57,  58,  and  59  specify  the  circular  store  shape 
and  are  required  only  for  MSHAPE (J)£50 .  They  are,  respectively: 
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Item  57: 


NSPOLJ  number  of  polynomials  specifying  store  shape; 

1  <  NSPOLJ  <  7 

Item  58: 

SXNDJ(J)  x/l  of  end  points  of  polynomials  specifying  store 

shape;  NSPOLJ  values 

Item  59 : 

SCOFJ(J,K)  coefficients  of  polynomials  specifying  shape 

These  data  specify  the  radius  distribution  of  the  store  and  are 
used  in  the  calculation  of  the  source-sink  distribution  which 
represents  the  store  volume  and  the  doublet  distribution  which 
represents  the  store  angle  of  attack  effects.  Up  to  seven 
polynomials  may  be  used.  The  polynomial  programmed  is  given  by 
Equation  (1) .  The  polynomials  must  be  input  for  a  shape  which 
is  made  dimensionless  by  the  store  length  since  the  trajectory 
program  is  written  assuming  this  to  be  the  case. 

Item  number  57  specifies  the  number  of  polynomials.  Item 
number  58  consists  of  one  card  which  contains  the  NSPOLJ  values 
of  the  end  points  of  the  polynomials  describing  the  shape. 

Item  number  59  is  a  set  of  NSPOLJ  cards  specifying  the  values 
of  the  coefficients  of  the  polynomials,  Equation  (1) .  All  seven 
coefficients  are  input  even  though  some  of  them  may  be  zero. 

For  circular  store  shapes,  this  is  the  last  item  required 
for  a  given  store  shape. 

Items  60  through  76  are  included  here  to  describe  the  store 
body  geometry  for  each  of  the  elliptic  store  shapes,  MSHAPE (J) >50. 
These  input  items  are  defined  similar  to  items  13  through  32  used 
to  describe  the  fuselage.  Some  of  the  same  input  subroutines  are 
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shared  between  the  fuselage  and  the  store.  The  generalizations 
to  nonelliptic  store  configurations  have  not  been  made. 

Item  60  consists  of  one  card  with  the  alphanumeric  descrip¬ 
tion  for  the  elliptic  store  shape  definition.  In  the  remaining 
cards  columns  73-80  in  items  61  through  76  have  been  left  open  for 
user  comments  and  will  be  used  here  to  designate  card  names. 

Item  61  defines  some  of  the  optional  print  and  execution 
parameters  associated  with  the  source  paneling  analysis  as 
follows. 

Item  61: 

IXZSYM  store  XZ  paneling  symmetry  option.  For  store  using 

paneling  option  both  sides  of  body  must  be  paneled 
due  to  nonuniform  flow  field.  For  the  store,  input 
IXZSYM=1 . 

IXZSYM=1,  configuration  is  symmetric;  panel  both 
halves  of  body  from  data  of  first  half 
IXZSYM=-1,  configuration  is  asymmetric;  data  for 
both  halves  are  input 

IPRT  optional  print  control  parameter 

IPRT ( 1 ) =0 ,  do  not  print  copy  of  input  data 
IPRT (1 ) =1 ,  print  formatted  copy  of  input  data 

IPRT (2) =0,  do  not  print  supplementary  geometry  data 
IPRT (2) =1,  print  panel  geometry,  angles,  and  areas 

IPRT (3) =0,  do  not  print  modified  shock  wave  shape 
IPRT (3) =1,  print  modified  shock  wave  shape  in  BSHOCK 

IPRT (4 ) =0 ,  not  used 

IPRT (5) =0,  do  not  print  velocity  and  source  strength 
output 

IPRT (5) =1,  print  source  strengths  and  velocities 
used  in  boundary  condition  associated 
with  each  panel 
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IUVW 


component  velocity  influence  coefficient  calculation 
option.  For  the  store  shapes,  input  IUVW=0. 

IUVW=0,  save  coefficient  arrays 
IUVW=1,  do  not  save  coefficient  arrays 

NSHOCK  controls  number  and  where  modified  shock  wave  shape 

traverses  are  made, 

NSHOCK=0,  a  single  traverse  is  calculated  at  0°, 

the  bottom  of  the  store;  circular  symmetry 
is  assumed 

NSHOCK=positive ,  NSHOCK  equally  spaced  shock 

traverses  are  computed  between  0°  and  90° 
NSHOCK=negative ,  NSHOCK  traverses  are  computed  at 
input  values  of  PHIS(K),  K=1,|NSH0CK|  in 
item  76 

MAXSHK  maximum  number  of  integration  steps  in  radial 

direction  used  to  generate  nonlinear  shock  shape; 

1  £  MAXSHK  <  100/NSHOCK 

NINLET  number  of  open  inlet  panels  (not  applicable  to 

stores) ;  input  NINLET=0 

NINBLK  number  of  blocked  inlet  panels  (not  applicable  to 

stores) ;  input  NINBLK=0 

N1NVEL  number  of  panels  used  in  inlet  velocity  calculation 

(not  applicable  to  stores) ;  input  NINVEL=0 

The  parameter  IXZSYM  is  to  allow  the  user  to  model  the  full 
body  from  the  geometry  of  only  the  right  hand  side,  or  a  full 
nonsymmetric  configuration.  Since  the  store  is  immersed  in  a 
nonuniform  flow  field,  both  halves  must  be  paneled.  The  print 
controls,  IPRT,  are  used  to  control  output  from  each  phase  during 
execution.  The  output  from  the  input  phase  is  used  to  echo  the 
data  read.  The  geometry  output  includes  the  panel  coordinates, 
control  points,  areas  and  inclination  angles.  The  shock  output 
includes  the  nose  shock  shape.  The  velocity  print  control  is  used 
only  to  examine  the  influence  coefficient  matrix.  The  last  print 
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control  displays  the  panel  strength  and  boundary  conditions. 

The  parameter  IUVW  is  used  to  suppress  the  storage  of  the 
component  velocity  coefficient  arrays.  These  are  only  required 
during  the  calculation  of  panel  pressures  and  loads.  This 
information  is  required  only  for  the  force  calculations  on  the 
separated  store.  Therefore  the  arrays  need  be  saved  only  for 
the  separated  store  shape. 

The  number  of  nose  shock  shape  traverses,  NSHOCK,  allows  the 
user  to  specify  the  traverses  which  best  describe  the  nonlinear 
shape  of  the  shock  in  the  meridional  direction  or  allow  the 
program  to  specify  an  even  spacing.  For  the  elliptic  store 
with  NSHOCK=0  the  nose  shock  shape  is  generated  at  $  =  0°,  the 
bottom  of  the  store,  and  assumes  both  right-left  and  top-bottom 
symmetry.  Meridional  traverses  are  required  between  <p  =  0°  and 
90°  so  NSHOCK  should  be  positive.  MAXSHK  limits  the  number  of 
points  in  each  radial  traverse  to  minimize  computing  time  and 
also  acts  as  the  internal  first  dimension  in  the  X  and  R  arrays 
containing  the  shape.  It  should  be  sufficient  to  allow  enough  points 
to  be  computed  to  describe  the  shape.  The  number  of  searches 
NSHOCK*MAXSHK,  used  to  define  the  shock  shape  is  among  the  most 
time  consuming  operations. 

The  quantities  NINLET,  NINBLK,  and  NINVEL  should  always  be 

zero. 


The  suggested  values  of  the  above  parameters  to  use  are 
IXZSYM=1,  IPRT (1 ) =IPRT ( 2 ) =IPRT { 3 ) =IPRT ( 5 ) =1 ,  IPRT(4)=0,  NSH0CK=4, 
and  MAXSHK=10.  The  parameter  IUVW  should  be  set  equal  to  zero 
for  the  separating  store.  It  may  be  set  to  one  for  any  shape  not 
to  be  separated  from  the  parent  aircraft  in  Program  II.  The  card 
is  designated  the  IOPTS  in  columns  73-80. 

Item  62  specifies  several  variables  which  locate  the  body 
interference  shell  and  modified  shock  wave  shape. 
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Item  62: 


XBIP  x-station  at  which  the  Y,Z  geometry  used  in  the 

definition  of  the  body  interference  shell  is 
obtained;  not  used  so  should  be  input  as  zero  or 
left  blank 

XSHLDR  x-station  which  is  considered  to  be  the  "shoulder" 

of  the  body  nose.  If  set  equal  to  zero,  the  value 
is  computed  from  the  location  of  the  first 
occurrence  of  a  maximum  in  the  cross-sectional  area 

EALPHA  angle  of  attack  correction  used  in  rotating  the 

body  shock  wave  shape  emanated  from  the  nose  about 
the  store  by  ac (1-EALPHA) 

Since  the  same  input  routine  is  being  used  for  both  the 
fuselage  and  stores,  XBIP  appears  here.  However,  no  value  for 
XBIP  is  required  for  the  stores  in  Program  I,  and  XBIP  should  be 
input  as  zero  or  left  blank.  Section  information  for  the  inter¬ 
ference  shells  on  the  separating  store  empennage (s)  are  generated 
from  input  in  the  second  program.  The  specification  of  the 
shoulder  x-location,  XSHLDR,  is  only  meant  to  provide  a  reasonable 
estimate  of  the  point  where  the  nose  stops  increasing  in  area. 

It  is  used  to  estimate  the  aft  point  where  the  modification  of 
the  shock  wave  shape  returns  to  the  original  linear  theory  value. 
The  parameter  EALPHA  defines  the  amount  of  deviation  from  a 
straight  rigid  rotation  of  the  store  nose  shock  wave  about  the 
nose.  The  angle  from  the  axis  of  symmetry  of  the  nose  shock  to 
the  body  centerline  is  thus  defined  as  ag  =  ac (1-EALPHA) .  EALPHA 
is  limited  to  values  from  0.0  to  1.0  representing  the  range  of 
rigid  shock  rotation  to  no  shock  rotation  with  angle  of  attack. 

The  value  of  EALPHA  for  circular  cones  may  be  determined  as  a 
function  of  Mach  number  and  nose  semi-vertex  angle  from  Table  14 
in  Reference  10.  A  plot  of  the  data  in  Table  14  is  also  found 
in  Figure  4  of  this  report.  In  lieu  of  better  information  for 
noncircular  shapes  use  the  vertex  angle  of  an  equivalent  circular 
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cone.  A  value  of  EALPHA=0.0  is  equivalent  to  the  methods 
assumed  by  the  circular  store  model.  This  card  is  designated 
XBIP  in  columns  73-80. 


The  next  eight  items,  numbers  63  through  70,  define  the 
external  shape  of  the  elliptic  store.  They  specify  the  store 
shape  in  as  great  a  detail  as  the  user  desires.  The  subdivision 
of  the  store  surface  described  here  into  discrete  panels  is 
performed  in  items  71  through  75. 

Item  63  contains  control  integers  used  to  specify  the  type 
and  amount  of  the  input  used  to  define  the  store  external  geometry. 
The  variables  defined  here  and  in  the  next  seven  items  control 
the  specification  of  the  external  body  shape  and  not  the  panel 
layout  as  follows. 


Item  63 ; 
JO 

J2 


read  reference  area  card;  for  stores  input  J0=0 
J0=0,  do  not  read  item  64 
J0=1,  yes,  read  item  64 


integer  flag  indicating  type  of  data  describing 
store  body  geometry  for  external,  shape. 

J2=-l,  data  for  circular  body  in  form  of  cross- 

sectional  areas  versus  XFUS  (item  65)  to  be 
entered;  include  section  data  in  item  67 

J2=-2,  data  for  circular  body  in  form  of  radius 

versus  XFUS  (item  65)  to  be  entered;  include 
section  data  in  item  68 

J2=-3,  data  for  elliptic  body  to  be  entered  in  form 
of  semi-axes  in  y-  and  z-directions;  include 
items  68  and  69 

J2=-4,  data  for  elliptic  body  to  be  entered  in  form 
of  vertical  semi-axis,  A,  and  elliptic 
ratio,  B/A;  include  items  69  and  70 
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J2=-5,  data  for  elliptic  body  to  be  entered  in 

form  of  horizontal  semi-axis,  B,  and  elliptic 
ratio,  A/B;  include  items  68  and  70 

J6  integer  indicating  whether  body  centerline  is 

c ambe red .  For  stores  input  J6=-l . 

J6=0,  body  is  cambered;  if  J2<0,  include  item  66 
J6=-l,  body  is  uncambered;  omit  item  66 

NFUS  number  of  body  segments;  for  store  only  one  body 

segment  is  allowed,  NFUS=1 

NRADX  number  of  points  used  to  specify  cross  section  of  the 

store.  If  the  configuration  is  symmetric  (IXZSYM=1) 
NRADX  is  input  for  the  half  section.  If  the  entire 
configuration  is  input  (IXZSYM=-1)  NRADX  is  input 
for  the  full  section.  If  the  body  is  circular  or 
elliptic,  the  program  computes  NRADX  Y-  and 
Z-ordinates  about  the  entire  section,  NRADX  33 

NFORX  number  of  axial  stations  on  body  at  which  external 

geometry  is  specified;  NFORX  <  51 

The  reference  area  read  under  J0=1  may  also  be  read  later  for 
K0=1  in  item  73.  For  the  store  to  be  separated,  reference  areas 
and  lengths  will  be  required.  The  K0  option  should  be  used  to 
input  tnese  items. 

The  quantities  NRADX  and  NFORX  specify  the  number  of  points 
in  the  meridional  and  axial  direction  to  be  input  or  computed  to 
define  the  external  store  shape.  The  value  of  NRADX  remains  con¬ 
stant  for  all  sections  along  the  store.  NRADX  represents  the 
number  of  pairs  of  Y  and  Z  values  calculated  at  each  section. 

NRADX  and  NFORX  should  be  greater  than  or  equal  to  KRADX  and 
KFORX  defined  later  to  avoid  paneling  that  lies  inside  the  actual 
external  geometry  shape.  Long  constant  geometry  sections  need 
only  be  defined  at  the  beginning  and  end.  When  used  to  define 
the  actual  panel  layout  they  should  follow  the  rules  for  KRADX 
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and  KFORX  in  defining  panels  later.  The  choice  of  NFORX  must 
further  satisfy  the  constraint  that  it  be  less  than  51.  This 
card  is  designated  JCARD  in  columns  73-80. 

Item  64  contains  the  reference  area  of  the  store.  It  may 
be  redefined  by  item  73.  This  item  is  not  included  if  JO  in 
item  63  is  input  as  J0=0.  If  not  defined,  REFA  is  set  equal  to 
1.0.  This  card  is  designated  REFA  in  columns  73-80. 

The  values  in  items  65  through  70  define  the  external  shape 
of  the  store  as  dictated  by  the  previous  control  variables.  This 
complete  set  of  items  is  included  only  once.  All  descriptions 
must  be  given  in  order  of  increasing  values  of  x  (nose  to  tail) . 

Items  65  and  66  specify  the  x-  and  z-values  of  each  axial 
station  defining  the  store  as  follows. 

Item  65: 

XFUS(I)  axial  station  at  which  body  cross  section  data  is 

included,  feet;  1=1,  NFORX 

Item  66; 

ZFUS(I)  input  only  if  J6=0  in  item  63;  z-value  of  camber 

offset,  feet;  1=1, NFORX 

The  values  of  XFUS  define  the  axial  stations  at  which  subsequent 
section  data  is  presented.  The  first  value  of  XFUS  (=0.0)  must 
be  at  the  body  nose.  All  values  must  increase  monotonically  to 
the  last  value,  XFUS (NFORX).  These  values  are  shown  pictorially 
in  Figure  6.  ZFUS  is  used  to  specify  curvature  in  the  centerline 
of  the  body.  The  cards  are  designated  the  XFUS  and  ZFUS  cards 
in  columns  73-80. 

For  J2  less  than  zero,  user  requested  combinations  of  items 
67  through  70  are  used  to  define  the  values  of  standard  sections 
as  follows  (1=1, NFORX). 
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Item  67 : 


FUSARD ( I )  for  J2=-l,  cross-sectional  area  of  circular  body, 

square  feet;  omit  this  item  for  other  values  of  J2 

I tem  6  8 ; 

FUSBY(I)  for  J2=-2,  radius  of  circular  body,  feet;  for 

J2=-3  or  -5,  length  of  horizontal  semi-axis,  B,  of 
an  ellipse,  feet;  omit  this  item  if  J2=-l  or  -4 

Item  69; 

FUSAZ(I)  for  J2=-3  or  -4  length  of  vertical  semi-axis,  A, 

of  an  ellipse,  feet;  omit  this  item  for  other  values 
of  J2 

Item  70: 

ERATIO(I)  for  J2=-4,  elliptic  ratio  B/A;  for  J2-5,  elliptic 
ratio  A/B;  omit  this  item  for  other  values  of  J2 

These  combinations  of  variables  define  the  possibilities  of 
specifying  the  circular  and  elliptic  shapes  at  each  axial  station. 
When  defining  an  elliptic  shape  (J2=-4  or  -5)  using  one  of  the 
semi-axes  and  the  elliptic  ratio,  only  the  first  value  has  to  be 
defined.  All  subsequent  values  equal  to  zero  (not  defined)  are 
set  equal  to  the  previous  value.  A  value  or  a  blank  must  be  set 
aside  for  each  value,  with  any  number  of  changes  in  elliptic 
ratio  in  between.  These  cards  are  designated  FUSARD,  FUSBY,  FUSAZ 
or  A/B  in  columns  73-80,  respectively. 

Item  71  contains  any  desired  identifying  information  relating 
to  the  paneling  sequence  from  the  previous  external  geometry  shape 

Item  72  defines  the  control  variables  used  to  specify  the 
revision  of  the  configuration  into  panels.  This  card  is  similar 
to  item  63  with  the  difference  that  it  defines  the  layout  of  the 
source  panels  themselves  as  follows. 
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Item  72: 

KO  read  reference  area  and  length  card; 

K0=0 ,  no  card  included 

K0=1,  yes,  include  item  73 
for  the  store  to  be  separated,  input  K0=1 

KRADX  number  of  meridian  lines  used  to  define  panel  edges. 

There  are  three  options  for  specifying  the  number 
of  panel  edges,  -33£KRADX<33 

KRADX=0 ,  the  number  of  meridians  is  set  equal  to 
NRADX 

KRADX=positive ,  the  meridian  lines  are  calculated 
at  KRADX  equally  spaced  angles  about  the 
body 

KRADX=negative ,  the  locations  of  meridian  lines 

are  input  at  specified  values  of  PHIK  in 
item  74 

KFORX  number  of  axial  stations  used  to  define  leading  and 

trailing  edges  of  panels.  Three  options  are 
available 

KFORX=0,  the  number  of  axial  stations  is  set  equal 
to  NFORX 

KFORX=positive ,  panel  edges  are  defined  at  KFORX 
values  of  XJ  defined  in  item  75 

KFORX=negative ,  NFORX  equally  spaced  panel  edges 
are  defined 

The  purpose  of  the  three  K-option  parameters  on  this  card  is 
to  specify  the  options  used  to  revise  the  external  body  shape 
defined  in  items  63  through  70.  This  allows  the  user  to  easily 
modify  the  panel  layout  without  changing  the  original  shape  input. 
For  the  store  to  be  ejected  the  reference  area  and  length  card 
must  be  read  (K0=1) .  These  dimensions  are  passed  to  Program  II 
and  are  the  values  used  in  normalizing  all  store  load  calculations. 
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For  symmetric  configurations  (IXZSYM=1),  KRADX  is  the  number 
of  meridians  on  the  half-section.  For  full  configurations 
(IXZSYM=-1) ,  KRADX  is  the  number  of  meridians  on  the  full  section 
including  the  meridians  at  0°  and  360°.  However,  for  IXZSYM=1 
after  generating  the  panel  geometry  KRADX  is  redefined  equal  to 
the  number  of  meridians  on  the  full  section  for  internal  program 
use.  The  value  of  KFORX  for  the  store  must  be  less  than  or 
equal  to  51.  The  number  of  panels  in  the  circumferential  and 
axial  directions  in  each  segment  is  thus  one  less  than  KRADX 
and  KFORX,  respectively.  For  elliptic  stores  KRADX  should  be  five 
or  greater  on  a  half  body  for  practical  results.  In  choosing  the 
number  of  axial  stations  for  a  store  body  the  resolution  of  the 
load  distribution  axially  must  be  of  concern.  If  the  store  is 
in  the  free  stream  or  is  not  to  be  separated,  a  ring  length  to 
diameter  ratio  as  high  as  one  should  be  sufficient.  If  the 
store  is  in  the  presence  of  shocks  propagating  from  other  bodies 
including  its  own  reflected  shock,  the  axial  spacing  should  be 
sufficiently  small  to  resolve  the  pressure  spike  of  the  shock  with 
at  least  two  to  three  rings  of  panels,  considering  both  the 
"smearing"  of  the  shock  pressure  rise  and  the  inclination  of  the 
shock  to  the  body.  This  ultimately  should  be  arrived  at  by  con¬ 
vergence  of  the  load  distribution.  The  total  number  of  panels  used 
around  the  body  interference  shell,  NBDCR,  is  also  derived  from 
KRADX  as:  NBDCR=KRADX-1 .  The  program  is  thus  restricted  to  a 
maximum  number  of  source  panels  NBODY= (KRADX- 1) (KFORX-1) £1600  and 
indirectly  by  the  number  of  panels  on  the  interference  shell 
NBDCR*NCWB£l00  in  the  second  program.  This  card  is  designated 
the  KCRAD  in  the  data  deck  in  columns  73-80. 

Item  73  defines  the  reference  lengths  and  areas  used  in 
force  and  moment  coefficient  definitions.  Omit  this  item  if 
K0  in  item  72  is  K0=0. 

Item  73: 

REFAR  reference  area,  square  feet.  If  greater  than  zero, 

redefines  REFA  in  item  64 
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REFD  reference  length  or  diameter  used  to  nondimen- 

sionalize,  feet.  If  blank,  set  equal  to  1.0 

REFL  body  length,  feet.  If  blank,  set  equal  to 

XFUS (NFORX) -XFUS (1 ) 


REFX 

Xg  -location 

of 

moment 

center. 

feet; 

see  Figure 

6 

REFZ 

z£  -location 

of 

moment 

center , 

feet; 

see  Figure 

6; 

input  REFZ=0 . 0 

Item  74  defines  the  meridian  angles  for  the  edges  of  panels 
as  optionally  requested  by  KRADX .  Omit  this  item  if  KRADXj>0  in 
item  72. 

Item  74: 

PHIK(I)  meridian  angle  of  panel  edges  expressed  in  degrees. 

Included  only  if  KRADX  is  negative;  1=1, | KRADX | 

The  convention  is  observed  that  PHIK=0°  is  at  the  bottom  of  the 
body  and  PHIK=180°  is  at  the  top  of  the  body.  This  option  should 
be  used  to  ensure  that  panel  edges  in  th£  presence  of  fins  on 
either  of  two  sets  of  empennages  match  exactly  the  meridian  angle 
at  the  fin-body  junction (s).  These  angles  are  also  used  in 
Program  II  to  lay  out  the  empennage  body  interference  panels. 

This  card  is  designated  PHIK  in  columns  73-80. 

Item  75  defines  the  axial  stations  for  the  edges  of  panels 
as  follows.  Omit  this  item  if  KFORX£0  in  item  72. 

Item  75; 

XJ ( I )  x-stations  along  body  defining  panel  edge  between- 

rings,  feet;  I=1,KF0RX 

These  values  are  used  to  define  both  the  leading  and  trailing 
edges  of  adjacent  rings  of  panels.  If  KFORX  is  zero,  these 
stations  are  set  equal  to  XFUS  in  item  65.  When  selecting  values 
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in  the  vicinity  of  the  body  interference  panels,  the  XJ  should 
be  specified  which  corresponds  exactly  to  the  leading  and 
trailing  edge  of  the  empennage-body  interference  shell.  This 
card  is  designated  XJ  in  columns  73-80. 

Item  76  defines  the  location  of  the  traverses  which  form 
the  table  describing  the  modified  shock  wave  shape.  This  item  is 
omitted  if  NSHOCK>0  in  item  61. 


PHIS (I) 


specified  angles  at  which  integration  of  modified 
shock  wave  shape  is  calculated,  degrees; 

1=1, |nshock | 


The  convention  is  used  such  that  PHIS(l)  is  zero  at  the  centerline 
on  the  bottom,  and  increases  counterclockwise  viewed  from  the 
rear.  For  those  stores  which  are  not  separated,  PHIS(NSHOCK) 
need  only  span  the  region  in  which  the  separated  store  may  travel. 
For  the  store  to  be  separated,  the  PHIS  values  can  be  specified 
from  0  to  90  degrees,  which  by  symmetry  will  provide  a  description 
of  its  shock  throughout  360  degrees.  If  NSHOCK  >  0,  the  program 
generates  values  from  0°  to  the  angle  formed  by  a  radial  line 
through  the  wing  tip.  Only  values  which  highlight  nonlinearities 
in  the  Rg^oc^  versus  4>  need  be  defined.  For  the  store  all 
modified  shocks  are  calculated  at  zero  degree  angle  of  attack  and 
with  PHIS (NSHOCK) =90 0 .  The  resulting  table  interpolation  assumes 
symmetry  in  the  shock  wave  shape  in  both  the  XY-  and  XZ-planes. 


3.2.3  Sample  input  data 


A  sample  input  data  deck  will  now  be  presented.  It  will 
utilize  the  wind  tunnel  model  components  shown  in  Figure  9  except 
that  they  have  been  scaled  up  by  a  factor  of  20  and  converted 
from  inches  to  feet  in  order  to  approximate  a  full-scale  aircraft 
and  store.  The  models  were  used  in  the  wind  tunnel  test  program 
described  in  Reference  11. 


The  configuration  for  the  sample  calculation  is  shown  in 
Figure  10.  It  consists  of  a  noncircular  wing-fuselage-inlet 
combination  of  Figure  9(a),  the  pylon  of  Figure  9(d),  and  the 
elliptic  store  of  Figure  9(e).  The  store  is  shown  in  its  initial 
position  at  the  carriage  position  on  the  pylon. 

The  input  data  deck  is  tabulated  in  Figure  11.  The  item  num¬ 
bers  indicated  on  the  figure  correspond  to  those  of  Figure  2. 

The  first  item  on  Figure  11  is  item  number  1  which  contains  the 
value  of  NCARDS,  in  this  case  8.  This  is  followed  by  item  number 
2,  which  consists  of  these  8  cards  of  identifying  information. 

Item  number  3  specifies  the  aircraft  flight  conditions.  Only 
the  angle  of  attack  and  Mach  number  of  the  parent  configuration 
are  required  in  the  first  program.  The  angle  of  attack  is  5.0° 
and  the  Mach  number  is  1.5. 

The  aircraft  components  which  are  present  are  specified  by 
item  4.  There  is  a  noncircular  fuselage,  NFU=2;  there  is  a 
pylon,  NPY=1;  and  there  is  one  elliptic  store,  NSTRS=1. 

Items  11  through  32  define  the  noncircular  fuselage  input 
data.  Item  11  contains  the  fuselage  length,  maximum  fuselage 
radius,  maximum  nose  shock  angle,  and  body  interference  shell 
length.  Item  12  contains  the  number  of  rings  of  u-velocity 
panels  on  the  interference  shell.  The  remaining  items,  13 
through  32,  define  the  fuselage  external  geometry  and  source 
panel  layout.  The  body  length,  maximum  radius,  and  interference 
shell  length  are  shown  in  Figure  10.  The  maximum  nose  shock 
angle  is  determined  from  Chart  5  of  Reference  9  for  a  nose  cone 
semivertex  angle  of  16°.  Since  the  wing  trailing  edge  is  super¬ 
sonic  and  since  no  rearward  motion  of  the  store  is  expected. 


77 


BODYPL  is  input  as  the  length  of  the  wing  root  chord  as  shown 
in  Figure  10.  Nine  rings  of  panels  are  placed  on  the  interference 
shell,  NCWB=9 ,  to  correspond  to  the  number  of  rows  of  panels  on  the 
wing  chord.  The  number  and  spacing  of  panels  around  the  circum¬ 
ference  is  set  identical  to  the  Y,Z  values  computed  at  station  XBIP 
indicated  in  item  15. 

Item  13  through  21  define  the  external  shape  of  the  fuselage 
using  the  arbitrary  cross  section  option.  Item  13  is  the 
descriptive  title  for  the  fuselage  model  configuration.  Item  14 
defines  optional  print  and  execution  parameters  controlling 
the  source  panel  calculations.  The  fuselage  layout  is 
defined  symmetric,  IXZSYM=0;  the  print  parameters  IPRT(1)=1, 

IPRT (2) =1 ,  IPRT(3) =1,  IPRT (4 ) =0 ,  and  IPRT(5)=1  control  the 
printing  of  input  data,  panel  geometries,  nose  and  inlet  shock 
shapes,  and  panel  strengths  and  supression  of  print  of  velocity 
coefficient  calculations.  The  parameter  IUVW=1  suppresses  the 
writing  of  source  panel  influence  coefficient  arrays  onto  TAPE12. 
The  parameter  NSHOCK«3  generates  equally  spaced  traverses  defining 
the  nonlinear  nose  shock  shapes  at  PHIS=0°,  45°,  and  90°.  MAXSHK=8 
limits  the  number  of  radial  steps  in  the  shock  traverses  to  eight. 
The  inlet  parameters,  NINLET,  NINBLK,  and  NINVEL,  identify  that 
there  are  4  open  inlet  panels,  2  blocked  inlet  panels  and  19 
additional  panels  to  be  used  in  the  velocity  calculations  used  in 
the  definitions  of  the  inlet  shock  shape.  This  panel  layout  on 
the  inlets  is  specified  in  order  to  match  an  inlet  mass  flow  ratio 
of  0.87. 


Item  15,  the  XBIP  card,  defines  the  axial  station  at  which 
the  body  interference  shell  cross  section  is  defined,  the  x- 
location  of  the  fuselage  shoulder  and  the  shock  angle  of  attack 
correction.  The  interference  shell  geometry  is  generated  at 
XBIP=28.0  as  shown  in  Figure  10.  XSHLDR  is  set  equal  to  zero  to 
allow  the  program  to  compute  the  value  from  the  first  nonexpanding 
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area  location.  The  value  EALPHA=0.315  is  obtained  from  Figure  4 
for  a  semivertex  cone  angle,  c,  of  16°  and  Mach  number  of  1.5. 

Item  16,  the  JCARD,  defines  the  form  and  number  of  external 
geometry  values  to  be  read  for  the  noncircular  configuration. 

J0=0  indicates  that  no  reference  area  (item  17)  will  be  read. 

J2=l  specifies  that  the  arbitrary  Y-Z  cross  section  option  will 
be  used  to  define  the  body  shape.  Any  chamber  in  the  body  will 
be  handled  through  the  YZ  input  so  J6=-l.  NFUS=5  identifies 
that  there  will  be  five  body  segments  used  to  define  the  geometry. 
The  five  pairs  of  KRADX(I)  and  KFORX(I)  define  the  number  of 
meridional  and  axial  stations  to  be  input  for  each  body  segment. 

Items  18,  20,  and  21  are  the  axial  stations  of  the  body 
cross  sections,  and  the  Y  and  Z  coordinates  of  the  meridional 
points  defining  the  fuselage  external  shape.  Item  18,  the  XFUS 
cards,  specify  the  19  axial  stations  in  the  first  segment  at 
which  Y-Z  pairs  are  input.  They  are  followed  by  19  pairs  of  Y,Z 
cards  with  10  values  of  Y  or  Z  per  card.  They  are  identified 
by  Y-N3B2WA4 ,  Z-N3B2WA4  for  the  wind  tunnel  model  configuration 
used.  The  next  three  repetitions  of  items  18,  20,  and  21  define 
the  cross  sections  at  the  inlet,  while  the  last  set  defines  the 
fuselage  aft  of  the  inlet.  The  values  of  Y  and  Z  are  computed 
from  the  values  of  R  and  X  in  Figure  9 (b)  at  each  of  the  fuselage 
stations.  The  cross-sectional  shape  is  subdivided  to  get  a  total 
number  of  nine  panels  on  the  half  body:  two  on  the  bottom 
straight,  one  on  the  lower  rounded  corner,  two  on  the  vertical 
side,  and  four  on  the  circular  arc  on  the  top.  If  KRADX(I)  and 
KFORX(I)  in  item  27  are  set  to  zero,  the  input  values  of  Y  and  Z 
will  become  panel  corners.  The  sets  of  values  for  the  inlet  are 
input  at  the  exact  corners  to  be  used  for  the  inlet  panels.  The 
Y  and  Z  values  proceed  counterclockwise  around  each  station.  Where 
a  triangular  panel  is  formed,  identical  Y  and  Z  values  are  input 
twice . 


Items  26  through  32  define  the  subdivision  of  the  external 
geometric  shape  of  the  fuselage  into  panels.  Item  26  is  the 
descriptive  title  containing  additional  information  particular  to 
the  panel  layout.  Item  27,  designated  the  KCARD,  contains  the 
reference  area  option  and  the  number  of  divisions  to  be  used  to 
partition  the  external  geometry  into  panels.  The  option  K0=0 
designates  that  no  reference  length  card,  item  30,  will  be  read. 

It  is  followed  by  five  sets  of  KRADX(I)  and  KFORX(I),  one  for 
each  segment.  In  the  pair  for  the  first  segment,  KRADX(1)=0  and 
KFORX (1) =7  specify  the  same  meridional  spacings  as  the  external 
geometric  shape,  while  seven  axial  sections  are  used  to  partition 
the  fuselage  forward  of  the  inlets  into  six  rings  of  panels.  The 
next  three  pairs  of  KRADX  and  KFORX  define  the  inlet  section 
panel  rings  to  use  the  same  coordinates  input  in  the  external 
shape.  The  sample  case  shows,  for  these  three  segments,  KRADX=0 
and  KFORX=2.  KFORX  could  also  be  input  as  KFORX=0  in  which  case 
the  second,  third,  and  fourth  cards  in  item  32  in  Figure  11  would 
be  omitted.  With  KF0RX=2  the  same  axial  layout  is  specified  as 
was  input  in  items  18,  20,  and  21.  The  last  segment  is  partitioned 
into  four  axial  sections  and  uses  the  same  coordinates  as  the 
input  external  shape  to  define  the  panel  edges.  The  panel  layout 
is  shown  pictorially  for  the  fuselage  in  Figure  12.  A  three- 
dimension  isometric  view  of  the  fuselage  showning  the  body 
segments  and  inlets  is  presented  on  the  first  page,  Figure  12(a). 

Items  28  and  29  define  the  properties  of  the  inlet  including 
the  number  of  traverses  used  to  define  the  inlet  shock  shape, 
the  inlet  panel  numbers,  the  panel  numbers  used  in  the  shock  cal¬ 
culation,  and  the  special  inlet  properties.  The  first  value  in 
item  28,  the  INLETA  card,  NIS=3  specifies  that  three  radial 
traverses  will  be  used  to  define  the  inlet  shock  shape.  The  next 
six  values  on  the  same  card  are  the  panel  numbers  of  the  source 
panels  modeling  the  surface  of  the  inlet  face.  The  first  four 
panel  numbers,  60,  61,  74,  and  75,  designate  the  open  inlet  panels. 
Similarly,  the  next  two  panel  numbers,  88  and  89,  designate  the 
blocked  inlet  panels.  The  layout  of  these  panels  is  depicted  in 
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the  fuselage  isometric  projection  in  Figure  12  (a)  and  the  cor¬ 
responding  rear  view  in  Figure  12 (b) .  The  remaining  values  on 
the  cards  of  item  28  are  the  values  of  JINLT  which  designate  the 
remaining  panel  numbers  of  the  set  of  panels  to  be  used  in  the 
inlet  velocity  computations  for  the  shock  shape.  These  additional 
inlet  panels  are  indicated  on  the  inlet  rear  and  side  views  in 
Figures  12(b)  and  12(c),  respectively.  The  value  BTINLT=.58  in 
item  29  is  selected  to  be  about  3-4%  less  than  the  cotangent  of 
the  angle  of  the  panel  to  the  flow  (cot  59°=. 60). 

The  five  cards  of  item  32,  designated  XJ  cards,  define  the 
axial  stations  bounding  the  panel  rings  in  each  of  the  five 
segments.  Each  card  contains  the  axial  stations  dividing  the 
rings  of  panels.  The  last  value  in  each  segment  must  also  be  the 
first  value  in  the  following  segment.  The  fuselage  is  only 
paneled  aft  to  body  station  45.845. 

The  next  eight  items,  items  34  through  38  and  40  through  42, 
are  the  wing  input  data.  Item  34  gives  the  position  of  the  wing 
root  chord  leading  edge  relative  to  the  fuselage  and  the  wing 
incidence  angle  and  item  35  specifies  the  root  chord  length  and 
the  semispan.  These  quantities  are  shown  in  Figure  10.  Items 
36  and  37  are  data  required  by  the  program  to  lay  out  the  constant 
u-velocity  panels.  There  are  to  be  9  panels  in  each  chordwise 
row  and  8  of  these  rows  across  the  semispan.  This  requires  that 
the  spanwise  locations  on  the  left  wing  panel  of  nine  side  edges 
and  the  sweep  angles  and  dihedral  angle  to  the  right  of  these 
points  be  specified.  These  data  are  contained  on  the  nine  cards  of 
item  37.  Note  that  the  first  side  edge  coincides  with  the  wing- 
fuselage  juncture.  Item  38  indicates  that  the  wing  has  neither 
twist  nor  camber,  and  thus  item  39  is  omitted  from  the  input  data 
deck.  Items  40,  41,  and  42  specify  the  wing  thickness  distribution. 
The  three  indices  on  item  40  indicate  that  50  thickness  panels 
are  to  be  placed  in  one  chordwise  row  (MSWS=1)  and  that  the  wing 
has  a  similar  thickness  distribution  at  all  spanwise  stations 
(NUNIS=1) .  The  single  chordwise  row  of  thickness  panels  is 
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bounded  by  the  two  spanwise  stations  with  leading  and  trailing 
edge  sweep  angles  and  dihedral  angle  given  in  item  41.  Item  42 
consists  of  seven  cards  with  the  50  values  of  the  slope  of  the 
thickness  distribution.  The  airfoil  section  specification  is 
shown  in  Figure  9 (a) .  The  wing  leading  edge  in  this  case  is 
subsonic.  The  maximum  value  of  the  slope  has  been  selected  in 
accordance  with  the  discussion  in  Section  3.4  of  Reference  2. 

Since  there  is  a  pylon,  items  43  through  48  are  input  and 
contain  the  pylon  data.  The  pylon  is  the  one  shown  in  Figure  9(d) , 
except  that  it  has  been  scaled  up  by  a  factor  of  20.  The  pylon 
as  shown  in  Figure  9(d)  extends  through  the  hole  in  the  body  duct 
assembly  in  Figure  9(c)  to  attach  directly  to  the  body.  The 
height  dimension  here  is  the  exposed  portion  of  the  pylon.  Item 
43  of  the  input  data  specifies  that  the  pylon  is  located  below 
the  fuselage  centerline  (IP=0) ,  the  leading-edge  and  trailing- 
edge  sweep  angles  are  0.0°,  the  root  chord  length  is  4.44  feet, 
the  exposed  height  is  0.75  feet,  and  the  leading  edge  is  4.53  feet 
behind  the  local  wing  chord  leading  edge  defined  to  be  the  wing 
root  chord  for  IP=0  (see  Figure  10) .  Items  44  and  45  are  the 
pylon  constant  u-velocity  panel  data.  There  are  to  be  two  panels 
in  a  chordwise  row  and  two  of  these  rows  spanwise.  The  three  side 
edge  locations  are  given  by  item  45.  Items  46  and  47  specify  the 
pylon  thickness  distribution.  There  are  to  be  40  thickness  panels 
in  one  chordwise  row  and  the  thickness  distribution  is  similar  at 
all  spanwise  stations.  Item  48  contains  the  slopes  of  the  panels 
at  the  centroid  of  each  of  the  40  thickness  panels. 

No  rack  is  present  so  items  49  through  53  are  omitted. 

Item  54  begins  the  store  data.  The  store  is  shown  in 
Figure  9 (e)  but  has  been  scaled  up  by  a  factor  of  20  and  the 
dimensions  converted  from  inches  to  feet.  Item  54  contains  the 
number  assigned  to  the  store  and  its  shape  number.  The  value  of 
NSHAPE ( 1 )  greater  than  50  indicates  that  this  store  is  to  be 
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modeled  using  the  elliptic  store  option.  The  next  seven  quantities 
in  this  item  are  the  store  length,  maximum  radius,  location  of  the 
store  nose,  and  attachment  incidence  and  roll  angles.  This  informa¬ 
tion  is  shown  in  Figure  10. 

Item  number  55  contains  the  value  of  NSHPT .  Only  one  store 
shape  is  input  here.  Items  56  through  75  contain  the  shape  data 
for  the  store.  Item  56  contains  the  value  of  MS H APE ,  which  must 
be  equal  to  the  value  of  NSHAPE(l)  read  in  as  part  of  item  54. 

The  value  of  MSOR  is  MSOR=0  since  the  elliptic  store  option  is 
being  used.  It  also  specifies  the  store  limiting  nose  shock  angle 
STHSHK=54°.  For  the  elliptic  store,  the  value  STHSHK  is  determined 
from  Chart  5  of  Reference  9  for  an  area  equivalent  nose  semivertex 
angle  of  24.4°. 

Items  60  through  75  define  the  store  using  the  same  methods 
and  routines  used  to  specify  the  fuselage  panel  layout.  The 
first  card,  item  60,  is  the  descriptive  title  for  the  store  shape. 
Item  61  is  the  same  IOPTS  card  as  in  item  14.  The  first  value, 
IXZSYM=1,  specifies  the  store  configuration  as  symmetric  with  both 
halves  to  be  paneled  with  only  the  positive-y  side  of  the  shape 
input.  The  print  control,  IPRT,  prints  the  input,  layout  geometry, 
shock  shape,  and  source  distribution  and  the  velocity  coefficient 
data.  IUVW=0  requests  that  the  u,v,w  velocity  coefficients  be 
saved  for  transmission  to  Program  II.  The  value  NSHOCK=3  specifies 
that  three  radial  traverses  be  used  to  define  the  store  shock 
shape  meridionally  at  4)  =  0°,  45°,  and  90°.  MAXSHK=8  limits  the 
number  of  radial  steps  in  each  traverse  to  eight.  No  inlet  panels 
are  included  or  permitted. 

Item  62  is  the  XBIP  card.  No  value  is  specified  for  XBIP  in 
Program  I.  All  values  for  the  interference  shell  will  be  computed 
in  Program  II.  The  shoulder  location  XSHLDR=4 . 167  is  used  to 
define  the  aft  end  of  the  store  nose.  The  value  EALPHA=0 . 575  is 
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obtained 
the  area 
elliptic 
an  early 


from  Figure  4  by  interpolation  of  the  M=1.5  line  to 
equivalent  nose  semi vertex  cone  angle  of  24.4°  for  the 
store.  This  value  of  EALPHA=0.575  was  obtained  from 
version  of  Figure  4 .  A  more  correct  value  would  be 


0.61. 


Items  63  through  70  define  the  external  shape  of  the  elliptic 
store.  Item  63,  the  JCARD,  defines  the  options  used  or  read  in 
specifying  the  store  shape.  The  first  value,  J0=0,  requests  that 
no  reference  area  card  (item  64)  be  read.  J2=-5  specifies  that 
the  elliptic  body  option  requesting  the  horizontal  semi-axis,  BY, 
and  the  ratio,  A/B  be  read.  The  value  J6=-l  identifies  that 
top-bottom  symmetry  exists  and  that  no  Z-value,  item  66,  be  read. 
The  value  NFUS=1  specifies  that  only  one  segment  (a  program 
restriction)  will  be  used  to  define  the  store  body.  NRADX=19 
and  NFORX=23  subdivide  the  surface  geometry  into  18  meridional 
divisions  on  the  half  body  and  uses  23  axial  body  sections  to 
define  the  external  store  shape.  The  meridional  angles  on  the 
right  half  of  the  store  are  mirror  imaged  to  the  left  half  to 
model  the  full  configuration.  The  meridional  angles  are  computed 
internally  at  equally  spaced  intervals  to  generate  the  Y,Z 
coordiantes  on  the  elliptic  shape.  Twenty-three  values  of  X,  BY, 
and  A/B  are  input  in  items  65,  68,  and  70. 

Item  71  contains  the  title  associated  with  the  store  paneling. 
Item  72,  the  store  KCARD,  has  three  indices.  K0=1,  requests  that 
the  reference  area  and  length  card  (item  73)  be  read  since  forces 
and  moments  will  be  computed  later  in  Program  II.  KRADX=-5 
indicates  from  the  negative  sign  that  the  new  meridional  angles 
will  bt  read,  and  from  the  magnitude  that  five  values  will  be 
input  to  define  the  right  half  of  the  store  body.  The  value, 
KFORX=19 ,  specifies  that  19  x-stations  will  be  used  to  define  18 
rings  of  eight  panels. 
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Item  73  contains  the  reference  area  and  lengths  used  in  the 

force  and  moment  calculations.  REFAR  uses  the  maximum  cross- 

2 

sectional  area  of  2.4544  ft  .  The  reference  length,  REFD=1.7678  ft, 
is  obtained  for  an  equivalent  circular  body  diameter.  The  third 
value,  REFL=10.0,  is  the  store  body  length.  It  is  followed  by 
the  x- location  of  the  moment  center,  REFX=5.0,  set  to  the  center 
of  the  store.  The  last  value,  REFZ=0.0,  indicates  that  there  is 
no  z-offset  of  the  moment  center. 

Item  74  contains  the  meridional  angles  used  to  subdivide  the 
store  into  panels  about  the  circumference.  Five  values  are  input 
for  the  store  half  body.  The  panel  layout  defines  panel  edges 
at  merdians  corresponding  to  the  angles  of  attachment  of  the 
interdigitated  tail  fins.  The  first  and  last  values,  PHI=0°  and 
180°,  must  lie  on  the  XZ  plane  of  symmetry.  The  last  item,  the 
two  XJ  cards,  contain  the  19  values  for  the  axial  stations  sub- 
dividing  the  panel  layout.  This  number  was  chosen  to  provide 
sufficient  axial  resolution  of  the  load  distribution  in  the 
presence  of  reflected  shocks.  It  is  noted  that  values  XJ=5.833, 
8.333,  8.75,  and  10.0  correspond  to  leading  and  trailing  edges 
of  the  two  sets  of  fins.  These  values  locate  the  leading  and 
trailing  edges  of  the  two  empennages  to  be  defined  in  Program  n 
for  the  ejected  store  as  depicted  in  Figure  13. 

3.3  Description  of  Output  from  Program  I 

Figure  14  presents  the  output  from  the  computer  program  for 
the  sample  case,  the  configurations  of  Figure  10,  and  the  data 
deck  of  Figure  11. 

The  first  page  of  output.  Figure  14(a),  tabulates  the  input 
identifying  information,  the  aircraft  flight  conditions,  and  the 
fuselage  data  input  in  item  11.  The  remaining  fuselage  data 
describing  the  noncircular  fuselage  including  inlets  is  continued 
over  the  next  five  pages.  It  is  followed  by  the  panel  geometry, 
singularity  strengths  and  nose  and  inlet  shock  shapes. 


Figures  14<b) ,  14  (c) ,  14  (d) ,  and  14 (e)  repeat  the  input  data 
for  the  external  geometric  shape  data  for  the  scaled  version  of 
experimental  model  configuration  N3-B2-A4 .  It  starts  with  the 
program  options  and  is  followed  by  three  pages  of  Y,Z  coordinates 
for  each  of  the  fuselage  cross  sections  in  each  of  the  five  body 
segments.  At  the  end  of  the  Y,Z  input  are  the  computed  maximum 
cross  sectional  area,  shoulder  area  and  x-locations,  and  body 
length  between  the  first  and  last  input  stations.  If  XSHLDR  is 
input  the  shoulder  area  is  not  computed.  With  these  exceptions, 
all  of  the  above  variables  are  input  quantities  whose  descriptions 
are  found  under  the  definitions  of  the  input  variables.  All 
fuselage  data  printed  in  this  section  are  controlled  by  print 
control  IPRT(1). 

Figure  14(f)  repeats  the  input  for  the  control  parameters  and 
variables  which  describe  the  layout  of  the  source  panels  on  the 
noncircular  fuselage.  They  include  the  special  properties  of  the 
inlet  panels  and  the  x-stations  at  which  Y,Z  data  is  interpolated 
in  the  external  geometric  shape  above. 

Figures  14(g)  through  14(1)  contain  the  output  of  the  source 
panel  corner  points,  control  points,  orientation  angles  and  panel 
areas  as  computed  for  each  of  the  Individual  panels.  This  output 
is  controlled  by  IPRT(2)  and  represents  the  primary  panel  property 
data  passed  to  the  second  program.  It  typically  may  be  helpful 
in  checking  geometric  input  and  in  Identifying  panel  indices  for 
the  inlet  panel  input. 

The  first  part  of  Figure  14 (m)  contains  the  information  for 
the  fuselage  interference  shell  constant  u-velocity  panels.  It 
includes  the  number  and  location  of  the  panels  in  wing  coordinates 
and  the  table  of  Y, Z  coordinates  obtained  from  the  source  panel 
data  to  be  used  as  the  noncircular  cross  section  defining  the 
interference  shell.  This  section  is  generated  from  the  external 
geometry  at  axial  station  XBIP.  The  remainder  of  Figure  14 (m) , 
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and  Figures  14 (n)  and  14  (o)  tabulate  the  source  panel  singularity 
strengths  and  the  boundary  conditions  used  in  their  confutation . 

Figures  14 (p)  and  14 (g)  tabulate  the  nonlinear  shock  shapes 
computed  for  the  fuselage  nose  and  for  the  ramp  inlet.  These 
tables  form  a  two-dimensional  representation  of  x- location  of  the 
shocks  as  functions  of  R  and  a.  Tabulated  along  each  radial 
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traverse  in  the  shape  integration  are  the  axial  location#  X#  for 
the  prescribed  distance  R#  the  axial  velocity  along  the  wind  axis# 
U#  maximum  radial  component  of  velocity#  W#  the  confuted  local 
Mach  number#  M^,  the  value  of  6  -  /M^  -  1,  the  local  shock  slope# 
DX/DR#  the  increment  in  Prandtl-Meyer  turning  angle  due  to  the 
local  flow#  DNU#  and  the  total  velocity  in  the  wind  axis#  UVM8Q. 

Figures  14 (r) ,  14 (s) ,  and  14 (t)  tabulate  most  of  the  wing 
and  pylon  input  data.  Figure  14 (r)  tabulates  the  wing  data  exclu¬ 
sive  of  the  twist  and  camber  and  thickness  distributions. 

Figure  14 (s)  tabulates  the  pylon  data  except  for  the  thickness 
distribution.  The  input  thickness  distributions  for  the  wing 
and  pylon  are  tabulated  in  Figure  14 (t) . 

Figures  14 (u)  through  14 (ii)  tabulate  the  input  and  confuted 
data  for  the  elliptic  store.  Figure  14 (u)  contains  the  store 
input  numbers#  length#  radius#  location#  and  orientation  relative 
to  the  local  wing  chord.  For  the  store  located  on  the  centerline# 
these  locations  are  relative  to  the  projection  of  the  root  chord 
through  the  body. 

Figure  14 (v)  contains  the  input  control  parameters  and 
elliptic  store  external  shape  variables.  As  requested  by  the 
option  J2— 5 ,  the  distributions  of  horisontal  semi-axis  and 
elliptic  ratio  are  input  versus  axial  station.  Since  the  elliptic 
ratio  if  constant  for  the  entire  store#  only  the  first  value  has 
to  be  input.  The  computed  vertical  semi-axis  is  also  shewn. 

Figure  14 (w)  contains  the  tabulated  values  input  for  the  control 
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parameters  which  are  used  with  the  input  geometry  to  divide  the 
store  into  discrete  panels.  The  reference  lengths  and  areas, 
the  meridional  partitioning  of  circumference  and  the  axial 
stations  subdividing  the  rings  are  shown. 

Figures  14 (x)  through  14 (cc)  contain  the  store  source  panel 
corner  points  and  computed  geometric  properties  as  controlled 
by  the  store  print  control  IPKT(2).  Two  singularity  strength 
solutions  are  computed  when  the  store  is  at  angle  of  attack.  The 
two  are  given  in  Figures  14  (cc)  through  14  (ii) .  The  singularity 
strengths  and  boundary  conditions  for  the  store  at  ac  ■  0*  are 
shown  in  Figures  14 (cc)  through  14  (ff ) .  The  shock  shape  tabulated 
for  the  store  at  zero  degrees  angle  of  attack  is  presented  in 
Figure  14  (ff)  for  the  three  radial  traverses.  This  shock  shape 
is  rotated  about  the  store  later  as  needed  for  use  at  angle 
of  attack.  Second,  the  store  singularity  strengths  computed  at 
the  input  angle  of  attack  are  shown  in  Figures  14 (gg)  through 
14 (ii) . 

The  next  three  pages  of  output.  Figures  14 (jj)  through  14(11), 
tabulate  quantities  associated  with  the  constant  u- velocity  panel 
layout  on  the  wing,  pylon,  and  fuselage  and  the  boundary  condi¬ 
tion  at  the  control  points  of  these  panels.  The  x,y,z  coordinates 
are  those  of  the  control  points  in  the  wing  coordinate  system, 
see  Figure  3.  The  next  three  columns,  U/VINP,  V/VINF,  and  W/VINF, 
are  the  sums  of  the  dimensionless  perturbation  velocities  in  the 
xw,  yw,  and  zw  directions,  respectively,  induced  at  the  control 

points  by  the  other  aircraft  components.  These  include  fuselage 
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on  the  wing  and  pylon,  wing  thickness  on  pylon,  wing  thickness 
on  wing  when  wing  has  dihedral,  pylon  thickness  on  wing,  and 
both  wing  and  pylon  thickness  on  the  fuselage.  VXNF  is  the  free- 
stream  velocity.  The  next  to  last  column  tabulates  the  twist 
and  camber  distribution  input  for  the  wing.  The  last  column 
tabulates  the  singularity  strengths  of  the  constant  u-velocity 
panels  laid  out  on  the  wing,  pylon,  and  fuselage. 
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,  3.4  External  Dataset  Generated  by  Program  I 

At  the  end  of  the  execution  of  the  first  program  to  compute 
the  parent  aircraft  flow  model,  the  data  required  to  be  passed 
to  the  second  program  is  written  on  the  external  file  TAPE12  by 
routine  WRFILE.  This  data  contains  some  of  the  input  parameters 
and  computed  quantities  from  Program  I. 

It  is  the  user's  responsibility  to  provide  the  job  control 
language,  JCL,  necessary  to  manipulate  and  save  the  data  on 
TAPE12  to  transmit  it  to  Program  II.  Multiple  runs  with  Program  II 
may  then  be  made  with  the  same  dataset.  If  the  programs  are  to  be 
run  sequentially  in  the  same  job  stream,  no  additional  JCL  other 
than  to  see  that  the  data  is  not  lost  is  required.  Tape  rewinds 
are  taken  care  of  internally  in  the  program. 

There  are  nine  sections  of  data  as  written  from  routine 
WRFILE.  The  data  consists  of  one  section  for  each  of:  program 
indices,  program  constants,  u-velocity  panel  data,  thickness 
panel  data,  either  line  source  or  source  paneling  fuselage  data, 
pylon  data,  rack  data,  and  circular  and  elliptic  store  data. 

For  the  details  of  the  variables  written  see  the  program  listing 
in  Figure  A-l  of  Volumne  III.  All  variables  written  were  either 
residing  in  blank  or  labeled  common  at  some  point  in  the  program 
execution.  The  individual  descriptions  of  the  variables  may  be 
found  in  Appendix  B  of  Volume  III. 

When  using  both  noncircular  options  for  store  and  fuselage 
this  file  may  be  quite  voluminous.  For  the  sample  case  in 
Figures  11  using  137  source  panels  on  the  fuselage  and  144  source 
panels  on  the  elliptic  store  the  dataset  required  650  sectors 
(640  characters/sector)  of  external  storage  on  a  CYBER-175.  The 
largest  single  items  saved  are  the  aerodynamic  influence  coef¬ 
ficient  and  U,V,W,  component  velocity  matrices  for  the  store  to 
be  separated.  Each  of  these  four  contain  about  N(N+l)/2  elements. 
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When  two  elliptic  store  shapes  exist,  use  of  the  IUVW  option 
in  item  61  can  avoid  the  unnecessary  saving  of  these  as 
appropriate.  If  only  one  shape  is  to  be  separated,  the  index 
IUVW  can  be  set  to  IUVW=1  for  the  other  shape. 


3.5  Program  Error  Messages 

The  possible  occurrence  of  certain  fatal  errors  and 
irresolvable  conflicts  of  data  during  program  execution  has  been 
anticipated  in  Program  I.  If  such  errors  should  occur,  the 
program  has  been  designed  to  print  a  diagnostic  error  message 
and  halt  when  appropriate.  Each  message  with  the  routine  in 
which  the  corresponding  error  is  detected  is  presented  in  the 
following  table.  The  program  tests  which  resulted  in  the  section 
are  also  included. 


Program 

Routine  Action  Error  Message 

LDCALC  STOP  **  NSTRS  input  as  _ ,  program  limited 

to  seven  stores. 

(Test:  NSTRS  >  7) 


BDYGEN  STOP  **  At  base  of  body  radial  distance  to 

Mach  cone  emanating  from  body  nose  is 
less  than  maximum  body  radius  check 
input  data  and  Mach  number. 

(Test:  BETA*RADIUS<BODYL*XEND (NSEG) ) 


STOP  **  Mach  cone-body  meridian  intersection 

not  found  after  100  trials. 

(Test:  ITRY  >  100;  loop  from  40  to  39) 


BOOVEL  STOP  100 


(Test:  NER  >0,  indicates 
|ACNER,NER)  |<10’20  in  PAS001) 


90 


Routine 

Program 

Action 

Error  Message 

GEOM 

STOP  320 

Error  -  total  number  of  external  geometry 

axial  stations  exceeds  50  -  GEOM 

(test:  NX  >  50) 

STOP  330 

Error  - 

exceeds 

(Test: 

maximum  number  of  meridional  cuts 

33  -  GEOM. 

MAXKR  >  33) 

STOP  340 

Error  -  maximum  number  of  segments  exceeds 

5  -  GEOM. 

(Test:  NPUS  >  5) 

Warning 

only 

Warning 

GEOM  - 

(Test: 

-  Blank  common  required  for 

__  exceeds  dimension  •  _ . 

NATOT  >  NADIM) 

INVER1  STOP  **  Matrix  is  singular 

(Test:  A(I,I)«0.0) 

S OR PAN  STOP  230  **  Mach  number  less  than  1,  SORPAN 

(Test:  AMACH  <  1) 

STOP  220  Error  -  body  panel  slope  exceeds  Mach 

angle,  J  »  _  -  **  SORPAN  ** 

(Test:  1.-BETAL2*TAND2<0.) 

STORIO  STOP  707  Shape  data  not  input  for  all  stores  - 

**  STORIO  ** 

(Test:  NCOUNT  ^  NSTRS) 

WITHIN  STOP  **  Slope  of  wing  thickness  envelope  at 

leading  edge  not  positive  for  some 
chord vise  row. 

(Test:  THBTAL(JLE)  <  0.) 
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If  one  of  the  above  messages  is  printed  out  during  program 
execution,  the  user  should  consult  the  input  data  first  for 
possible  error.  The  program  variable  descriptions  in 
Appendices  A  and  B  of  Volume  III  may  also  be  helpful  for  the 
meanings  of  the  above  test  quantities. 

3.6  Program  I  Running  Times 

The  first  program  described  in  this  report  has  been  run  on 
both  the  CDC  6600  and  CYBER-175  machines.  Because  of  machine 
differences,  the  running  time  varies  from  one  machine  to  the 
other.  As  a  consequence,  only  an  approximate  running  time  can 
be  given.  The  first  program  can  be  run  only  once  for  multiple 
trajectory  calculations,  thus  its  percentage  of  the  total 
trajectory  calculation  time  will  vary. 

The  running  time  for  the  calculation  of  the  singularity 
strengths  associated  with  the  parent  aircraft  and  its  components 
is  a  function  of  a  number  of  factors,  some  of  which  are; 

(a)  Number  of  sources  and  doublets  required  to  represent 
circular  fuselage,  rack,  or  store  bodies. 

(b)  Number  of  source  panels  used  to  represent  the  non¬ 
circular  fuselage  or  elliptic  stores. 

(c)  Number  of  constant  u-velocity  panels  on  the  wing  and 
pylon. 

(d)  Number  of  traverses  used  to  define  the  nonlinear  shock 
shapes  associated  with  source  panel  modeling  of 
fuselage  or  stores. 

All  of  the  above  should  be  kept  to  a  minimum  required  to  resolve 
the  separating  store  loads.  For  the  sample  case  for  the  first 
program  in  Figure  14,  the  run  time  was  about  55  seconds  on  the 
CYBER-175. 


t: 
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4.  DETAILS  OF  THE  USE  OF  PROGRAM  II 


The  second  program  to  integrate  the  trajectory  of  the 
separated  store  consists  of  a  main  program  and  106  subroutines. 
Appendix  C  of  Volume  IV  presents  a  detailed  description  of  the 
flow  of  calculations  in  each  of  these  routines.  Table  C-l  in 
Appendix  C  .lists  these  subroutines  in  alphabetical  order  and 
gives  a  one-sentence  description  of  what  each  subroutine  does. 

A  listing  of  the  routines  in  Program  II  is  presented  in 
Figure  C-l.  A  general  flow  chart  of  the  main  program  for 
Program  II  (TRJTRY)  is  presented  in  Figure  15  of  this  report. 

The  program  as  written  in  Figure  C-l  adheres  to  ANSI  FORTRAN 
standards.  Only  the  first  card  is  specifically  for  the  CYBER-175 
series  machines.  The  following  sections  outline  the  flow  of 
calculations  as  presented  in  the  flow  chart  of  the  main  program, 
TRJTRY,  the  preparation  of  input  cards,  the  descriptions  of  the 
output,  and  any  special  conditions  and  messages  encountered  within 
the  program. 


4.1  General  Flow  Chart  of  Program  I I 

The  purpose  of  the  second  program  is  to  define  the  properties 
of  the  separating  store  including  any  empennages  and  to  integrate 
the  six-degree-of-freedom  equations  of  motion  which  are  derived 
in  Appendix  B  of  Reference  1.  The  separating  store  may  be  circular 
with  a  single  empennage  or  elliptic  in  shape  with  one  or  two  sets 
of  arbitrarily  oriented  fins.  The  elliptic  option  may  also  be 
used  for  a  circular  store. 

A  general  flow  chart  of  the  second  program  is  presented  in 
Figure  15.  Page  1  of  the  flow  chart  describes  the  sequence  of 
reading  the  external  dataset  containing  the  parent  aircraft  data 
and  of  reading  the  input  defining  the  separating  store  and  any 
empennages  on  it.  Constants  are  first  defined  and  additional 
heading  information  is  read  and  printed.  The  external  file 
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containing  the  parent  aircraft  and  store  body  information 
generated  in  Program  I  is  then  read  by  a  call  to  RDF1LE.  Any 
elliptic  store  data  or  noncircular  fuselage  data  are  copied  by 
that  routine  onto  TAPE10  and  TAPEll,  respectively.  The  air¬ 
craft  flight  conditions  are  input  and  compared  with  those  on  the 
dataset  from  Program  1;  if  a  mismatch  occurs  the  job  is  terminated. 

If  a  store  is  present,  the  program  reads  in  the  additional 
data  required  to  describe  the  store  to  be  separated.  If  an 
ejector  is  used,  the  force  histories  are  read.  The  store  mass 
and  inertia  properties  are  input  as  are  the  indices  specifying 
the  ejected  store  options.  If  the  ejected  store  is  circular  in 
shape  the  body  geometry  is  read.  If  a  noncircular  fuselage  or 
elliptic  stores  exist,  their  data  are  copied  sequentially  into 
blank  common  in  STRDAT. 

If  the  store  to  be  separated  is  circular  and  has  an  empennage, 
the  data  required  to  describe  the  empennage  are  input  and  sub¬ 
routine  SEMPIN  is  called  to  initialize  the  force  and  moment 
calculation.  If  the  store  to  be  separated  uses  the  elliptic 
model,  one  or  two  sets  of  arbitrarily  oriented  fins  can  be  input. 
The  locations  of  the  finned  sections  on  the  body  are  read  and  the 
appropriate  indices  of  body  panels  computed  in  IXBOD.  The  empen¬ 
nage  data  are  input  and  u-velocity  panel  influence  coefficient 
matrix  is  computed  in  CRFWBD.  If  two  empennages  exist,  all  data 
in  labeled  commons  necessary  to  restart  the  calculations  and  the 
influence  coefficients  are  written  on  TAPE 3  for  later  reuse. 

As  shown  on  the  Second  page  of  the  flow  chart,  Figure  15(b), 
the  program  next  reads  in  polynomials  specifying  the  thrust  time 
history  if  this  option  is  being  used. 

The  provision  is  made  in  the  program  for  prescribing  initial 
translational  and  rotational  velocities  relative  to  the  parent 


94 


aircraft.  These  are  next  read  and  the  initial  values  of  the  12 
dependent  variables  in  the  trajectory  calculation  computed.  These 
are: 

(1)  The  three  coordinates  of  the  store  center  of  moments 
relative  to  the  fuselage  nose. 

(2)  The  three  translational  velocity  components  of  the  store 
center  of  moments  relative  to  the  fuselage. 

(3)  The  three  store  rotational  velocities. 

(4)  The  three  angles  giving  the  store  orientation  relative 
to  the  fuselage. 

Reference  positions  of  the  store  nose,  center  of  moments,  and  base 
are  next  calculated  and  the  initial  and  final  trajectory  times,  as 
well  as  the  Integration  interval,  are  input.  If  the  initial  time 
is  not  zero,  then  the  trajectory  is  being  restarted  from  a  previous 
run  and  the  current  values  of  the  12  dependent  variables  obtained 
from  that  run  are  read  in.  The  last  steps  in  the  initialization  are 
to  compute  the  ejector  force  at  t_tINiTIAL  and  to  *nit*alize 
subroutine  ADAMS,  the  integration  routine. 

The  remainder  of  the  flow  chart.  Figure  15,  is  the  integration 
loop  of  the  program.  The  first  steps  are  to  calculate  the  aero- 
dynamic  forces  and  moments  acting  on  the  body  and  the  empennage(s) , 
if  present,  of  the  separating  store.  Different  methods  are  used 
for  each  of  these  calculations  depending  on  whether  the  circular 
or  elliptic  store  model  is  used.  If  the  store  is  circular,  the 
body  forces  and  moments  are  determined  by  using  the  equations 
presented  in  Appendix  A  of  Reference  1.  The  empennage  forces 
and  moments  are  determined  as  discussed  in  Section  5.3  and 
Appendix  I  of  Reference  6. 
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Tf  the  elliptic  store  model  is  used,  the  body  forces  and 
moments  are  computed  in  SF0RC2  using  the  body  source  panels 
described  in  References  12  and  13.  The  store  body  solution 
includes  the  presence  of  the  parent  aircraft  and  any  image  store 
resulting  from  shock  reflections  off  fuselage  or  wing  surfaces. 

If  a  single  empennage  is  present  the  forces  and  moments  are  com¬ 
puted  for  the  u-velocity  panels  in  DEM0N2  according  to  the  methods 
described  in  Reference  14.  This  solution  includes  the  presence 
of  the  parent  aircraft  and  any  body  nose  shock  reflections.  If  a 
second  empennage  is  present,  the  trailing  edge  vorticity  generated 
by  the  first  empennage  may  be  integrated  aft,  for  certain  con¬ 
figurations,  to  the  leading  edge  of  the  second  set  of  fins.  Due 
to  the  size  of  the  data  arrays  involved,  the  fuselage  data  must 
be  saved  temporarily  on  TAPE7  while  the  empennage  calculations 
are  performed.  The  data  for  each  set  of  empennages  must  be 
brought  in  from  external  storage  during  its  computation. 

One  of  the  options  in  the  computer  program  is  to  calculate 
a  wind-tunnel  captive-store  trajectory  as  opposed  to  a  free-flight 
trajectory.  It  is  customary  in  the  wind  tunnel  to  change  the 
store  orientation  relative  to  the  parent  aircraft  while  measuring 
the  forces  and  moments  in  order  to  approximately  account  for  the 
store  translational  motion  relative  to  the  aircraft.  The  computer 
program  also  does  this  during  the  force  and  moment  calculation. 
Thus,  if  a  captive-store  trajectory  is  being  calculated,  the 
next  step  in  the  program  is  to  put  the  store  back  to  its  correct 
orientation  and  call  subroutine  DIRCOS  in  order  to  calculate  the 
free-flight  direction  cosines  between  the  store  body  coordinate 
system  and  the  inertial  coordinate  system  fixed  to  the  fuselage. 

The  next  series  of  steps  determines  the  store  translational 
and  rotational  accelerations.  This  involves  solving  the  set  of 
six  simultaneous  equations  given  by  Equations  (B-16)  through 
(B-18)  and  (B-41)  through  (B-43)  of  Appendix  B  of  Reference  1 
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making  use  of  the  relationships  given  in  Section  J  of  that 
reference.  The  coefficient  matrix  is  first  calculated  and  then 
the  right-hand  sides  'hr e  determined.  Subroutine  INVER2  is  called 
to  solve  the  set  of  six  equations  for  the  accelerations.  The 
rates  of  change  of  the  orientation  angles  are  next  determined 
from  Equation  (B-l)  of  Appendix  B  of  Reference  1. 

A  check  is  next  made  to  see  if  output  is  to  be  printed  at 
the  end  of  an  integration  step.  If  output  is  not  required,  the 
integration  continues  by  calling  subroutine  ADAMS.  If  it  is 
required,  subroutine  SOUTPT  is  called.  Upon  return  from  this 
subroutine  a  check  is  made  to  see  if  the  time  is  equal  to  or 
greater  than  the  final  time  which  was  read  in  and,  if  so,  the 
trajectory  is  stopped.  If  it  is  not,  the  integration  is 
continued. 


4 . 2  Input  Data 

This  section  of  the  report  will  describe  in  detail  the 
preparation  of  the  input  data  deck  for  Program  II.  Only  the 
data  read  from  the  input  data  deck  is  described  here  and  not 
the  data  generated  in  Program  I  which  is  transferred  on  TAPE12. 
Program  II  reads  the  data  required  to  define  the  separating  store 
and  its  empennages  and  to  initialize  the  store  trajectory 
calculation. 

4.2.1  Input  formats 

The  format  for  the  input  data  for  Program  II  is  shown  in 
’’igure  16.  Four  lines  of  information  are  shown  for  each  item. 

The  first  line  gives  the  item  number,  how  many  and  when  the  cards 
are  read,  and  the  routine  in  the  program  in  which  the  data  is 
actually  read.  The  second  line  gives  the  program  variable  names, 
the  third  line  shows  the  card  column  fields  into  which  the  data 
are  to  be  punched,  and  the  fourth  line  shows  the  FORTRAN  format 
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type.  Data  punched  in  I  and  E  formats  are  right  justified  in 
the  fields  whereas  data  in  P  format  can  be  punched  anywhere  in 
the  field.  A  decimal  point  should  be  included  in  both  E-  and 
F-type  data. 

4.2.2  Input  descriptions  for  Program  II 

Two  input  files  are  read  by  Program  II.  The  first  is  read 
from  cards  via  unit  5  on  the  computer.  The  second  file  is  read 
from  an  unformatted  disk  file  on  unit  12  that  was  previously 
generated  by  Program  I.  The  descriptions  of  variables  here  are 
for  the  input  cards  read  on  unit  5. 

Item  number  1  has  an  index  NCARDS  which  indicates  how  many 
cards  of  information  are  to  follow  to  identify  the  run,  item 
number  2.  The  value  of  NCARDS  must  be  one  or  greater. 

Item  number  2  is  a  set  of  NCARDS  cards  containing  hollerith 
information  identifying  the  run.  This  information  may  start  and 
end  anywhere  on  the  card.  The  cards  are  reproduced  in  the  output 
just  as  they  are  read  in. 

Item  number  3  consists  of  one  card  containing  the  following 
flight  condition  information. 

Item  3: 

ALFA  fuselage  angle  of  attack,  degrees 

GAMF  fuselage  flight  path  angle,  degrees 

CMACH  Mach  number 

RHO  air  density  at  flight  altitude,  slugs  per  cubic  foot 

VINF  aircraft  free-stream  velocity,  feet  per  second 
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The  values  of  ALFA  and  CMACH  must  be  identical  to  those  values 
of  ALFAC  and  FMACH  read  in  item  3  of  Program  I.  The  values  of 
ALFA  and  CMACH  are  compared  with  ALFAC  and  FMACH  on  TAPE12.  If 
either  value  fails  to  compare,  program  execution  is  terminated. 
The  aircraft  is  assumed  to  be  flying  in  a  straight  line;  however, 
it  can  be  climbing  or  diving.  For  climbing  flight,  GAMF  is 
positive.  The  Mach  number  should  be  between  1.2  and  3.0  and  the 
angle  of  attack  should  not  exceed  10°  for  the  flow  models  to  be 
valid. 

Item  number  4  is  one  card  which  contains  eleven  indices. 

They  are: 

Item  4: 

NEJECT 
NSEG 

NSORCE 


NGAM 


NPOLY 


number  of  the  store  being  separated 

number  of  equal  length  segments  the  body  is  to  be 
broken  into  for  the  force  calculation;  NSEG  <  40; 
NSEGs0  for  elliptic  store  option 

number  of  line  sources  and  doublets  to  be  used  to 
model  the  store  volume  and  angle  of  attack  effects; 
NSORCE  <  100;  NSORCE=0  for  elliptic  store  option 

trajectory  to  simulate  wind-tunnel  captive-store 
trajectory? 

NGAM=0 ,  no 
NGAM=1 ,  yes 

number  of  polynomials  specifying  store  shape; 

1  <  NPOLY  <  7;  NPOLY=0  for  elliptic  store  option 


NROLL 


rolling  moment  to  be  calculated? 
NROLL=0,  no 
NROLL=l,  yes 


NEMP 


empennage  present? 

NEMP-0 ,  no 

NEMP=1 ,  yes,  single  set  of  fins 
NEMP=2,  yes,  for  elliptic  store  option  only, 
include  two  sets  of  fins 

NDAMP  damping  to  be  included  in  force  calculation? 

NDAMP=0,  no 
NDAMP=1 ,  yes 

NTHRUS  thrust  time  history  to  be  specified? 

NTRHUS=0 ,  no 
NTRHUS=1,  yes 

NTHETA  number  of  theta  stations  around  the  full  circumfer¬ 

ence  of  store  body  used  in  pressure  calculations; 
NTHETA  <  37.  NTHETA  must  further  be  an  odd  value. 
NTHETA=0  for  elliptic  store  option 

NJECTR  ejector  present? 

NJECTR=0,  no 
NJECTR=1,  yes 

The  index  NEJECT  must  correspond  to  one  of  the  store  numbers 
that  were  read  into  NUMSTR  in  item  54  of  Program  I;  that  is, 
NEJECT*NUMSTR(J) .  The  shape  number  of  this  store  is  NSHAPE (NEJSTR) 
where  NEJSTR=J.  If  NSHAPE (NEJSTR) £50  the  circular  store  option 
will  be  used.  The  NSHAPE  array  was  input  in  item  54  of  Program  I. 
The  variables  NSEG,  NSORCE,  NPOLY,  and  NTHETA  are  only  required 
by  the  circular  body  modeling  of  the  store. 

In  the  circular  store  body  force  and  moment  calculations, 
the  body  is  divided  into  NSEG  equal  length  segments.  Experience 
in  using  the  program  has  shown  that  20  body  segements,  NSEG*20, 
usually  yields  converged  forces  and  moments.  This  can  only  be 
checked  for  a  specific  case  by  varying  NSEG  and  comparing  results. 
To  minimize  calculation  time  NSEG  should  be  kept  as  small  as 
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possible.  For  elliptic  store  bodies,  NSEG  is  internally  set 
equal  to  the  number  of  rings  of  panels  as  determined  from  KFORX 
for  the  store,  input  item  72  of  Program  I. 

NSORCE  defines  the  same  information  as  MSOR  in  item  56  of 
the  first  program  for  the  store  to  be  separated.  Though  the 
variables  are  identical,  there  is  no  requirement  that  they  have 
the  same  value.  More  or  less  sources  and  doublets  may  be  used 
as  necessary  for  the  accuracy  of  the  predictions. 

The  index  NGAM  is  included  as  input  for  the  purpose  of  allowing 
the  program  to  be  used  to  compare  with  captive  store  data 
obtained  in  the  wind  tunnel.  Since  the  wind  tunnel  cannot  produce 
a  flow  where  the  store  sees  a  free- stream  velocity  coming  from  a 
different  direction  than  that  seen  by  the  parent  aircraft  the 
captive  store  case  must  be  handled  differently  by  the  computer 
program. 

The  number  of  polynomials  specifying  the  circular  store  shape, 
NPOLY,  is  the  number  required  to  specify  the  shape  from  the  store 
nose  to  its  base.  These  polynomials  are  used  in  the  force  and 
moment  calculation.  The  polynomials  are  of  the  form  given  by 
Equation  (1) .  The  value  of  NPOLY  may  be  equal  to  the  value  of 
NSPOLJ  read  under  item  57  in  Program  I.  It  is  used  only  by  the 
circular  store  option. 

The  next  index,  NROLL,  indicates  whether  or  not  the  rolling 
moment  for  a  store  with  an  empennage  is  to  be  calculated.  NEMP 
specifies  whether  there  are  empennages  present.  For  the  circular 
body  option,  only  one  empennage  may  be  modeled  using  the  slender- 
body  methods  detailed  in  Reference  6.  This  option  is  restricted 
to  monoplane  and  cruciform  fin  configurations.  For  the  elliptic 
store  option,  there  may  be  one  or  two  sets  of  empennages.  Each  of 
these  are  modeled  using  constant  u-velocity  panels  and  consist  of 
an  interference  shell  and  from  one  to  four  arbitrarily  oriented 


fins.  The  index,  NDAMP,  is  used  by  the  program  to  determine 
whether  or  not  aerodynamic  damping  in  pitch,  yaw,  and  roll  is 
to  be  included  in  the  force  and  moment  calculation.  The  index, 
NTHRUS,  specifies  whether  a  thrust  time  history  is  to  be  read. 

If  present,  items  29,  30,  and  31  are  read. 

The  index,  NTH ETA,  seta  the  number  of  stations  used  in  the 
pressure  integration  around  the  circumference  of  the  body  for 
the  circular  store  option.  This  includes  the  requirements  that 
the  first  and  last  points  be  the  same  and  that  the  number  be 
odd  to  maintain  body  symmetry.  The  last  index,  NJECTR,  determines 
whether  an  ejector  is  used  during  the  store  separation.  If  a 
blank  or  zero  is  input  the  option  is  not  used.  However,  setting 
NJECTR-1  causes  the  ejector  model  to  be  activated.  In  this  case, 
the  additional  items  5,  6,  7,  and  8  are  read  to  specify  the 
ejector  forces  and  moments  versus  time  or  distance.  The  value 
of  one  is  used  to  initiate  reading  of  the  input  values  in  routine 
EJECTR.  On  return  it  is  internally  reset  to  2  to  next  compute 
the  ejector  foot  location  and  forces  and  moments  at  time  equals 
zero.  It  is  again  internally  reset  to  NJECTR- 3  for  the  remainder 
of  the  trajectory  calculations. 

Item  numbers  5  through  8  are  used  to  specify  the  ejector 
time  or  displacement  histories.  If  NJECTR  in  item  4  is  input 
as  NJECTR- 0  these  items  are  omitted.  Item  number  5  is  used  to 
specify  the  ejector  options  as  follows: 

Item  5: 

NFEET  number  of  ejector  feet;  N7EET  <  2 

nstrke  switch  specifying  whether  time  or  stroke  will  be 

independent  variable: 

NSTRXE-0,  time  or  stroke  exceeding  limits  terminate 
ejector  forces  and  moments 
NSTRXE-1,  ejector  stroke  is  sole  criterion  for 
termination  of  ejector  forces 
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Exceeding  either  of  the  termination  criteria  will  cause  the 
ejector  forces  and  moments  to  be  set  to  aero.  If  time  is  to 
be  the  sole  criterion  for  ejector  termination,  use  the  option 
NSTRKE-0  and  set  STROKE (I)  in  item  6  equal  to  a  large  number. 

The  ejector  model  has  the  capability  of  modeling  1  or  2 
ejector  feet.  The  ejector  is  assumed  to  act  in  a  plane  perpen¬ 
dicular  to  the  fuselage  longitudinal  axis  (the  XB  axis,  see 
Figure  5) .  It  is  also  assumed  that  the  ejector  acts  in  the  plane 
containing  the  store  longitudinal  axis  when  the  store  is  in  the 
carriage  position.  The  ejector  foot  has  no  volume:  i.e.,  the 
ejector  foot  is  assumed  to  be  a  line.  The  rack/aircraft 
structure  is  rigid:  i.e.,  the  rack  transmits  100%  of  the  force 
to  the  store. 

The  next  two  items,  items  6  and  7,  are  repeated  as  a  set 
NFEET  times;  one  for  each  ejector  foot.  The  first  card,  item 
number  6,  specifies  the  ejector  location  and  movement  as  follows 
(1*1 ,NFEET) : 

Item  6: 

NEPOLY (I) 

XE(I) 


THETAE ( I ) 


STROKE (I) 


number  of  5th  order  polynomials  used  to  specify  the 
Ith  ejector  foot's  force  as  a  function  of  the 
independent  variable 

longitudinal  distance  along  the  store  centerline 
to  the  Ith  foot;  feet,  measured  positive  aft  from 
moment  center. 

angle  between  the  vertical  and  the  line  of  action 
of  the  ejector  foot;  measured  positive  counter¬ 
clockwise  viewed  from  rear,  degrees 

stroke  length  of  the  Ith  foot;  feet 
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The  polynomial  independent  variable  may  be  time  or  distance  as 
specified  by  NSTRKE.  The  sense  of  XE  is  positive  along  the  xg 
axis  in  Figure  17,  with  the  origin  shifted  to  the  store  moment 
center;  i.e.,  XE(I)  is  negative  if  the  foot  is  between  the  center 
of  gravity  and  the  store  nose.  The  positive  sense  of  XE  and 
THETAE  are  shown  in  the  sketch  below. 


The  second  card,  item  number  7,  contains  the  end  points  of 
the  NEPOLY ( I )  polynomials  as  follows: 

Item  7 : 

TEEND(I,J)  time  or  stroke  end  points  of  ejector  force  poly¬ 
nomials  for  Ith  ejector  and  Jth  polynomial, 

1  <  J  <  NEPOLY (I) 

The  dimensions  of  TEEND  are  of  time  in  seconds  or  feet  depending 
on  the  independent  variable  specified  by  the  option  NSTRKE .  Both 
variables  are  assumed  to  start  at  zero  at  time  equals  zero.  Only 
NEPOLY (I)  values  need  to  be  input  on  each  card. 

Item  number  8  consists  of  NEPOLY (I)  cards  for  each  ejector 
foot,  each  containing  the  coefficients  of  a  single  polynomial  as 
follows: 
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Item  8: 


AKE (I ,  J)  each  card  contains  coefficients  Cq,C^,C2/C3,C4,Cj 
for  a  single  5th  order  polynomial 

The  units  are  force  F(I)  in  pounds  with  the  independent  variable 
X  of  time  or  distance  in  seconds  or  feet.  The  coefficients  are 
for  a  polynomial  of  the  form: 

F (I)  =  CQ  +  CjX  +  C2X2  +  C3X3  +  C4X4  +  C5X5 


where  X  is  the  independent  variable  of  time  or  stroke  length. 

NFEET 

This  group  contains  E  NEPOLY(I)  cards.  The  cards  begin  with 

1=1 

the  coefficients  for  the  polynomials  consistent  with  the  input 
for  1=1.  Coefficients  for  the  second  foot  (1=2)  if  NFEET=2 
follow  immediately  after  those  for  the  first  foot. 


Item  number  9  is  one  card  and  specifies  the  store  mass  and 
inertia  characteristics.  The  quantities  are: 


Item  9: 

SMASS 

FIXX 

FIYY 

FIZZ 

FIYZ 

FIXZ 

FIXY 


store  mass,  slugs 

2 

I  moment  of  inertia,  slug-ft 

Aa 

2 

Iyy  moment  of  inertia,  slug-ft 

2 

I  moment  of  inertia,  slug-ft 
z  z 

2 

I  product  of  inertia,  slug-ft 

2 

I  product  of  inertia,  slug-ft 
xz 

2 

I  product  of  inertia,  slug-ft 
xy 


The  equations  defining  the  moments  and  products  of  inertia  are 
given  by  Equation  (B-36)  of  Appendix  B  of  Reference  1. 
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The  one  card  of  item  number  10  contains: 

Item  10: 

"  t 

XMOM  location  along  store  axis  about  which  the  pitching 

and  yawing  moments  are  to  be  taken,  negative  behind 
nose,  feet;  same  point  about  which  moments  of 
inertia  are  taken;  for  elliptic  store  XMOM  must  be 
the  same  as  REFX  input  in  item  73  of  Program  I 

XBAR  x  location  of  store  c.g.  measured  from  moment 

center,  feet;  positive  forward  toward  the  nose 

YBAR  y  location  of  store  c.g.  measured  from  store  axis, 

feet;  positive  to  the  right  looking  forward 

ZBAR  z  location  of  store  c.g.  measured  from  store  axis, 

feet;  positive  down  viewed  from  the  rear 

The  next  two  items  of  input,  items  11  and  12,  describe  the 

shape  of  a  separated  circular  store.  They  are  omitted  if 

NSHAPE (NEJSTR) >50.  The  quantities  contained  on  the  cards  are: 

Item  11; 

XEND(J)  x/£  of  end  points  of  polynomials  specifying  shape 

of  ejected  store,  NPOLY  values 

Item  12: 

COEF ( j , K)  coefficients  of  polynomials  specifying  shape 

These  data  are  a  repeat  of  the  data  input  as  items  58  and  59  in 

Program  I. 

Item  number  13  contains  one  quantity  which  is 

CA  store  axial-force  coefficient;  reference  area  is 

store  maximum  cross-sectional  area  for  the  circular 
store  option  and  equal  to  REFAR  of  item  73  of 
Program  I  for  elliptic  store  option 
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The  store  axial-force  coefficient  is  not  calculated  by  the 
computer  program  so  it  is  required  input. 


The  next  two  items  of  input  data,  items  14  and  15,  are 
included  in  the  input  data  deck  for  the  circular  store  option 
if  an  empennage  is  present;  NEMP=1,  in  item  4  and  NSHAPE (NEJSTR) 
less  than  or  equal  to  50. 

Item  14 ; 

IPLNR  IPLNR=0,  cruciform  empennage 

IPLNR=1 ,  planar  empennage 

MSF  number  of  spanwise  control  points  on  each  fin; 

must  be  odd  and  5  <  MSF  <  11 

For  most  stores,  MSF=5  has  been  found  to  give  accurate  results. 
The  larger  the  fin  span  to  body  radius  ratio,  the  more  points 
required. 

Item  15 ; 

XTAIL  x  location  at  which  empennage  forces  are  to  act 

measured  from  store  nose,  feet;  negative  nuftber 

RADAV  average  store  body  radius  in  empennage  regidn, 

feet;  positive  number 

FINSS  tail  fin  semispan,  measured  from  body  longitudinal 

axis,  feet;  positive  number 

PHIROL  initial  fin  orientation,  degrees;  0°  <  PHIROL  <  90°; 

PHIROL  =  0°  if  fins  vertical  and  horizontal  for  a 
cruciform  empennage  or  horizontal  for  a  planar 
empennage 

CLALPH  lift-curve  slope  of  two  exposed  panels  joined 

together,  per  radian;  reference  area  is  store 
maximum  cross-sectional  area 
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The  location  at  which  the  empennage  forces  are  assumed  to  act, 

XTAIL,  is  arbitrary.  It  can  be  estimated  using  Chart  10  of  Refer¬ 
ence  15.  The  lift-curve  slope,  CLALPH,  can  be  estimated  using  the 
method  described  as  follows.  The  value  of  the  tail-fin  lift-curve 
slope,  CL^,  can  determined  using  Reference  15;  however,  the  value 
of  CLa  from  Chart  8  of  that  reference  must  be  modified  in  order  to 
use  the  empennage  force  calculation  method  incorporated  into  the 
computer  program.  This  method  assumes  full  lift  carryover  onto  the 
body  due  to  the  tail  fins.  At  supersonic  speeds  this  is  not  the 
case,  particularly  when  the  body  ends  at  or  near  the  trailing  edge 
of  the  fins. 

A  method  for  determining  the  value  of  the  lift-curve  slope 
which  accounts  for  the  reduced  lift  carryover  will  now  be  presented. 
The  modified  lift-curve  slope  is  calculated  using  Charts  1,  4,  and 
8  of  Reference  15  and  the  following  equation: 

^(B).  +  KB(T). 

- ±_ - i  (10) 

KT(B)  +  KB(T) 


The  value  of  Cl  is  obtained  from  Chart  8.  The  aspect  ratio 
used  is  that  of  two  of  the  exposed  panels  joined  together.  The 
reference  area  (ST)  used  in  forming  this  coefficient  is  the  area 
of  the  two  fins  joined  together.  Thus,  it  is  multiplied  by  the 
ratio  of  this  area  divided  by  the  reference  area  (SR)  used  in 
the  computer  program,  the  store  maximum  cross-sectional  area. 
This  is  the  ratio  ST/SR. 

The  quantities  in  the  bracketed  tern  are  lift  ratios  and  are 
obtained  from  Charts  1  and  4  of  Reference  15.  The  number  sub¬ 
scripts  in  Equation  (10)  refer  to  the  chart  number.  The 
denominator  in  this  term  is  the  sum  of  the  lift  ratios  for  full 
lift  carryover  onto  the  body.  The  numerator  accounts  for  only 
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partial  lift  carryover  onto  the  body.  In  determining  the  value 
of  Kb  (T)  4  from  Chart  4,  (CLa)w  is  C^St/Sr,  the  quantity 
multiplying  the  bracketed  term  in  Equation  (10) . 

It  is  recommended  that  Chart  4  be  used  for  low-aspect-ratio 
as  well  as  high-aspect-ratio  tail  fins.  Reference  15  recommends 
that  Chart  4  only  be  used  when  the  aspect  ratio  parameter  defined 
there  is  greater  than  4.0.  The  use  of  this  chart  in  determining 
KB(T)4  for  vaiues  of  the  aspect  ratio  parameter  of  3.0  and  2.25 
has  resulted  in  good  agreement  with  experimental  data  for  two 
different  stores  in  uniform  flow.  Items  14  and  15  are  omitted 
for  elliptic  store  shapes,  NSHAPE (NEJSTR)  greater  than  50. 

The  next  13  items  of  input  are  used  to  define  the  empennage (s) 
for  the  elliptic  store  option.  Omit  items  16  through  28  if  NEMP 
equals  zero  or  NSHAPE (NEJSTR)£50.  Item  number  16  is  used  to 
partition  the  store  into  body  alone  and  body-fin  sections  and 
contains  the  following  information  on  one  card. 

Item  16: 

NBOD  number  of  body  and  body-fin  sections;  NBOD  <  5 

XBOD(I)  trailing  edge  of  Ith  body  or  body-fin  section, 

feet;  1=1, NBOD 

LFIN(I)  logical  variable  defining  Ith  section  to  be  body- 

alone  or  body-fin  section;  1=1, NBOD 
LFIN(I) =T(true) ,  section  contains  both  body  and 

fins 

LFIN (I) =F (false) ,  section  contains  only  body 

The  data  on  this  card  are  used  to  partition  the  body  into  inter¬ 
fering  and  noninterfering  sections.  The  values  of  XBOD  define 
the  body  sections.  XBOD(l)  is  the  end  of  the  first  body  section. 
The  XBOD  values  are  compared  with  the  XJ  values  of  the  source 
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panel  boundaries  to  determine  those  panel  rings  which  interfere 
with  the  fins.  They  must  correspond  to  panel  edges. 

The  next  12  items  define  an  empennage  consisting  of  from 
one  to  four  fins  and  a  body  interference  shell.  Include  one  set 
of  items  17  through  28,  as  required,  for  each  of  NEMP  empennages. 
They  define  the  layout  of  the  constant  u-velocity  panels  and  the 
geometric  parameters  not  available  in  the  store  description. 

Item  number  17  consists  of  one  card  with  the  alphanumeric 
description  of  the  empennage. 

Item  number  18  consists  of  one  card  with  the  following 
indices.  The  fin  number  system  is  shown  in  Figure  18. 


Item  18 ; 
MSWR 

MSWL 

MSWU 

MSWD 

NCW 


number  of  constant  u-velocity  panels  in  spanwise 
direction  on  right  (R)  or  first  (1)  fin; 

1  <  MSWR  <  19 

number  of  constant  u-velocity  panels  in  spanwise 
direction  on  left  (L)  or  second  (2)  fin; 

0  <  MSWL  <  19 

number  of  constant  u-velocity  panels  in  spanwise 
direction  on  upper  (U)  or  third  (3)  fin; 

0  <.  MSWU  <  19 

number  of  constant  u-velocity  panels  in  spanwise 
direction  on  lower  (D)  or  fourth  (4)  fin; 

0  <  MSWD  <  19 

number  of  constant  u-velocity  panels  in  chordwise 
direction,  same  for  all  four  fins 
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number  of  constant  u-velocity  panels  in  the 
longitudinal  direction  on  the  surface  of  the  body 
over  the  body  interference  length ,  BIL.  BIL  is 
input  in  item  23. 

number  of  constant  u-velocity  panels  on  the  full 
circumference  of  the  body.  If  NBDCR=0,  NBDCR  is 
set  to  KRADX-1  from  the  store  body  description. 

spanwise  breaks  in  fin  option; 

LVSWP=0,  no  breaks  in  fin  leading  or  trailing  edge 
sweep  angles;  equal  spanwise  spacings  of 
panels 

LVSWP=1,  up  to  19  breaks  in  fin  leading  or  trailing 
edge  sweep  angles;  that  is,  up  to  19 
unequal  spanwise  spacings 

The  spanwise  number  of  panels  per  lifting  surface  is  used  to 
determine  the  presence  of  each  of  the  four  fins.  Only  the  first 
fin,  corresponding  to  MSWR,  must  have  a  nonzero  value.  A  zero 
value  for  the  remaining  spanwise  panel  numbers  eliminates  the 
lifting  surface.  The  spanwise  number  of  panels  may  differ  from 
one  fin  to  another  but  the  chordwise  number  NCW  is  the  same  for 
all.  The  specification  of  a  nonzero  number  of  body  interference 
panels,  NBDCR,  on  the  circumference  may  be  used  to  override  the 
internal  default  value  which  is  set  equal  to  the  number  of  source 
panels  around  the  circumference.  The  default  value  is  recommended. 
The  number  of  body  interference  panels  in  the  axial  direction, 

NCWB,  should  be  set  equal  to  the  number  of  chordwise  panels  on 
the  fins  except  where  the  interference  shell  length,  BIL,  spanned 
by  the  lifting  surfaces  differs  from  the  fin  root  chord.  In 
this  case,  NCWB  follows  the  rules  for  defining  NCWB  for  the  wing 
in  item  9  of  Program  I.  The  control  index  LVSWP  indicates 
whether  breaks  in  leading-edge  and/or  trailing-edge  sweeps  are  to 
be  input.  If  LVSWP=0  the  sweep  angles,  SWLEP,  SWTEP,  SWLEV,  and 
SWTEV  are  input  in  items  19  and  20  and  items  25  through  28  are 


NCWB 


NBDCR 


LVSWP 
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not  included  in  the  input  deck.  If  LVSWP=1,  the  angles  SWLEP , 
SWTEP,  SWLEV ,  and  SWTEV  in  items  19  and  20  need  not  be  specified, 
but  items  25  through  28  must  be  included  in  the  input  deck. 

The  numbers  of  spanwise  and  chordwise  panels  are  limited 
by  the  dimensions  of  the  program.  Three  dimensions  within  the 
program  limit  the  number  of  panels.  The  total  number  of  u-velocity 
panels  on  the  fins  and  interference  shell  is  limited  to  250  with 
a  required  blank  common  array  size  of  62,500  decimal.  The  total 
number  of  u-velocity  panels  on  the  interference  shell  is  limited 
to  100.  The  number  of  spanwise  panels  on  the  fins  is  limited 
to  19.  Thus,  the  program  limitations  are: 

(MSWR+MSWL+MSWU+MSWD) *NCW  +  NBDCR*NCWB  <  250 


with 


NBDCR*NCWB  <  100 


and 


MSWR 

MSWL 

MSWU 

MSWD 


>  <  19 


The  next  two  items  of  input,  items  19  and  20,  specify  the 
root  chords,  spans,  and  sweep  angles  for  the  two  pairs  of  fins, 
fins  1  and  2  and  fins  3  and  4,  shown  in  Figure  18.  If  there  are 
breaks  in  leading-  and  trailing-edge  sweep,  the  sweep  angles  in 
items  19  and  20  are  input  as  zero  and  items  25  through  28  as 
required  are  included  in  the  input  deck. 

For  the  constant  sweep  option,  item  number  19  defines  the 
planforms  of  both  the  right  and  left  (horizontal)  lifting  surfaces 
(fins  1  and  2).  It  consists  of  one  card  with  the  following 
information. 
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Item  19: 
CRP 

SWLEP 

SWTEP 


B2 


horizontal  fin  root  chord  at  the  fin-body  juncture, 
feet 

horizontal  fin  leading-edge  sweep  angle  measured  in 
fin  planform  plane,  positive  for  sweep  back,  degrees 

horizontal  fin  trailing-edge  sweep  angle  measured 
in  fin  planform  plane,  positive  for  sweep  back, 
degrees 

exposed  horizontal  fin  semispan,  feet 


For  the  constant  sweep  option,  item  number  20  defines  the 
planforms  of  both  the  upper  and  lower  (vertical)  lifting  surfaces 
(fins  3  and  4) .  It  is  included  in  the  input  deck  if  either  MSWU 
or  MSWL  of  item  18  are  greater  than  zero  and  consists  of  one 
card  with  the  following  information. 


Item  20: 
CRPV 

SWLEV 

SWTEV 

B2V 


vertical  fin  root  chord  at  the  fin-body  juncture, 
feet 

vertical  fin  leading-edge  sweep  angle  measured  in 
fin  planform  plane,  positive  for  sweep  back,  degrees 

vertical  fin  trailing-edge  sweep  angle  measured 
in  fin  planform  plane,  positive  for  sweep  back, 
degrees 

exposed  vertical  fin  semispan,  feet 


The  designation  of  horizontal  and  vertical  fins  used  here 
only  has  meaning  for  planar  and  cruciform  fins.  See  Figure  18 
for  the  fin  number  definitions  for  interdigitated  and  arbitrarily 
located  fins.  It  is  a  program  limitation  that  requires  that  fin 
geometry  must  be  input  in  pairs  with  this  description.  That  is, 
fins  1  and  2  have  the  same  geometry  as  do  fins  3  and  4.  The 
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root  chord  for  the  fin-body  combination  is  the  line  formed  by 
the  junction  of  the  lifting  surfaces  and  the  body.  The  semispan 
is  measured  from  the  root  chord.  If  different  fin  geometry  is 
required,  use  the  LVSWP=1  option. 

Item  number  21  consists  of  one  card  with  the  following  fin 
angular  information.  This  card  specifies  an  interdigitated  tail 
configuration  as  shown  in  Figure  18.  If  there  are  less  than 
four  fins  or  if  the  configuration  is  cruciform  with  the  fins 
vertical  and  horizontal,  PHIDIH  and  THETIT  are  input  as 
PHIDIH=THETIT=0. 0°  and  item  22  is  included  in  the  input  deck. 

The  fin  deflection  angles,  however,  are  included  in  item  21. 

Item  21: 

PHIDIH  dihedral  angle  associated  with  interdigitated  fin, 

measured  from  the  horizontal  axis  (see  Figure  18), 

0°  <  PHIDIH  <  90° 

THETIT  polar  angle  of  meridian  associated  with  intersection 

of  interdigitated  fin  with  store  body,  measured  from 
horizontal  axis  (see  Figure  18),  0°  <  THETIT  <.  90° 

DELR  6r,  deflection  angle  of  horizontal  right  fin  (fin  1) ; 

positive  trailing  edge  down,  degrees 

DELL  deflection  angle  of  horizontal  left  fin  (fin  2) ; 

positive  trailing  edge  down,  degrees 

DELU  6^,  deflection  angle  of  vertical  right  fin  (fin  3) ; 

positive  trailing  edge  down,  degrees 

DELD  6d,  deflection  angle  of  vertical  lower  fin  (fin  4) ; 

positive  trailing  edge  down,  degrees 
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If  case  involves  interdigitated  fins: 

6r  applies  to  right  upper  fin,  trailing  edge  down  is  positive 
$L  applies  to  left  lower  fin,  trailing  edge  down  is  positive 
applies  to  right  lower  fin,  trailing  edge  down  is  positive 
6d  applies  to  left  upper  fin,  trailing  edge  down  is  positive 


The  representation  of  the  angles  PHIDIH  and  THETIT  for  the 
interdigitated  tails  on  an  elliptic  cross  section  body  is  given 
in  Figure  18.  The  meaning  of  the  interdigitated  angles  for  each 
of  the  four  fins  is  defined  on  the  figure.  The  deflection 
angles  are  defined  only  for  a  rigid  rotation  of  the  entire  con¬ 
trol  surface.  These  angles  are  used  only  in  the  specification 
of  the  surface  boundary  conditions. 

If  an  interdigitated  tail  configuration  is  not  used 
(THETIT=0 . 0) ,  item  number  22  defines  the  dihedral  and  polar 
angles  which  locate  each  of  the  four  fins. 


Item  22; 

PHIFR 

PHIFL 

PHIFU 

PHIFD 

THETR 

THETL 


dihedral  angle  of  right  fin  (fin  1) ,  degrees 

dihedral  angle  of  left  fin  (fin  2) ,  degrees 

dihedral  angle  of  upper  fin  (fin  3) ,  degrees 

dihedral  angle  of  lower  fin  (fin  4) ,  degrees 

polar  angle  of  intersection  of  right  fin  (fin  1) 
with  body,  degrees 

polar  angle  of  intersection  of  left  fin  (fin  2) 
with  body,  degrees 

polar  angle  of  intersection  of  upper  fin  (fin  3) 
with  body,  degrees 


THETU 


THETD  polar  angle  of  intersection  of  lower  fin  (fin  4) 

with  body,  degrees 

All  angles  are  measured  from  the  horizontal  axis  (Y-axis) .  The 
right  and  upper  fin  angles  are  measured  counterclockwise  from 
the  positive  Y-axis.  The  left  and  lower  fin  angles  are  measured 
counterclockwise  from  the  negative  Y-axis.  This  designation 
is  based  on  laying  out  cruciform  fins  in  symmetric  pairs  where 
PHIFR«PHIFL=THETR=THETL=0°  and  PHIFU=PHIFD=THETU=THETD=90° . 

Item  numbers  23  and  24  define  some  additional  empennage 
parameters  as  follows. 

Item  23: 

BIL  length  of  body  influenced  by  fins  to  account  for 

interference.  For  fins  with  no  trailing  edge  sweep, 
BIL=CRP 

VRTMAX  maximum  magnitude  of  vortex  induced  velocities 

included  in  flow  tangency  condition;  use  0.35  as 
default  value 

Item  24 : 

NVRTIN  option  to  override  vortex  tracking  calculation  for 

multiple  empennages; 

NVRTIN=0,  vortex  from  forward  lifting  surfaces  are 
integrated  back  to  leading  edge  of  aft 
fins 

NVRTIN=1,  no  vortex  influence  is  computed  for  aft 
fins 

The  variable  BIL  follows  the  same  definition  rules  as  BODYPL  in 
item  10  of  Program  I.  See  item  2,  Appendix  J  of  Reference  14  for 
implications  of  use  of  variable  VRTMAX.  NVRTIN  has  no  meaning 
when  only  one  empennage  exists  or  for  integration  of  vortex  paths 
aft  of  second  empennage. 
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Item  numbers  25  through  28  define  the  optional  planform 
descriptions  for  fins  with  breaks  in  leading-edge  and/or  trailing- 
edge  sweep.  They  are  only  included  in  the  input  deck  when 
LVSWP-1  in  item  18.  The  following  quantities  are  defined  for 
each  of  the  four  fins. 

Item  number  25  is  the  optional  input  for  the  right  fin 
(fin  1).  It  is  included  if  LVSWP=1  in  item  18. 

Item  25: 

YRT(KJ)  distance  from  fin  root  chord  to  the  constant 

u-velocity  panel  outboard  side  edge  on  the  right 
fin;  1  <  KJ  <  MSWR+1;  (MSWR  <  19);  YRT(l)  =  0.0; 

YRT (MSWR+1)  =  B2  of  item  19 

VSWLER(KJ)  leading-edge  sweep  of  fin  between  YRT(KJ-l)  and 
YRT (KJ) ,  positive  for  sweep  back,  degrees; 

1  <  KJ  <  MSWR+1;  VSWLER(l)  =0.0 

VSWTER (KJ)  trailing-edge  sweep  of  fin  between  YRT(KJ-l)  and 
YRT(KJ),  positive  for  sweep  back,  degrees; 

1  <  KJ  <  MSWR+1;  VSWTER (1)  =0.0 

The  optional  input  of  this  item  is  associated  with  a  fin-body 
combination  with  breaks  in  leading-edge  and/or  trailing-edge 
sweeps.  This  input  should  be  used  also  for  configurations  if 
the  constant  u-velocity  panel  side  edges  are  to  be  laid  out 
with  user-determined  unequal  spacings. 

Item  number  26  is  the  optional  input  for  the  left  fin 
(fin  2).  It  is  included  if  LVSWP=1  and  MSWL>0  in  item  18. 
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Item  26 


YLT(KJ)  distance  from  fin  root  chord  to  the  constant 

u-velocity  panel  outboard  side  edge  on  the  left 
fin,  1  <  KJ  <  MSWL;  (MSWL  <  19);  YLT(l)  »  0.0; 

YLT (MSWL+1)  *  -B2  of  item  19 

VSWLEL(KJ)  leading-edge  sweep  of  fin  between  YLT(KJ-l)  and 
YLT (KJ) ,  negative  for  sweep  back,  degrees; 

1  <  KJ  <  MSWL+1;  VSWLEL(l)  =  0.0 

VSWTEL(KJ)  trailing-edge  sweep  of  fin  between  YLT(KJ-I)  and 
YLT { KJ ) ,  negative  for  sweep  back,  degrees; 

1  <  KJ  <  MSWL+1;  VSWTEL (1) 

This  optional  input  accompanies  item  25  and  is  associated  with 
the  left  fin  or  fin  2.  The  Y  spanwise  values  vary  from  0.0  at 
the  root  chord  in  the  negative  direction  to  the  last  value  for 
YLT  equal  to  the  negative  semispan,  -B2. 

Item  number  27  is  the  optional  input  for  the  upper  fin 
(fin  3).  It  is  included  if  LVSWP=1  and  MSWU>0  in  item  18. 

Item  27 ; 

ZUT(KJ)  distance  from  fin  root  chord  to  the  constant 

u-velocity  panel  outboard  side  edge  on  the  upper 
fin,  1  <  KJ  <  MSWU+1 ;  (MSWU  <  19);  ZUT(l)  «  0.0; 
ZUT (MSWU+1)  «  B2V  of  item  20 

VSWLEU(KJ)  leading-edge  sweep  of  fin  between  ZUT(KJ-l)  and 
ZUT (KJ) ,.  positive  for  sweep  back,  degrees; 

1  <  KJ  <  MSWU+1;  VSWLEU (1)  =0.0 

VSWTEU (KJ)  trailing-edge  sweep  of  fin  between  ZUT(KJ-l)  and 
ZUT (KJ) ,  positive  for  sweep  back,  degrees; 

1  <  KJ  <  MSWU+1;  VSWTEU (1)  =0.0 
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This  information  is  similar  to  items  25  and  26.  This  input 
should  be  used  for  breaks  in  sweep  or  unequal  spacings.  Though 
the  variable  ZUT  is  used  to  represent  the  distance  from  the  fin 
root  chord  to  the  outboard  constant  u-velocity  panel  edges  on 
the  upper  fin  (fin  3) ,  this  variable  is  used  in  exactly  the  same 
sense  as  YRT  to  define  the  spanwise  direction  in  the  plane  of 
the  fin. 

Item  number  28  is  the  optional  input  for  the  lower  fin 
(fin  4).  It  is  included  if  LVSWP=1  and  MSWD>0  in  item  18. 

Item  28: 

ZDT(KJ)  distance  from  fin  root  chord  to  the  constant 

u-velocity  panel  outboard  side  edge  on  the  lower 
fin,  1  <  KJ  <  MSWD+1 ;  (MSWD  <  19);  ZDT(l)  -  0.0; 

ZDT (MSWD+1)  =  -B2V  of  item  20 

VSWLED(KJ)  leading-edge  sweep  of  fin  between  ZDT(KJ-l)  and 
ZDT(KJ),  negative  for  sweep  back,  degrees; 

1  <  KJ  <  MSWD+1;  VSWLED(l)  *  0.0 

VSWTED(KJ)  trailing-edge  sweep  of  fin  between  ZDT(KJ-l)  and 
ZDT (KJ) ,  negative  for  sweep  back,  degrees; 

1  <  KJ  <  MSWD+1;  VSWTED (1)  =0.0 

This  optional  information  is  similar  to  items  25  through  26. 

The  variable  ZDT  again  only  represents  the  distance  from  the  fin 
root  chord  to  the  outboard  constant  u-velocity  panel  edges  in 
the  plane  of  the  lower  fin  (fin  4) .  The  Z  spanwise  values  vary 
in  the  negative  direction  from  0.0  at  the  -  oot  chord  to  the 
last  value  for  ZDT  equal  to  the  negative  r  aispan,  -B2V. 

The  next  three  items  of  input  data,  item  numbers  29,  30, 
and  31,  are  included  in  the  input  data  deck  if  a  thrust  time 
history  is  to  be  specified,  NTHRUS=1,  in  item  4.  The  input 
quantities  are: 
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Item  29: 


NTPOLY  number  of  polynomials  specifying  thrust  time 

history;  1  <  NTPOLY  <  5 

Item  30: 

TEND(J)  final  times  over  which  the  NTPOLY  polynomials 

apply;  J=l, NTPOLY 

Item  31: 

TC(J,K)  coefficients  of  the  polynomials;  K=l,6;  J=l, NTPOLY 

The  thrust  force  as  programmed  acts  along  the  store 
longitudina-  axis  and  is  specified  by  a  series  of  polynomials 
of  the  form 


FT  =  l  V 
1  k=l  * 


where  FT  is  the  thrust  force  in  pounds.  Up  to  five  polynomials 
can  be  used  to  specify  the  time  history  and  item  29  contains  the 
number  of  polynomials,  NTPOLY.  Item  30  is  one  card  and  contains 
the  NTPOLY  values  of  the  final  time  t  over  which  each  of  the 
polynomials  applies.  That  is,  the  first  polynomial  is  used  from 
t  *  0  to  t  *  TEND(l),  the  second  from  t  ■  TEND (1 )  to  t  ■  TEND (2 ) , 
and  so  on.  Item  31  is  a  set  of  NTPOLY  cards  which  specify  the 
six  a^  coefficients,  TC(J,K),  for  the  polynomials. 

Item  number  32  is  one  card  which  contains  the  initial  veloc¬ 
ities  of  the  store  relative  to  the  parent  aircraft.  The  six 
quantities  are: 
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Item  32: 


I 

i 


I 


VXZERO 

initial  forward  velocity  of  the  store  parallel 
store  longitudinal  axis,  ft/sec 

to 

VYZERO 

initial  lateral  velocity  of  the  store  perpendicular 
to  store  longitudinal  axis,  ft/sec 

VZZERO 

initial  downward  velocity  of  the  store  perpen¬ 
dicular  to  store  longitudinal  axis,  ft/sec 

VAR (4) 

store  initial  rolling  rate,  p,  radians/ sec 

VAR ( 5 ) 

store  initial  pitching  rate,  q,  radians/sec 

VAR (6) 

store  initial  yawing  rate,  r,  radians/sec 

At  the  beginning  of  the  trajectory,  the  store  is  oriented  such 
that  its  y-axis  shown  in  Figure  17  is  parallel  to  the  yB 
fuselage  axis  shown  in  Figure  5.  Thus,  the  x,z  store  plane  is 
parallel  to  the  xB,zfi  fuselage  plane  and  the  initial  vertical 
translational  velocity  of  the  store,  VZZERO,  is  in  the  x,z  store 
plane  perpendicular  to  the  store  x-axis.  The  velocities 

VXZERO,  VYZERO,  and  VZZERO  are  positive  in  the  positive  x,  y,  and 

z  directions. 

The 

next  card,  item  number  33,  contains  three  items. 

They 

are: 

Item  33: 

DTIME 

integration  interval,  seconds 

TIMEI 

initial  time,  seconds 

TIMEF 

final  time,  seconds 

The  first,  DTIME,  is  the  integration  interval  to  be  used  in  the 
integration  subroutine,  ADAMS.  The  version  of  this  routine  which 
is  included  in  the  present  program  does  not  adjust  this  interval 
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to  satisfy  certain  accuracy  requirements  but  uses  the  input  value 
of  DTIME.  A  value  of  0.025  or  0.01  seconds  should  probably  work. 
However,  different  values  should  be  tried  and  the  trajectory  results 
compared  to  determine  that  a  rmall  enough  value  is  being  used. 

Finned  stores  probably  require  a  smaller  value  than  finless  stores 
because  of  the  large  moments  produced  by  the  empennage.  The 
initial  time,  TIMEI,  must  be  input  as  0.0  unless  a  trajectory  is 
to  be  restarted  using  the  last  page  of  output  from  a  previous  run 
to  obtain  the  initial  conditions.  Then,  TIMEI  should  be  given  the 
value  that  appears  on  that  page  of  output.  The  final  time,  TIMEF , 
is  the  time  at  which  the  trajectory  calculation  is  to  be  terminated. 
Except  for  very  unusual  trajectories,  a  value  of  0.5  to  0.7  second 
is  normally  sufficient  to  determine  if  a  store  will  or  will  not 
clear  the  aircraft. 

If  TIMEI  and  TIMEF  are  both  input  as  zero,  no  trajectory  calcu¬ 
lation  will  be  performed.  However,  the  store  load  distributions 
and  forces  and  moments  will  be  calculated  with  the  store  in  its 
initial  position.  This  feature  can  be  useful  in  checking  out  the 
entire  input  data  deck  prior  to  running  a  trajectory  or  for  study¬ 
ing  store  loads  at  specific  points. 

Item  number  34,  the  last  item  of  input,  is  input  only  if  a 
trajectory  is  being  restarted,  TIMEI  ^  0.  This  item  consists  of 
two  cards  with  VAR(l)  through  VAR (8)  on  the  first  card  and  VAR (9) 
through  VAR(12)  on  the  second.  The  following  table  gives  the 
notation  used  to  identify  VAR(l)  through  VAR(12)  on  the  trajec¬ 
tory  program  output  which  will  be  discussed  in  Section  4.3. 
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Program  Notation 


VAR 

(1) 

VAR 

(2) 

VAR 

(3) 

VAR 

(4) 

VAR 

(5) 

VAR 

(6) 

VAR 

(7) 

VAR 

(8) 

VAR 

(9) 

VAR 

(10) 

VAR 

(11) 

VAR 

(12) 

Output  Notation 
DXF ,  ft/sec 
DYF,  ft/sec 
DZF,  ft/sec 

P,  radians/sec 

Q,  radians/sec 

R,  radians/sec 
XF  of  XMOM,  ft 
YF  of  XMOM,  ft 
ZF  of  XMOM,  ft 
PSI,  degrees 
THETA,  degrees 
PHI,  degrees 


4.2.3  Sample  input  data 

A  sample  input  data  deck  to  the  second  program  will  now  be 
presented.  It  will  utilize  the  wind  tunnel  model  components  shown 
in  Figure  9  except  that  they  have  been  scaled  up  by  a  factor  of 
twenty  in  order  to  approximate  a  full-scale  aircraft  and  store. 

The  data  will  either  be  input  directly  in  the  data  deck  or 
indirectly  through  the  calculations  performed  in  Program  I  and 
passed  through  the  intermediate  dataset.  The  models  were  used  in 
the  wind  tunnel  test  program  described  in  Reference  11. 

The  configuration  for  the  sample  calculation  is  shown  in 
Figure  10  as  input  to  the  first  program.  All  geometric  data  for 
that  configuration  is  input  to  the  first  program  except  that 
pertaining  to  the  layout  of  the  planar  and  interdig itated 
fins.  The  trajectory  calculations  start  with  the  store  in 
the  carriage  position  as  shown. 

The  input  data  deck  for  the  second  program  is  tabulated  in 
Figure  19.  The  item  numbers  indicated  on  the  figure  correspond  to 
those  of  Figure  16.  The  first  item  of  Figure  19  is  item  number  1 


ll'b 


which  contains  the  value  of  HOARDS,  in  this  case  14.  This  is 
followed  by  item  number  2,  which  consists  of  these  14  cards  of 
identifying  information.  It  is  noted  that  the  tabulated  data 
presented  in  this  sample  case  and  the  resulting  calculations  are 
intended  to  simulate  the  experimental  captive  store  trajectory 
of  Test  V6A,  Group  363  in  Reference  11. 

Item  number  3  specifies  the  aircraft  flight  conditions  used  in 
the  experimental  case.  The  angle  of  attack  is  5.0°,  the  flight 
path  angle  is  0.0°,  and  the  Mach  number  is  1.5.  The  free-stream 
air  density  is  0.0005852  slugs  per  cubic  foot,  which  corresponds 
to  a  flight  altitude  of  40,000  feet.  The  flight  velocity  is 
1,462  feet  per  second. 

Item  4  contains  the  ejected  store  indices.  Store  number  10, 
NEJECT=10,  is  the  store  being  separated.  Since  the  store  is 
modeled  using  the  elliptic  body  option  NSEG  and  NSORCE  are  set 
equal  to  zero.  The  parameter,  NGAM=1,  is  set  to  specify  that  a 
captive-store  trajectory  is  to  be  calculated.  No  polynomials, 
NPOLY=0,  will  be  read  since  it  is  an  elliptic  configuration. 

Though  the  store  is  mounted  on  the  centerline,  the  rolling 
moment  is  computed,  NROLL=l,  to  verify  the  symmetry  of  the  load 
calculations.  The  store  has  two  empennages,  NEMP=2,  corresponding 
to  the  two  planar  fins  and  the  four  fins  in  an  interdigitated 
configuration.  No  aerodynamic  damping  will  be  computed,  NDAMP=0, 
and  no  thrust  time  history  is  to  be  included,  NTHRUS=0.  No  theta 
locations  are  required,  NTHETA=0,  for  the  elliptic  shape.  Ejector 
force  input  calculations  will  be  included,  NJECTR=1. 

Since  an  ejector  is  present,  the  next  four  items  in  Figure  19 
contain  the  descriptions  of  the  ejector  force  versus  time  history. 
Item  number  5  specifies  that  one  ejector  foot  is  present,  NFEET=1 , 
and  that  it  is  input  as  a  function  of  time,  NSTRKE=0.  The  next 
card,  item  number  6,  specifies  that  one  polynomial  will  be  used  to 
describe  the  ejector  time  history,  NEPOLY=l,  that  it  is  located 
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one  foot  aft  of  the  store  moment  center,  XE=1.0,  and  that  it 
is  pointing  straight  down,  THETAE=0°.  The  stroke  over  which 
the  force  occurs  is  set  to  a  large  number,  STROKE=10.0  feet,  to 
allow  the  force  to  be  terminated  when  exceeding  the  final  time 
specified  for  the  polynomial  time  history,  TEEND=0.075  seconds, 
in  item  number  7.  Item  number  8  contains  the  5th  order  poly¬ 
nomial  coefficients  for  the  ejector  force  time  history  from  the 
following  equation  from  Reference  11: 

F  =  1.41448  x 105T  +  2 . 68614  x 106T2  -  2.2209  * 108T3 
+  3.52968  x  109T4  -  1 . 70086  x  1010T5 

The  next  two  items  contain  the  mass  and  inertia  character¬ 
istics  of  the  separating  store.  Item  number  9  contains  the 

store  mass,  SMASS=20.0  slugs,  moments  of  inertia,  I  =20., 

2 

I  =250. ,  and  I  =250.  slug-ft  ,  and  products  of  inertia, 

yy  zz  ^ 

I  =0.0,  I  =0.0,  and  I  =0.0  slug-ft  .  Item  number  10  sets 
yz  xz  xy  ’ 

the  moment  center  (XMOM)  at  5.0  feet  and  specifies  that  there 
is  no  center  of  gravity  offset  from  the  moment  center. 

The  next  card  is  item  number  13.  It  contains  the  axial 
force  coefficient,  CA=0.73  obtained  from  the  experimental  data. 

The  next  fifteen  cards  are  used  to  describe  the  two  sets  of 
elliptic  store  empennages.  Item  number  16  is  used  to  partition 
the  elliptic  store  body  into  body-alone  and  body-fin  sections. 

The  first  value,  NBOD  =  4,  indicates  that  the  body  is  to  be  split 
into  four  sections  for  load  calculations.  The  variables  XBOD(I) 
and  LFIN(I)  indicate  the  length  over  which  the  section  extends  and 
whether  it  is  finned  body  section.  The  first  pair,  XBOD (1 ) =5. 833 
and  LFIN(l)=F{false) ,  say  that  the  unfinned  body  nose  extends  from 
the  nose  at  X=0.0  aft  to  5.833  feet.  The  second  pair,  XBOD(2)«8. 333 
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and  LFIN{2) =T (true)  ,  indicate  that  a  finned  body  section 
extends  from  5.833  feet  to  8.333  feet.  The  next  unfinned  section 
extends  from  the  forward  fin  trailing  edge  at  8.333  to  8.75  feet. 
The  last  section  is  finned  and  extends  from  8.75  to  10.0  feet. 

Two  sets  of  item  numbers  17  through  24  follow.  The  first  is 
for  the  planar  fins.  It  contains  the  input  parameters  to  describe 
the  fin  and  interference  shell  layo  it.  Item  number  17  is  the 
descriptive  title.  For  descriptive  purposes  here  the  1st 
empennage  has  been  designated  the  "monoplane  wing."  Item  number 
18  contains  the  control  indices  for  the  type  and  number  of  panels 
on  each  fin.  The  first  two  values,  MSWR  =  MSWL  =  4,  place  four 
rows  of  panels  in  the  spanwise  direction  on  each  of  the  horizontal 
fins.  The  next  two  values  are  zero  indicating  no  third  or  fourth 
fins  present.  The  number  of  constant  u-velocity  panels  in  each 
chordwise  row,  NCW  =  5,  and  the  number  of  panels  axially  on  the 
interference  shell,  NCWB=5,  are  set  equal  to  the  number  of  rings 
of  body  source  panels  used  to  span  that  region.  NBDCR  is  set 
to  zero  to  allow  the  program  to  use  the  default  values  correspond¬ 
ing  to  the  KRADX  for  the  body  panels.  The  last  value,  LVSWP=0, 
indicates  that  there  are  no  planform  breaks  in  sweep  and  the 
parameters  in  item  19  will  be  used  to  describe  the  fins.  Item 
number  19,  then,  defines  the  root  chord  (CRP=2.5  feet),  the 
leading  and  trailing  edge  sweeps  (SWLEP=53 . 13°  and  SWTEP=0°) , 
and  the  exposed  semispan  (B2=1.25  feet)  for  the  monoplane  wings. 
Item  number  20  is  not  included  since  there  are  only  two  fins. 

Item  number  21  states  there  are  no  interdigitated  angles  and 
that  there  is  no  deflection  of  any  of  the  fins.  If  the  polar 
angle  of  the  interdigitated  tail  is  zero  (THETIT  *  0°) ,  the 
individual  angles  for  each  of  the  wings  are  read  in  item  number 
22.  PHIFR  -  THETR  *  0°  indicate  the  right  fin  is  horizontal. 

PHIFL  *  THETL  =  0°  similarly  indicate  that  the  second  fin  is 
horizontal  but  to  the  left.  The  interference  shell  body  length, 
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BIL=2.5  feet,  in  item  number  23  is  the  same  as  the  root  chord. 

The  remaining  value  VRTMAX=0.35  is  the  recommended  default. 

NVRTIN  in  item  24  is  set  to  zero  for  the  forward  set  of  fins 
since  no  vortices  have  been  computed. 

Items  numbers  17  through  24  are  repeated  as  before  with  the 
geometry  for  the  interdigitated  tails  given  in  Figure  9(e).  In 
item  number  21  the  interdigitated  angles  (PHIDIH=THETIT=35°) 
are  used  to  lay  out  the  orientation  of  all  four  fins  in  lieu  of 
inputting  each  angle  in  item  22.  For  the  purpose  of  these  cal¬ 
culations  NVRTIN  is  set  to  1  in  item  24  so  that  the  vortices 
from  the  forward  empennage  do  not  affect  the  aft  set  of  fins. 

Item  number  32  specifies  initial  axial  and  lateral  transla¬ 
tional  velocities  of  the  store  of  0  feet  per  second  in  all 
directions.  This  corresponds  to  the  initial  condition  of  the 
store  fixed  at  the  carriage  position.  Similarly,  the  initial 
rotational  velocities  are  all  set  to  0.0  radians  per  second. 

The  last  card  input  in  Figure  19,  item  33,  provides  the 
integration  interval  and  the  initial  and  final  times.  Since 
the  initial  time  is  0.0,  item  34  is  not  included. 

4.3  Description  of  Output  from  Program  II 

Figure  20  presents  the  output  from  Program  II  for  the  sample 
case,  the  separation  of  the  elliptic  store  from  the  configuration 
in  Figure  10. 

The  first  page  of  output.  Figure  20(a),  tabulates  the  input 
identifying  information  and  the  parent  aircraft  atmospheric  flight 
conditions  of  item  3.  The  first  eight  title  lines  contain  the 
heading  information  read  in  Program  I  and  transferred  through  the 
external  file.  The  next  14  lines  of  title  information  contain 
the  additional  heading  information  read  in  Program  II  which 


127 


contain  information  specific  to  the  present  run.  The  aircraft 
flight  conditions  follow. 

The  information  on  the  next  seven  pages  contains  the 
remaining  information  describing  the  separating  store,  the 
ejector  force,  and  the  resulting  layout  of  the  two  empennages. 

In  Figure  20(b),  the  ejected  store  is  identified  and  the  options 
in  input  item  4  specified.  They  are  followed  by  the  ejector 
force  history  input  and  polynomials.  The  additional  input 
specifying  the  store  moments  of  inertia,  mass,  center  of  gravity 
location,  and  axial  force  coefficient  are  presented  next. 

Following  this,  the  finned  body  layout  output  identifies  the 
store  body  stations  which  partition  the  source  panels  into  body 
alone  or  empennage  sections  for  force  calculation  purposes.  The 
first  two  columns  repeat  the  input  of  item  16.  The  last  column 
identifies  the  x-station  index  of  input  item  75  of  Program  I 
closest  to  the  input  x-station  of  the  first  column  which  will 
actually  be  used  to  partition  the  body  into  sections.  This  is 
followed  by  the  empennage  leading-edge  locations  of  the  elliptic 
store  fins. 

The  third  page  of  output,  Figure  20(c),  contains  the  output 
of  the  input  data  describing  the  first  empennage  on  the  elliptic 
store,  designated  the  monoplane  wing.  Fin  properties  are  input 
for  only  two  fins.  Of  the  body  properties,  only  the  body  inter¬ 
ference  shell  length  and  its  number  of  streamwise  body  panels  are 
input  at  this  point  in  the  program.  The  remaining  quantities, 
the  elliptic  axes  and  the  number  of  panels  around  the  body  were 
obtained  from  the  source  panel  layout  from  Program  I.  The 
interdigitated  angles  are  not  used  for  this  layout.  The  reference 
area  and  length  are  the  reference  area  and  diameter  Refined 
for  the  elliptic  store  body  in  the  first  program.  The  vorticity 
options  are  set  to  allow  the  program  to  compute  the  trailing-edge 
vorticity  off  the  monoplane  wing  for  integration  aft  to  the  second 
set  of  surfaces. 
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The  next  two  pages,  Figures  20(d)  and  20(e)  contain  the 
constant  u-velocity  panel  corner  coordinates  for  the  two  fins  and 
the  interference  shell,  respectively.  These  are  provided  for  use 
in  verifying  the  layout  of  the  panel  geometry.  The  coordinates 
are  in  the  empennage  coordinate  system,  xWl ,  yw^,  z^ ,  shown  for 
the  first  store  empennage  in  Figure  13. 

The  next  three  pages,  Figures  20(f),  20(g),  and  20(h)  repeat 

the  input  for  the  second  empennage  on  the  elliptic  store, 

designated  the  interdigitated  tail.  For  this  empennage  all  four 

fins  are  input  and  the  interdigitated  tail  angles  are  used  to 

specify  their  orientation.  The  panel  coordinates  in  the 

empennage  coordinate  system,  x  ,  yw  ,  z  ,  shown  in  Figure  13 

" 2  w2  w2 

follow. 

The  last  three  parts  of  Figure  20,  parts  (i) ,  ( j) ,  and  (k) , 
are  output  for  three  points  in  the  trajectory,  the  first  point, 
an  intermediate  point,  and  the  point  at  which  the  trajectory  is 
terminated.  The  complete  output  from  the  program  will  contain  a 
page  like  this  for  each  integration  step,  in  this  case  every 
0.025  seconds. 

At  the  top  of  each  page  is  the  trajectory  time  in  seconds. 
Following  this,  the  forces  and  moments,  components  as  well  as 
totals  are  printed.  These  are  followed  by  the  ejector  produced 
forces  and  moments.  Next,  the  axial  distributions  of  body-alone 
normal  force  and  side  force  are  tabulated.  These  were  integrated 
to  determine  the  "BODY"  CN  and  CY  listed  at  the  top  of  the  page. 
The  force  and  moment  calculation  is  discussed  in  detail  in 
Section  6  of  Reference  1. 

The  following  table  relates  the  program  output  variables  to 

the  x  ,  y  ,  z  coordinate  system  and  positive  directions  shown 
s  s  s 

in  Figure  17. 
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Program  Notation 


Notation  of  Figure  17 


CN 

CN 

CY 

CY 

CLM 

Cm 

CLN 

C 

n 

CLL 

X,  FT 

xs,  feet 

X/L 

xs/£s 

DCN/DX 

dCN/dxs , 

per 

DCY/DX 

dC y/ dxs , 

per 

foot 

foot 


As  the  store  pitches,  yaws,  and  rolls  during  the  trajectory,  the 
x _ ,  y„,  z_  coordinate  system  pitches,  yaws,  and  rolls  with  it.  The 

S  S  S 

forces  are  always  calculated  in  this  coordinate  system. 


The  remaining  quantities  tabulated  on  each  page  of  the  trajectory 
output  specify  the  store  location,  orientation,  velocities,  and 
accelerations  relative  to  the  parent  aircraft  at  that  particular 
time.  Before  discussing  these  quantities,  the  coordinate  systems 
must  be  mentioned.  Figure  17  shows  another  coordinate  system,  x, 
y,  and  z,  which  is  fixed  in  the  store  and  moves  with  the  store  as 
it  yaws,  pitches,  and  rolls.  The  origin  of  this  system  is  fixed 
at  the  store  moment  center.  This  coordinate  system  is  also  shown 
in  Figure  21  along  with  another  system,  £,  ri ,  5.  This  latter 
system  is  an  inertial  coordinate  system  whose  origin  is  fixed  in 
the  fuselage  nose  and  is  parallel  to  the  x_,  y_,  zQ  system  of 

ODD 

Figure  5.  At  any  given  time,  the  two  coordinate  systems  are 
orientated  with  respect  to  each  other  by  a  system  of  angles.  The 
angles  are  those  shown  in  Figure  21  and  consist  of  three  rotations 
in  the  yaw,  y,  pitch,  0,  and  roll,  $,  sequence.  The  positive 
senses  of  the  three  store  rotational  velocities  about  the  x,y,z 
axes  are  also  shown  in  the  figure. 
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Following  the  body  load  distribution  output  on,  for 
example.  Figure  20{i),  the  location  of  the  store  in  the  fuselage, 
or  inertial,  coordinate  system  is  tabulated.  The  locations  of 
the  store  nose,  NOSE,  moment  center,  XMOM,  and  base,  BASE,  are 
tabulated  relative  to  the  fuselage  nose  and  also  relative  to  the 
position  of  the  store  at  time  t  =  0,  the  initial  position.  In 
this  tabulation  XF  is  xg  of  Figure  5  or  £  of  Figure  21.  Likewise 
YF  is  yB  or  n  and  ZF  is  zB  or  z, . 

The  next  output  are  the  translational  velocities  and  accelera 
tions  of  the  store  relative  to  the  moving  aircraft.  For  example, 

DXF  -  dT  at  '  ft/sec 

^  XR  2 

D2XF  =  - £  or  — 2-  ,  ft/sec 

dt^  dt^ 

Next  the  rotational  velocities  shown  in  Figure  21  are  listed  as 
are  the  rotational  accelerations.  The  notation  is: 

P  =  p,  radians/sec 

Q  =  q,  radians/sec 

R  =  r,  radians/sec 

j  _  2 

PDOT  =  ,  radians/sec 

QDOT  =  ^  f  radians/sec2 

j  „  2 

RDOT  *  ,  radians/sec 
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The  last  output  printed  at  each  integration  step  are  the  values 
of  the  three  orientation  angles  shown  in  Figure  21  and  their  time 
rates  of  change.  The  notation  is 

PSI  ■  4*,  deg. 

THETA  =  0,  deg. 

PHI  =  4>,  deg. 
d¥ 

DPSI  =  ,  radians/sec 

DTHETA  =  ,  radians/sec 

cl  $ 

DPHI  =  ^  ,  radians/sec 


4.4  Program  Error  Messages 

The  possible  occurrence  of  certain  fatal  errors  and  irresolvable 
conflicts  of  data  during  program  execution  have  been  anticipated 
in  Program  II.  If  such  errors  should  occur,  the  program  has  been 
designed  to  print  a  diagnostic  error  message  and  halt  when 
appropriate.  Each  message  with  the  routine  in  which  the  cor¬ 
responding  error  is  detected  is  presented  in  the  following  table. 

The  program  test  which  resulted  in  the  action  is  aijo  included. 

The  table  is  arranged  in  alphabetical  order  by  the  program  or 
subroutine  from  which  it  was  written. 


TRJTRY  STOP  712  Angle  of  attack  and/or  Mach  number  read 

in  does  not  agree  with  value  from  file. 
(Test:  ALFAR  /  ALFACR  or  CMACH  FMACH) 
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TRJTRY  STOP  normal  termination 

(test:  TIME  >  TIMEF-0. 00001  or 
NDIFEQ  =  1) 

BDCOEF  message  warning  -  NAFLD  exceeds  blank  common, 

**BDC OEF** ,  NTAP7= _ ,  IUB* _ ,  NAFLD= _ . 

(Test:  NAFLD  >  NADIM) 

CRFWBD  STOP  001  (Test:  NER  >  1;  equation  solution 

singular) 

INVER2  STOP  **Matr ix  is  singular,  INVER2 

(Test:  A(I,I)=0.0) 

NUMACH  STOP  **Axial  location  of  traverse  above  wing 

chordal  plane,  **  NUMACH  **. 

(Test:  Z  >  ZLOC) 

NUMACH  STOP  101  (Test:  X+BETA*Z+2*CHRD-WLEX+BETA*Z  *  0; 

search  for  maximum  turning  angle  failing) 

SORPAN  STOP  230  **  Mach  number  less  than  1,  SORPAN 

(Test:  AMACH  <  1.) 

STOP  220  Error  -  body  panel  slope  exceeds  Mach 

angle,  J= _ ,  **  SORPAN  ** 

(Test:  1.-BETAL2*TAND2<0.) 

If  one  of  the  above  messages  is  printed  out  during  program 
execution,  the  user  should  consult  the  input  data  first  for 
possible  error.  The  program  variable  descriptions  in  Appendices 
C  and  D  may  also  be  helpful  for  the  meanings  of  the  above  test 
quantities. 
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4.5  Program  II  Running  Times 


The  second  program  described  in  this  report  has  been  run  on 
both  the  CDC  6600  and  CYBER  175  machines.  Because  of  machine  differ¬ 
ences,  the  running  time  varies  from  one  machine  to  another.  As  a 
consequence,  only  an  approximate  running  time  can  be  given.  The 
second  program  can  be  executed  multiple  times  for  the  single 
execution  of  the  first  program. 

The  running  time  for  the  calculation  of  the  trajectory  of  the 
separated  store  is  a  function  of  a  number  of  factors,  some  of 
which  are: 

(a)  The  number  and  method  used  to  model  the  fixed  stores, 
rack  and  fuselage  bodies. 

(b)  Number  of  sources  and  doublets  required  to  represent 
circular  fuselage,  rack,  or  store  bodies. 

(c)  Number  of  source  panels  used  to  represent  the  noncircular 
fuselage  or  elliptic  stores. 

(d)  Number  of  constant  u-velocity  panels  on  the  wing,  pylon, 
and  elliptic  store  empennages. 

(e)  Number  of  thickness  panels  on  the  wing  and  pylon. 

(f)  Number  of  body  segments  and  tail  fin  control  points  used 
in  force  and  moment  calculations  for  the  circular  store. 

(g)  Duration  interval  of  separated  store  in  presence  of 
reflected  shocks. 
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(h)  Integration  interval. 

(i)  Real  time  duration  of  trajectory. 

All  of  the  above  should  be  kept  to  the  minimum  required  to  resolve 
the  separated  store  loads.  For  the  sample  case  for  the  second 
program  in  Figure  19,  the  run  time  was  about  1,310  second  on  the 
CYBER  175. 


# 
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(b) 

Figure  1.-  Concluded 


!pu*  PROGRAM  I  -  L.I'CAl, 


Figure  2.-  Program  I  input  formats. 


Continued 


Figure  2.-  Continued 


Figure  2.-  Continued 


ZBWO 


(a)  Definition  of  certain  variables 
Figure  7.-  Wing  input  variables. 


(b)  Example  wing, 
igure  7.-  Concluded. 


(a)  Pylon  under  fuselage  centerline. 


(b)  Pylon  under  wing. 


Figure  8.-  Variables  describing  and  locating 
pylon,  input  data  item  number  43. 
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Figure  9.-  Wind-tunnel  models. 


Continued 


I 


All  dimensions  in  inches 


(d)  Pylon,  (P2) c 
Figure  9.-  Continued. 
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Figure  11.-  Input  data  deck  for  sample  case  for  Program  I. 
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Figure  11.-  Continued. 


165 


Open  inlet  panels 
Blocked  inlet  panels 


★  Additional  panels  used 
in  inlet  velocity 
calculation 


Fuselage  frontal  view  looking  upstream 
showing  inlet  panels 

Figure  12.-  Continued. 


168 


9$  ^  *101,102 

* -Additional  panels  used  in 
inlet  velocity  calculations 


Inlet  layout  side  view  showing  panels  used 
in  inlet  velocity  calculations 


Figure  12.-  Concluded. 
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Figure  14.-  Program  I  output  for  sample  case 
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Figure  14.-  Continued. 
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Figure  14.-  Continued. 
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Figure  14.-  Continued. 
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Figure  14.-  Continued 
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Figure  14.-  Continued. 
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Figure  14. -  Continued. 
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Figure  14.-  Continued. 
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Figure  14.-  Continued. 
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Figure  14.-  Continued. 
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Figure  14.-  Continued. 
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Figure  14.-  Continued. 
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Figure  14.-  Continued. 
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Figure  14.-  Continued. 
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Figure  14. -  Continued. 
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Figure  14.-  Continued. 
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Figure  14.-  Concluded. 
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Input  PROGRAM  2  -  TRJTRY 


Figure  16.-  Program  II  input  formats 
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Figure  16.-  Continued 
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Figure  16.-  Concluded 


219 


-**  *  V»  -  *K, i-l 
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these  functions  are  used  in 
the  transformation  (rotation 
cf  the  body  interference 
panel  to  the  reference 

< *B  -  >'  b-  *B  1  or  (XV’ 
syster* 


;F  •  fir.  dihedral  angle,  PH1D1H 

*  •  dihedral  angle  of  right  upper 
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r,  .  ■  polar  angle  of  left  lower 
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Body  interference  panels : 

W.j  *  THTI  -  • 


*-,,p  V“T 


Figure  18.-  Geometrical  angles  associated  with  case  involving 
interdigitated  fins  on  body  with  elliptical  cross  section 
and  force  coefficients  associated  with  fins  and  body 
interference  panel  (view  looking  upstream) . 
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Figure  19.-  Program  II  sample  input  case 
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Figure  20.-  Program  II  output  for  sample  case. 
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Figure  20.-  Continued. 
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Figure  20.-  Continued. 
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Figure  20.-  Continued 
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